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We  have  made  substantial  progress  In  understanding  the  behavior 
of  multiplled-Polsson  processes,  and  In  particular,  an  Important  special 
case  of  such  processes,  the  shot-nolse-driven  doubly  stochastic  Poisson 
point  process  (SNDP) .  The  SNDP  turns  out  to  provide  an  appropriate  model 
for  Important  effects  Inherent  in  the  generation  and  detection  of  optical 
and  ionizing  radiation.  Including: 

(1)  the  detection  of  scintillation  photons  produced  by  ionizing 
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radiation; 

(11)  the  generation  of  luminescence  noise  by  ionizing  radiation*  -and 

(ill)  cathodo luminescence  emission. 

The  SNDP  can  be  quite  important  in  properly  interpreting  the  behavior  of 
optical  detectors,  fibers  and  electro-optic  devices  and  in  determining  the 
performance  of  communications,  surveillance,  navigation,  guidance  and  control 
systems  incorporating  these  devices.  The  SNDP  will  also  be  useful  for  studying 
image  intensification,  cosmic-ray  noise  in  electro-optic  systems,  and  cosmic 
ray  soft  errors  in  electronics  and  computer  systems . 

Substantial  progress  has  also  been  made  toward  the  analysis  and  implemen¬ 
tation  of  the  first  wideband  Incoherent  frequency-shift  keyed  (FSK)  fiber¬ 
optic  coonunicatlon  link.  Results  were  obtained  in  three  specific  areas: 

(1)  optimal  (single-threshold)  processing  and  receiver  design;  (ii)  new 
technique  for  system  performance  analysis;  (ill)  FSK  system  demonstration. 

In  the  area  of  receiver  design,  a  simplified  (single-threshold)  processor 
was  derived  that  is  optimal  for  a  broad  variety  of  fiber-optic  systems 
using  PIN  or  avalanche  photodiodes.  A  new  technique  was  developed  for  cal¬ 
culating  fiber-optic  system  performance  using  a  normalizing  transform.  Whereas 
previously,  Gaussian  approximations  were  resorted  to  for  such  calculations, 
the  new  technique  is  as  simple,  but  more  accurate.  Finally,  a  10  Mbit/sec 
wideband  FSK  fiber-optic  link  was  implemented  with  center  wavelength  827  nm, 
using  a  direct  detection  avalanche  photodiode  receiver.  The  performance 
was  shown  to  be  3-4. 5  dB  better  than  a  comparable  on-off  keyed  system.  It 
is  expected  that  the  FSK  fiber-optic  link,  when  combined  with  wavelength 
multiplexing,  will  be  important  in  high  -data-rate,  high-performance  digital 
communications . 

We  have  developed  a  general  theoretical  model  of  grain  boundary 
(GB)  recombination,  carrier  transport  and  electrostatics  under  assump¬ 
tions  of  Gaussian  energy  distributions  of  GB  interface  states,  and 
unequal  capture  cross-section  of  these  GB  states  for  electrons  and  holes 
in  polyctystalline  materials.  Calculations  have  been  performed  of  the 
recombination  current  density  and  the  recombination  velocity  at  grain 
boundaries  in  polycrystalline  silicon  for  four  different  energy  dis¬ 
tributions.  It  appears  that  the  experimental  data  of  photoconductivity 
variation  with  light  intensity  are  consistent  with  the  multi-level  trap 
distribution  proposed  although  detailed  correlation  requires  further 
studies . 

We  report  here  data  of  the  capacitance-voltage  and  current- 
voltage  measurements  on  cw  laser  processed  Al-nSi  diodes.  The  range  of 
the  laser  power  is  between  12  and  20W.  Significant  changes  in  the 
ideality  factor,  reverse  saturation  current,  thermal  activation  energy, 
doping  concentration,  and  effective  Schottky  barrier  height  have  been 
observed.  These  results  appear  to  be  consistent  with  the  mechanisms 
of  aluminum  diffusion  and  defect  generation  with  high-temperature 
stress  under  high-power  laser  irradiation. 

The  grain  boundary  recombination  velocity  of  polycrystalline  Si  Is 
separated  from  bulk  effects  and  its  value  is  determined  with  practically  no 
assumption  of  numerical  values  for  any  parameter  other  than  the  minority 
barrier  diffusion  coefficient.  The  method  used  utilizes  photoconductance 
measurements  taken  by  the  scanning  laser  spot  technique.  It  is  found  that. 
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Block  20  continued  -  Abstract 

£or  the  illumination  levels  used,  there  is  a  transition  from  low  level  to 
high-level  injection  of  minority  carriers,  and  this  transition  is  used  to 
estimate  the  quasi-Fermi  level  separation.  The  recombination  velocity 
S,  ranging  in  value  from  5000  to  17000  cm/sec.  is  thus  presented  both  as 
a  function  of  the  illumination  level  as  well  as  of  the  quasi-Fermi  level 
separation.  It  is  found  that  S  increases  monotonically  with  illumination 
and  does  not  saturate,  even  at  the  high  light  concentrations  used. 

A  physical  theory  using  a  charge  scattering  model  is  proposed  to 
interpret  the  experimental  data  of  grain  boundary  transport  in  poly¬ 
crystalline  semiconductors.  The  calculated  result  explains  the  need 
of  an  attenuation  factor  as  an  added  coefficient  in  the  thermionic  emission 
current . 

The  behavior  of  Al-nSi  diodes  heat  treatment  above  the  eutectic 
temperature  (600°C)  is  studied  by  varying  the  annealing  time  and  cooling 
rate.  We  have  used  the  capacitance-voltage  characteristic  as  a  measure 
of  the  influence  of  the  doping  density  and  thickness  of  the  recrystallized 
silicon  layer  at  the  interface.  Poisson's  equation  is  solved  to  obtain 
the  C-V  characteristic  and  to  relate  it  to  the  effective  barrier  height. 

We  find  that  the  barrier  height  is  independent  of  the  annealing  time,  but 
it  is  a  function  of  the  cooling  rate.  It  appears  that  the  diffusion  of 
alluminum  into  silicon  does  not  play  a  significant  role,  but  fast  quenching 
tends  to  prevent  aluminum  from  precipitation,  so  that  the  recrystallized 
silicon  layer  is  highly  doped. 

A  new  laboratory  for  studying  novel  processing  and  diagnostic 
techniques  for  microelectronics  materials  has  been  established.  One  major 
research  area  of  this  laboratory  is  the  application  of  laser  surface 
chemistry  to  microfabrication.  In  conjunction  with  this  experimental 
research,  a  theoretical  study  has  shown  the  importance  of  surface  composi¬ 
tion  and  morphology  in  controlling  laser-initiated  photochemistry. 

We  have  determined  the  relaxation  dynamics  of  electronically  excited 
bromine  molecules  imbedded  in  rare  gas  matrices  of  Ar,  Kr,  and  Xe  at  temper¬ 
atures  ranging  from  10  -  40°K.  A  restriced  relaxation  mechanism  leads 
to  emission  for  all  electronic  states  (B,  A,  and  A')  energetically  accessible 
to  the  laser  photons  used  for  excitation.  Preliminary  results  for  iodine 
molecules  indicate  that  emission  is  only  observed  from  the  A  electronically 
excited  state.  The  relaxation  mechanism  for  iodine  will  provide  information 
useful  in  determining  the  upper  limits  of  efficiency  for  chemical  iodine 
lasers . 

The  vibrational  relaxation  of  oxalyl  fluoride  has  been  studied  and 
found  to  be  very  rapid.  The  presence  of  a  low-lying  torsional  frequency 
in  this  molecule  enhances  the  relaxation  efficiency  by  orders  of  magnitude 
over  rigid  molecules  of  roughly  the  same  size.  Vibrational  relaxation  of 
0CS  in  He  has  also  been  investigated  to  determine  if  an  unusual  cyclic 
relaxation  path  occurs  in  this  molecule.  Definitive  evidence  for  this 
process  is  still  lacking.  — * 
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Block  20  continued  -  Abstract 

Photofragmentation  of  HBr  in  the  presence  of  CO2  and  CO  has  been 
investigated.  Vibrational  excitation  of  CO2  and  CO  by  fast  H  atoms  produced 
in  the  photodissociation  of  the  diatomic  has  been  observed.  In  CO2  collisions 
with  fast  hydrogen  atoms  produce  predominantly  bend/stretch  excitation  while 
in  CO  very  high  levels  of  rotational  excitation  are  observed  (AJ  of  approximately 
20).  Photofragmentation  of  triazine  (HCN) 3  has  also  been  studied.  This 
molecule  cleanly  produces  three  HCN  molecules  when  exposed  to  193  nm  excimer 
laser  light.  We  have  observed  strong  vibrational  excitation  of  both  the 
bend  (V2)  and  C-H  stretch  (V3)  modes  of  HCN  produced  via  photodissociation. 
Population  inversions  may  exist  for  these  molecules. 

A  new  capillary  ion  source  has  been  constructed  and  coupled  to  a 
standard  molecular  beam/mass  spectrometer  system.  Preliminary  results  with 
this  device  show  that  hot  ions  having  a  significant  spread  in  ion  energy 
are  produced  at  discharge  voltages  of  a  few  kilovolts.  Emission  spectra  of 
the  discharge  indicate  that  substantial  quantities  of  metastable  neutrals 
are  produced  by  this  source  in  addition  to  fast  ions. 

The  origin  and  dynamics  of  an  anomalous  dual  fluorescence  in  molecules 
such  as  p-dimethylamino  benzonitrile  have  been  investigated  using  picosecond 
laser  sources.  The  key  issues  of  the  relationship  between  the  ultraviolet 
and  visible  emissions  and  the  role  of  the  solvent  in  this  photolnduced  charge 
separation  in  molecules  are  elucidated. 

We  have  found  in  our  studies  of  small  chain  motions  in  molecules  of 
the  type  anthracene-(CH2) 3-N ,N ’-dimethyl  andine,  using  picosecond  laser  methods 
that:  (1)  the  dynamics  of  end  to  end  reactive  encounters  are  influenced  by 
different  starting  ground  state  conf ormations ,  (2)  the  lifetime  of  chain  relaxa¬ 
tion  scales  as  (viscosity) 1/2  and  (3)  that  a  single  chain  relaxation  mode, 
perhaps  the  large  scale  end  to  end  relaxation  mode,  is  dominant  for  the  three 
methylene  chain  in  the  picosecond  to  nanosecond  time  domain. 

We  have  used  picosecond  spectroscopy  to  study  the  pathway  involved 
in  the  generation  of  singlet  oxygen,  from  photoexcited  anthracene 

endoperoxide.  We  have  determined  the  dynamics  and  relative  importance  of 
the  various  channels  for  photofragmentation  of  these  molecules. 

A  centralized  laser  pump  facility  was  developed  to  allow  a  pair  of 
YAG  lasers  to  be  simulataneously  shared  at  three  separate  experimental 
stations.  The  Billiard  Ball  Echo  Theory  was  invented  to  explain  most 
optical  coherent  transient  phenomena  in  an  unconventional  but  simple  way. 

New  optical  coherent  effects  were  discovered  which  will  be  used  to  obtain 
physical  data.  A  theoretical  analysis  of  the  combined  effects  of  phase 
changing  and  velocity  changing  collisions  in  gases  was  completed.  Measurements 
were  made  to  determine  the  relative  orientation  of  the  principle  axes  of 
the  ground  and  excited  state  nuclear  hyper fine  hamiltonians  of  Pr3+  in  LaF^. 
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I.  QUANTUM  DETECTION  AND  SENSING  OF  RADIATION 

A.  WAVE  AND  PARTICLE  NOISE  IN  OPTICAL  DETECTION  AND  COMMUNICATIONS* 

(M.  C.  Teich,  P.  R.  Prucnal) 

(Principal  Investigator:  M.  C.  Teich  (212)  280-3117) 

(JSEP  work  units  12,  13  and  14,  1979-1982) 

We  have  made  substantial  progress  in  understanding  the  behavior  of 
multiplied-Poisson  processes,  as  Indicated  below. 

For  a  Poisson  primary  process  generating  a  Poisson  number  of  secondaries 
for  each  primary,  we  arrive  at  an  important  special  case,  the  shot-noise-drlven 
doubly  stochastic  Poisson  point  process  (SNDP),  as  illustrated  in  Fig.  1. 

The  sum  of  the  functions  associated  with  the  primary  events  constitutes  a  filtered 
Poisson  point  process,  or  shot  noise  (Fig.  lb).  Since  the  primary  events  are 
stationary,  this  will  also  be  a  stationary  stochastic  process.  The  subsidiary 
events  now  form  a  doubly  stochastic  Poisson  point  process  whose  stochastic  rate 
is  the  shot  noise  X(t);  this  provides  the  rationale  for  the  name  SNDP.  The 
structure  of  the  model  is  schematically  illustrated  in  the  block  diagram  of 
Fig.  2,  which  clearly  shows  its  multiplicative  Poisson  nature.  Physical  systems 
that  are  describable  by  an  SNDP  model  will  exhibit  a  physical  structure  charac¬ 
terized  by  these  three  cascaded  blocks. 


SNDP 

f(t) 


Figure  2.  Model  for  generation  of  the  SNDP  process. 
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The  SNDP  model  turns  out  to  provide  an  appropriate  description  for  a 
variety  of  important  effects  involving  multiplication.  These  effects  are  often 
inherent  in  the  generation  and  detection  of  optical  and  ionizing  radiation.  This 
multiplication  noise  can  be  quite  important  in  properly  interpreting  the  be¬ 
havior  of  individual  electro-optic  devices  and  in  determining  the  performance 
of  communication,  surveillance,  navigation,  guidance,  and  control  systems  in¬ 
corporating  these  devices. 

We  have  carried  out  a  number  of  studies  to  ascertain  the  importance  of 
multiplication  noise. ^  ^  We  have  explicitly  demonstrated  the  applicability 
of  the  SNDP  model  for  describing  the  following  processes:  (i)  the  detection  of 
scintillation  photons  produced  by  ionizing  radiation,  e.g.,  high-energy  elec¬ 
trons  or  nuclear  particles;  (ii)  the  generation  of  luminescence  noise,  by 
ionizing  radiation,  in  optical  detectors,  optical  fibers,  and  electro-optic 
components;  and  (ill)  cathodoluminescence  emission.  We  pointed  out  that  the 
SNDP  model  will  also  be  useful  for  studying  other  important  effects,  such  as 
image  intensification,  cosmic-ray  noise  in  electro-optic  systems,  and  cosmic- 
ray  soft  errors  in  electronics  and  computer  systems.  In  all  of  these  cases, 
there  is  a  multiplication  process  involved,  in  which  a  primary  event  produces 
a  random  number  of  secondary  events,  with  a  certain  time  course.  The  SNDP  may 
be  useful  in  other  electronics  applications  as  well,  e.g.,  for  determining 
scaling  limitations  and  reliability  factors  in  submicron  structures  and  VLSI 
systems  resulting  from  statistical  effects.  If  dopant  impurities  cluster 
to  an  appreciable  degree,  as  has  been  reported, ^  this  can  be  accommodated 
in  the  framework  of  a  spatial  multiplicative  noise  description.  The  doubly 
stochastic  Thomas  distribution^  might  provide  a  first  approximation  to  such 
behavior,  when  there  is  such  clustering  and  when  there  is  a  simultaneous  un¬ 
certainty  in  the  impurity  profile.  We  now  proceed  to  discuss  multiplication 
noise  in  more  detail  for  the  three  processes  cited  above. 
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(i)  Scintillation  Photon  Counting 


The  detection  of  ionizing  radiation  is  often  accomplished  through  a 
radiation-matter  interaction  in  which  a  single  high-energy  particle  produces 
a  shower  of  lower  energy  particles.  A  case  in  point  is  the  scintillation 
detector  which  is  a  combination  of  a  scintillation  crystal  (e.g.,  NaI:Tl, 
plastic)  with  a  photomultiplier  tube.  Conditions  for  the  validity  of  the 
SNDP  in  describing  scintillation  detection  are  that  the  incident  primary  ion¬ 
izing  particles  (e.g.t  electrons,  gammas,  protons)  be  represented  as  a  homoge¬ 
neous  Poisson  point  process,  and  that  each  primary  event  have  associated  with 
it  an  impulse  response  function  h(t)  that  governs  the  rate  of  production  of 
Poisson  optical  photons. 

A  model  for  the  physical  process  is  illustrated  in  Fig,  3.  It  is 
seen  to  be  Identical  in  form  to  the  block  diagram  presented  in  Fig.  2,  so 
that  the  resulting  photon  emissions  form  an  SNDP,  The  function  h(t)  will 
often  be  able  to  be  roughly  represented  by  a  decaying  exponential  so  that 
the  counting  and  time  statistics  can  be  calculated  explicitly.  The  descrip¬ 
tion  is  completely  characterized  by  the  driving  rate  p  and  the  impulse-response 
function  h(t),  and  therefore,  in  this  case,  by  three  parameters:  u,  the  rate 
of  the  primary  process;  a,  the  multiplication  parameter;  and  t^,  the  lifetime 
of  the  process  of  secondary  event  generation.  If  we  perform  photon  counting, 
we  must  add  T,  the  counting  time.  In  the  limit  of  counting  times  long  in 
comparison  with  the  exponential  decay  time,  the  counting  distribution  will 
reduce  to  the  Neyman  Type-A  for  arbitrary  h(t) .  It  is  usually  assumed  in  the 
literature,  generally  tacitly,  that  scintillation  photon  counting  statistics 
are  described  by  the  fixed  multiplicative  Poisson  distribution;  this  provides 
an  appropriate  approximation  only  for  T  >>  tc  and  a  >>  1,  however, ^ 


PHOTON 

POINT 

PROCESS 


Figure  3.  Model  for  the  photon  point  process  generated  by  scintillation 

radiation,  cathodoluminescence ,  and  photoluminescence.  Observe 
the  relation  to  Fig.  2  which  is  the  mathematical  model  studied 
here. 

We  have  demonstrated  that  the  experimental  counting  distribution  for 

radioluminsecence  photons  produced  in  glass  can  be  represented  in  terms  of 
(2) 

the  Neyman  Type-A.  That  this  is  plausible  can  be  understood  from  the 
following  argument.  Consider  a  stream  of  Poisson-distributed  3  rays  irradia¬ 
ting  a  glass  sample.  If  each  high-energy  electron  produces  a  Poisson  dis¬ 
tributed  number  of  luminescence  photons  with  a  maximum  time  delay  that  is 

short  in  comparison  with  the  counting  time  (t  <<  T) ,  the  overall  photon  count 

c 

will  reflect  both  sources  of  independent  Poisson  randomness,  leading  to  the 
Neyman  Type-A  photon-counting  distribution. 

(2)  - 

A  series  of  such  experiments  was  carried  out.  3  particles  from  a 

90  90 

Sr-  Y  equilibrium-mixture  source  irradiated  the  Corning  7056  glass  faceplate 

of  an  EMR  type  541N-01  photomultiplier  tube  from  a  distance  of  about  11.5  cm. 

-  90  90 

The  maximum  3  energies  wer\  0.54  and  2,23  MeV  for  the  Sr  and  Y,  respec- 

3-2-1 

tively,  and  the  3  flux  was  ^8,2  x  10  cm  s  .  External  light  was  excluded. 
The  photomultiplier  anode  pulses  were  passed  through  a  discriminator  and  stan¬ 
dardized.  Unavoidable  system  dead  time  was  ^60  ns.  The  standardized  pulses 
were  counted  during  consecutive  fixed  counting  intervals  (T  *  40  us)  and  the 
counts  were  recorded.  The  experiment  was  performed  repeatedly  to  obtain  good 
statistical  accuracy,  and  a  histogram  representing  the  relative  frequency  of 


Che  counts  was  conscrucced.  The  cocal  duration  of  a  run  was  about  4  min.  In 
the  particular  experiment  we  Illustrated,  the  observed  mean  count  was  85.89 
(this  number  was  substantially  higher  than  the  mean  dark  count  which  could 
therefore  be  neglected)  and  the  observed  count  variance  was  429.58.  The  data 
are  shown  as  the  dots  in  Fig.  4.  The  solid  curve  represents  the  Neyman  Type-A 
theoretical  counting  distribution  with  the  count  mean  and  variance  fixed  at 
the  experimental  values.  It  is  clearly  in  accord  with  the  data.  The  experi¬ 
mental  multiplication  parameter  a*  4.0.  A  Poisson  distribution  with  mean 
85.89  (indicated  by  arrow)  is  plotted  as  the  dashed  curve  in  Fig.  4;  clearly 
it  bears  no  relation  to  the  data. 

Because  the  primary  process  in  this  case  consisted  of  high-energy 
charged  particles  (electrons) ,  Cerenkov  radiation  could  have  been  produced 
in  addition  to  luminescence  radiation.  However,  even  if  a  large  number  of 
photoelectrons  were  generated  by  the  Cerenkov  photons  arising  from  a  single 
particle,  they  would  nevertheless  appear  as  a  single  (large)  photoelectron 
pulse  since  the  Cerenkov  radiation  emission  time  is  much  shorter  than  the 
transit  time  In  the  photomultiplier  tube.  The  appropriate  counting  distribu¬ 
tion  has  been  worked  out  for  this  case  as  well  and  it  is  the  Thomas  distribu¬ 
tion,  which  is  very  similar  in  character  to  the  Neyman  Type-A  distribution, 
and  to  the  fixed  multiplicative  Poisson  distribution  for  large  multiplica¬ 
tion  parameter. ^ 

(ii)  Luminescence  Noise  in  Electro-Optic  Devices 

Consider  an  application  in  which  we  wish  to  observe  photon  arrivals  by 
using  a  photomultiplier  tube  or  solid-state  detector  in  an  ionizing  radiation 
environment.  Examples  might  be  in  surveillance,  high-altitude  astronomy,  the 
guidance  system  of  a  spacecraft,  or  in  the  aftermath  of  a  nuclear  detonation. 
In  these  cases,  the  description  provided  above  may  be  characteristic  of  the 
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COUNTING  DISTRIBUTION  p(n) 


Figure  4.  Photon  counting  distribution  p(n)  vs  number  of  photon  counts  n. 

Data  (dots)  represent  radioluminescence  photon  registrations 
from  the  glass  faceplate  of  a  photomultiplier  tube,  induced  by 
90Sr_90Y  g- 

rays.  The  counting  time  T  *  400  ps.  The  experimental 
count  mean  and  variance  are  85.89  and  429.58,  respectively. 

The  solid  curve  represents  the  Neyman  Type-A  theoretical  counting 
distribution  with  the  same  values  of  count  mean  and  variance 
(a  *  4.0).  The  Poisson  distribution  with  mean  85.89  (indicated 
by  arrow)  is  shown  as  the  dashed  curve. 
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noise  rather  than  of  the  signal.  Several  authors  have  discussed  sources  of 

noise  in  photomultiplier  tubes  in  the  ionizing  radiation  environment  of  space. 

v 

Such  noise  may  arise  from  several  mechanisms  such  as  luminescence  and  Cerenkov 
radiation  in  the  photomultiplier  window;  secondary  electron  emission  from  the 
window,  photocathode,  and  dynodes;  Bremsstrahlung,  in  turn,  causing  such 
secondary  electron  emission;  cosmic-ray  bursts;  and,  of  course,  thermionic 
emission  dark  current.  These  effects  clearly  degrade  both  the  dynamic  range 
and  the  photometric  accuracy  of  low-light-level  measurements,  and  therefore 
must  be  properly  modeled.  It  is  evident  from  the  experimental  results  reported 
in  the  previous  subsection  that  the  SNDP  provides  a  sound  point  of  departure 
for  modeling  a  number  of  these  sources  of  noise  (luminescence  will  be  the 
dominant  source  of  noise  in  many  applications).  Similar  effects  will  occur 
in  solid-state  devices. 

Photomultiplier  tubes  are  often  preferred  to  solid-state  detectors  in 
such  applications  precisely  because  of  their  relative  radiation  hardness.  The 
radiation  hardness  of  optical  fibers, in  particular,  and  of  optical  systems 
more  generally,  has  received  substantial  attention  recently.  Aside  from  their 

(Q) 

many  other  advantages,  optical  fiber  and  signal  processing  systems  have 
the  additional  allure  of  immunity  to  electromagnetic  pulse  (EMP)  interfer¬ 
ence,  As  is  clear  from  our  earlier  discussion,  an  additional  factor 

that  must  be  considered  in  using  electro-optic  systems  in  ionizing-radiation 
environments  is  the  induced  multiplicative  luminescence  noise. 
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(iii)  Cathodoluminescence  Noise 


Cathodoluminescence  is  an  important  process  in  which  a  beam  of  accelerated 
electrons  incident  on  a  luminescent  material  (e.g.,  a  phosphor)  induces  the 
emission  of  light.  It  is  the  process  responsible  for  the  radar, 

television  monitor,  and  oscilloscope  images.  Conditions  for  the  validity  of 
the  SNDP  in  describing  cathodoluminescence  emission  are  that  the  electron  beam 
reaching  the  sample  be  represented  as  a  homogeneous  Poisson  point  process,  and 
that  each  primary  event  have  associated  with  it  an  impulse  response  function 
h(t)  that  governs  the  rate  of  production  of  Poisson  optical  photons.  Cathodo¬ 
luminescence  may  therefore  be  represented  by  the  physical  model  in  Fig.  3 
if  we  simply  replace  the  block  labeled  "Poisson  ionizing  radiation"  by  a  block 
entitled  "Poisson  electron  beam." 

The  counting  and  time  statistics  can  then  be  obtained  by  using  techniques 
to  be  outlined  in  the  next  section.  The  description  is  completely  characterized 
by  u  and  h(t)  and,  again,  in  the  limit  of  counting  times  long  in  comparison 
with  the  time  scale  of  h(t),  the  counting  distribution  will  reduce  to  the  Neyman 
Type-A  for  arbitrary  h(t). 

If  the  function  h(t)  is  a  decaying  single-time-constant  exponential,  each 
electron  entering  the  sample  induces  an  exponentially  decaying  flux  of  lumines¬ 
cence  photons.  In  this  special  case,  the  results  will  be  the  same  as  those 
derived  by  van  Rijswijk^^  *  via  an  entirely  different  method.  It  appears 
that  van  Rijswijk's  result  can  be  used  only  for  luminescence  that  decays  pre¬ 
cisely  exponentially  whereas  our  general  result  will  be  valid  for  arbitrary 
h(t).  In  particular,  we  can  easily  account  for  the  finite  "decay"  (rise)time 
of  the  luminescence  emission. 

The  statistics  for  the  time  to  the  first  cathodoluminescent  photon  arrival 
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were  experimentally  measured  by  van  Rijswijk^^  for  a  YVO^:Eu^+  phosphor  and, 

more  recently,  by  Timmermans  and  Zijlstra^^  for  ZnS:Ag.  In  Fig,  5  we  present 

3+ 

the  normalized  forward  recurrence  time  data  (dots)  for  YVO^:  Eu  obtained  by 
(14) 

van  Rijswljk  with  an  incident  electron  energy  of  10.2  keV,  The  mean  photon 

4-1  3  ■ 

counting  rate  was  2.54  x  10  s  and  the  mean  dark  counting  rate  was  1.28  x  10  s 

The  solid  curve  is  a  plot  of  the  expected  theoretical  forward  recurrence  time 

probability  density  P^(t)  for  an  exponentially  decaying  shot-noise  pulse,  with 

t  -  39.7  ps,  a  *  9.6,  and  t  >0.5  ms.  The  effects  of  dark  noise  were  not 

c 

negligible  in  this  case,  and  were  accounted  for  in  the  solid  curve.  The  experi¬ 
mental  data  clearly  deviate  from  the  straight-line  behavior  representing  the 
exponential  interval  distribution  characteristic  of  a  Poisson  process,  and 
so  does  the  theoretical  result.  Clearly,  cathodolumlnescence  photons  are  not 
Poisson  distributed. 

Intensity  fluctuations  for  cathodolumlnescence  were  first  investigated 
by  means  of  the  spectral  properties  of  excess  noise  in  the  photomultiplier 
anode  current.  The  theoretical  expression  for  the  spectrum  in  the 

general  case  (viz.,  for  arbitrary  h(t)  and  arbitrary  T)  can  be  derived  using 
techniques  we  propose  to  develop.  The  result  for  an  exponentially  decaying 
shot-noise  pulse,  in  a  fora  appropriate  when  measurements  are  made  by  a  real¬ 
time  spectrum  analyzer,  turns  out  to  be  Lorentzian,  in  agreement  with  previously 
derived  results. 

(14) 

Recent  excess  noise  measurements  were  performed  by  van  Rijswljk  for 
YV0^:Eu^+  and  by  Timmermans  and  Zijlstra^^  for  ZnS:Ag.  In  Fig.  6  we  present 
the  ZnS:Ag  experimental  relative  excess  noise  as  a  function  of  frequency  (dots)1 
obtained  by  Timmermans  and  Zijlstra.  The  solid  curve  is  a  Lorentzian  which 
clearly  provides  an  excellent  fit  to  the  data,  lending  confirmation  to  the 
exponential  decay  shape  and  the  luminescent  center  lifetime.  Use  of  the  SNDP 
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NORMALIZED  FORWARD  RECURRENCE  TIME  r/<r> 


Figure  5.  Dots  represent  the  relative  frequency  of  the  time  to  the  first 
cathodoluminescent  photoelectron  event  for  a  YVO^:Eu^+  phosphor 
with  the  exciting  electron  energy  fixed  at  10.2  keV.  The  mean 
photon  counting  rate  was  2.54  x  10^s-1  and  the  mean  dark  counting 
rate  was  1.28  x  lO^s-*.  The  solid  curve  is  a  plot  of  the  SNDP 
theoretical  probability  density  function  <t>Pj^(t)  vs  the  nor¬ 
malized  forward  recurrence  time  t/<t>  for  an  exponentially 
decaying  shot-noise  pulse,  with  t  ■  39.7  ps,  a  ■  9.6,  and  tc  * 
0.5  ms.  Dark  noise  has  been  accounted  for  in  the  theoretical 
curve.  Data  adapted  from  Fig.  6  of  the  paper  by  van  Rijswijk. 

The  data  demonstrate  that  cathodoluminescence  photons  are  clearly 
not  Poisson  distributed. 
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Fig.  3  of  the  paper  by  Timmermans  and  Zljistra. 


would  permit  us  to  obtain  the  excess  noise  for  an  arbitrary  impulse  response 
function. 

Substantial  progress  has  also  been  made  toward  the  implementation  and 
analysis  of  a  frequency-shift  keyed  (FSK)  fiber-optic  communication  link.  Work 
accomplished  in  the  areas  of  optimal  (single-threshold)  processing  schemes, 
a  new  technique  for  computing  performance,  and  the  demonstration  of  a  wideband 
FSK  fiber-optic  link  are  summarized  below. 

(1)  Single-Threshold  Detection  for  an  FSK  Fiber-Optic  System. 

It  is  well  known  that  optimal  processing  in  photocounting  optical  com¬ 
munications  systems  consists  of  the  likelihood  ratio  test  (LRT)  applied  to  an 
observed  quantity  at  the  detector  output.  In  practice,  however,  the  LRT  may 
be  difficult  or  Impossible  to  implement.  We  have  derived  a  remarkably  simple 
optical  processor  that  replaces  the  complicated  LRT  for  a  broad  class  of  systems. 

The  results  were  applied  to  an  FSK  photocounting  fiber-optic  communica¬ 
tion  system  in  which  the  noise  was  considered  to  be  of  the  Poisson  or  Poisson- 
avalanche  types  (corresponding  to  a  PIN  or  avalanche  photodiode  detector) . 

For  the  PIN  detector,  single- threshold  processing  was  shown  to  be  optimal  for 
several  decision  rules.  For  the  avalanche  photodiode,  it  was  shown  that  single¬ 
threshold  processing  did  not  necessarily  apply.  Examples  of  avalanche  photodiode 
receivers  were  single-threshold  processing  was  not  optical  were  presented.  Other 
cases  for  which  the  noise  contained  a  point  of  inflection  were  also  investigated, 
such  as  the  Polya  density  under  certain  conditions,  the  distribution  of  photo¬ 
counts  from  a  sinusoidally  modulated  single-mode  laser,  and  the  Cauchy  density. 


(2)  Probability  of  Error  for  an  FSK  Fiber-Optics  System 


It  has  become  the  habit  of  those  engaged  in  photocounting  optical 
receiver  design  who  are  unable  to  compute  bit-error  rates  exactly  to  resort 
to  Gaussian  approximations,  series  expansions,  and  special  bounds.  These 
techniques  all  trade  off  between  simplicity  and  accuracy. 

A  new  technique  for  evaluating  system  performance  was  developed  which 
is  as  simple  as  the  Gaussian  approximation,  yet  is  more  accurate.  This  tech¬ 
nique  is  based  upon  a  point  transformation,  "the  normalizing  transform",  which 
is  applied  to  the  noisy  receiver  output,  and  renders  its  statistics  more  nearly 
Gaussian.  The  resulting  output  is  then  more  accurately  approximated  by  con¬ 
ventional  Gaussian  techniques. 

The  form  of  the  normalizing  transform  was  determined  for  a  number  of 
probability  density  functions.  Including  the  Poisson,  binomial,  negative- 
binomial,  lognormal,  and  non-central  chi-squared.  In  particular,  the  normalizing 
transform  for  the  Poisson  was  a  simple  square-root  point  transformation.  The 
form  of  the  normalizing  transform  was  also  determined  for  a  broad  class  of 
receiver  outputs  x  where  the  variance  var  x  is  a  known  function  of  the  mean 
<x> 

var  x  -  D2  (<x>),  (1) 

We  have  shown  that  the  normalizing  transform  solution  for  this  case  is 

g(x)  -  / (1/D(x))dx  (2) 

This  resulting  family  of  random  variables  g(x)  was  nearly  normally  distributed 
even  though  the  untransformed  family  x  might  exhibit  excessive  skewness. 

The  normalizing  transform  technique  was  used  to  evaluate  the  performance 
of  fiber-optic  systems  with  avalanche  photodiode  receivers.  The  output  voltage 
of  a  typical  binary  digital  receiver  employing  an  avalanche  photodiode  was 
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represented  by  a  filtered  randomly  multiplied  Poisson  process. 

It  was  found  in  many  examples  that  the  normalizing  transform  provided 
an  approximation  that  was  as  good,  or  better  than  the  Gaussian  approximation 
for  evaluating  system  performance,  with  little  additional  computation.  The 
improvement  was  more  pronounced  when  the  pdf  was  skewed  and  exhibited  a  strong 
signal-dependence.  The  improvement  was  less  pronounced  in  limits  where  the 
pdf  approached  Gaussian  and  was  relatively  signal-independent. 

(3)  Wideband  FSK  Fiber-Optic  Experiment 

A  series  of  experiments  have  been  performed  in  which  a  digital  fiber¬ 
optic  communication  link  has  been  implemented  using  wideband  frequency  modula¬ 
tion.  System  performance  was  measured  and  compared  to  a  system  using  on-off 
keyed  modulation. 

By  now  the  merits  of  frequency  modulation  in  radio  communication  are 
widely  recognized.  For  digital  FM  (frequency-shift  keying  or  FSK),  when 
compared  to  on-off  keying  (00K)  on  a  peak-power  basis,  the  SNR  is  increased 
in  3 dB  for  both  coherent  and  noncoherent  detection.  For  FSK,  optimal  processing 
is  straightforward  to  implement.  In  contrast,  00K  requires  knowledge  of  all 
key  system  parameters  to  set  an  optimal  threshold.  These  parameters  may  be 
impossible  to  measure  exactly  and  may  change  with  time  causing  severe  degrada¬ 
tion  of  00K  performance. 

The  advantages  of  FM  in  optical  communication  are  even  more  striking, 
since  optical  bandwidth  is  plentiful  at  the  present  time,  and  the  potential 
bandwidth  expansion  is  enormous.  However,  all  attempts  to  utilize  FM  in 

/  1  Q\ 

optical  communications  have  been  thwarted  by  the  lack  of  tunable  sources.  ' 

(19)  (20) 

Various  subcarrier  FM  intensity  modulation  schemes  have  been  resorted  to. 

These  subcarrier  FM.  IM  schemes  are  inherentlv  limited  by  the  excursion  of  the 
intermediate  frequency  carrier  (e.g,,  25+5  MHz)  and  therefore  do  not  exploit 
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the  advantages  of  wideband  optical  FM. 

Recent  developments  In  frequency  tunable  lasers  have  now  made  FM  optical 

(21)— (23) 

communications  viable.  v  '  Olsson  and  Tang  have  demonstrated  electronic 

(24) 

tuning,  at  100  MHz,  of  a  room-temperature  cw  laser  over  a  5  tun  spectral  range. 

/  /0/\ 

Gigabit  data  rates  can  be  obtained  with  a  double-wavelength  mode-locked  laser.  '  ’ 

This  paper  reports  the  first  experimental  demonstration  of  a  wideband 
FSK  optical  fiber  communication  system.  The  center  wavelength  was  827  nm  with 
an  excursion  of  3  nm.  Pseudorandom  data  was  transmitted  at  a  rate  of  10  Mbits/ 
sec.  A  block  diagram  of  the  FSK  system  is  shown  in  Fig.  7 .  A  comparable  00K 
system  is  implemented  by  removing  one  of  the  lasers  and  removing  the  narrowband 
filter. 

Shown  in  Fig.  8  is  the  bit-error  rate  versus  total  optical  attenuation 
for  FSK  and  00K.  The  curves  represent  the  theoretical  results  obtained  using 
the  appropriate  estimates  of  system  parameters.  The  theoretical  optical  loss 
is  given  in  dB  on  the  upper  scale.  The  dots  are  experimental  data  obtained 
using  the  systems  described  in  the  previous  section.  The  experimental  optical 
loss  is  given  in  dB  on  the  lower  scale.  Approximately  12.9  dB  attenuation  is 
introduced  in  Che  connectors,  couplers  and  focussing  elements  of  this  system. 

The  total  optical  loss  also  includes  3.3  dB  from  the  spool  of  graded  index  fiber, 
plus  36.3  -  41.8  dB  introduced  using  neutral  density  filters.  This  attenuation 
is  equivalent  to  13.2  -  15  km  of  fiber  with  3  dB/km  loss.  Negligible  inter¬ 
symbol  interference  will  occur  over  these  distances  at  the  bit  rates  used  in 
our  system.  As  seen  in  Fig.  8  ,  the  experimental  results  are  in  good  agreement 
with  the  theory  in  terms  of  the  slope  and  relative  position  of  the  data  points. 

A  3  dB  improvement  in  sensitivity  of  FSK  over  00K  is  exhibited  over  6  orders  of 
magnitude  of  bit-error  rates.  The  experimental  losses  are  approximately  4.2  dB 
higher  than  predicted  by  theory,  as  seen  by  the  difference  between  the  experimental 
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Figure  7.  Block  diagram  of  FSK  fiber-optic  experiment. 
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(lower)  and  theoretical  (upper)  scales  in  Fig.  8,  This  is  accounted  for,  in  part, 
by  the  fact  that  average  rather  than  actual  estimates  of  device  parameters  were 

used  in  the  theoretical  calculations. 

In  summary,  the  superior  noise  immunity  of  wideband  frequency-modulation 
was  demonstrated  in  a  fiber-optic  communication  system.  If  suitable  tunable 
lasers  are  used,  FSK  can  provide  a  high-performance  alternative  to  conventional 
00K  for  fiber-optic  transmission. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under  Grants 
NSF-ENG78-26498  and  NSF-ECS-8006881. 
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II.  PHYSICAL  AND  PHOTOCHEMICAL  PROPERTIES  OF  ELECTRONIC  MATERIALS 

A.  PHENOMENOLOGICAL  MODEL  OF  GRAIN-BOUNDARY  TRAPPING  STATES  IN 
POLYCRYSTALLINE  SILICON  UNDER  OPTICAL  ILLUMINATION  * 

(E.  Poon,  W.  Hwrng  and  E.  S.  Yang) 

(JSEP  work  units  16  and  17,  1979-82) 

(Principal  Investigator:  W.  Hwang  (212)  280-3115) 

Electronic  transport  properties  of  polycrystalline  silicon  have  been 
receiving  much  attention  because  of  their  applicability  to  VLSI  bipolar  and 
MOS  integrated  circuits  and  their  potential  for  low-cost  terrestrial  solar 
energy  applications^’  ^  Photovoltaic  devices  fabricated  on  polycrystalline 
silicon  films  have  been  studied  by  many  investigators,^’  ^  It  appears  that 
the  grain  boundary  (GB)  interface  states  play  a  dominant  role  in  determining 
the  electronic  and  photovoltaic  properties  of  polycrystalline  silicon  by  acting 
as  traps  and  recombination  centers.  However,  the  fundamental  understanding 
of  these  GB  states  is  still  not  clear  and  substantial  research  effort  in  this 
direction  is  required.  In  this  work,  we  develop  a  new  theoretical  model  of 
the  grain  boundary  recombination,  charge  carrier  transport  and  electrostatics 
under  the  assumptions  of  Gaussian  energy  distributions  of  interface  states, 
and  unequal  capture  cross-sections  for  electrons  and  holes.  We  also  demonstrate 
that  the  Gaussian  energy  distributions  can  be  reduced  to  other  distributions 
as  special  cases,  in  particular  to  exponential  energy  distribution,  uniform 
distributions  and  discrete  energy  levels.  Carrier  generation  is  expected  to 
introduce  another  degree  of  freedom  and  provide  valuable  additional  insight 
into  the  optical  dependence  of  the  electrostatic  effects  of  grain  boundary 
Interface  states. 

The  energy  band  diagram  of  a  polycrystalline  semiconductor  in  the  dark 
and  under  illuminated  conditions  in  the  neighborhood  of  two  grain  boundaries 
is  shown  in  Fig.  1,  At  the  grain  boundary  plane  are  localized  interface 
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states,  which  in  general  are  statistically  distributed  in  energy  over  the 
energy  gap.  These  states  have  a  distribution  determined  by  the  nature  of  the 
disorder  in  the  silicon  bond  structure  such  as  dangling  bonds,  and  by  the 
local  electronic  potential  fluctuations  producing  stress  fields  in  the  regions 
of  structural  irregularities.  They  may  also  arise  from  chemical  impurities 
between  adjacent  grains.  Together  dangling  bonds  and  local  structural  strain 
can  produce  an  inhomogeneous  broadening  of  the  discrete  levels  in  the  gap  if 
their  effects  are  sufficiently  severe.  As  a  feasible  approximation  of  the 
energy  spectra  for  these  localized  states,  Gaussian  energy  distributions  are 
proposed.  These  states  are  classified  into  two  types,  "donor-like"  and 
"acceptor-like"  interface  states.  The  energy  distribution  of  the  resultant 
donor-like  and  acceptor-like  states  htD(E)  and  htA(E)  are  expressed  by  the 

htAw  •  /(jitsA-  •** 

and 

HtD 

htD(E)  “  TTMy  “P 

where  HtA  and  are  the  total  number  of  the  localized  states  per  unit  area, 

SA  and  SD  are  distribution  parameters  and  E  A  and  EtD  are  the  mean-values  of 

the  acceptor-like  and  donor-like  distributions,  respectively. 

Figure  1(b)  shows  the  band  diagram  in  the  absence  of  illumination.  The 

GB  diffusion  potential  barrier  height,  V^q,  can  be  calculated  by  equating  the 

depletion- region  charge  accumulated  in  the  interface  states.  The  numerical 

value  of  VJrt  have  been  calculated  for  both  uniform  and  nonuniform  energy 
au 
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distributions  of  interface  states  by  several  workers.^  ’  ^ ^ 

The  presence  of  optical  lllimlnation,  as  shown  in  Fig.  1(c),  separates 

the  quasi-Fermi  levels  E,  (0)  and  E,  (0)  for  electrons  and  holes,  and  diffusion 

rn  tp 

potential  decreases  from  its  dark  value  V^q.  Between  ejq(0)  and  E^p(0)  the 

GB  interface  states  act  as  recombination  centers  and  the  occupancy  of  the 

recombination  center  is  determined  by  the  relative  magnitude  of  capture 

cross-sections  and  by  the  relative  concentration  of  electrons  and  holes  at 

GB,  according  to  Shockly-Read-Hall  (SRH)  statistics. ^  The  recombination 

centers  at  the  GB  are  either  neutral  or  charged  depending  on  whether  they 

are  occupied  donor-like  or  acceptor-like  states,  or  whether  they  are  empty 

acceptor-like  or  donor-like  states.  In  the  steady  state  a  flux  of  electrons 

and  an  equal  flux  of  holes  of  the  GB  is  sustained  by  the  photogeneration 

of  these  carriers  in  the  neighborhood  of  the  GB.  Thus  flux  is  equal  to  the 

steady  state  recombination  rate  at  the  grain  boundary. 

The  carrier  transport  in  polycrystalline  semiconductors  is  dominated  by 

the  thermionic  emission  of  majority  carrier  over  the  GB  diffusion  potential 

barrier.  The  conductivity  of  polycrystalline  films  is  given  by  Eq.  (9)  for 
kT  (5), (8), (11) 

V  <  — 

q 


o 


q2d N, 


(2irm*kT) 


1/2 


exp(- 


qv 


do. 
kT  ' 


(3) 


The  optical  illumination  modifies  the  GB  potential  barrier  V.,  as  discussed 

d 

above,  and  V  decreases  from  its  dark  value  V,  .  Hence  the  enhancement  of 
a  do 

conductivity  due  to  the  optical  illumination  (photoconductivity)  is  given  by 
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Equation  (4)  can 


where  ■  V^q  -  AV^ ,  and  it  is  assumed  that  o  >>  °(jarjt* 
also  be  rewritten  as 


where 


a  .  -  q  N.  u  -c 
ph  ’  d  eff 


eff 


SA. 


(2irm*kT) 


1/2 


qV. 

exp(-  -£-)  exP(— > 


(5) 


(6) 


In  these  expressions,  ^e£f  i3  the  effective  majority  carrier  nubility  which 
reflects  the  optical  modulation  of  the  GB  potential  barrier.  depends 

exponentially  on  AV^  which  in  turn  is  controlled  by  the  optical  illumination. 

In  other  words,  the  photoconductivity  in  polycrystalline  semiconductors  can 
be  a  "majority"  carrier  effect  due  to  optical  mobility  modulation  which 
differs  in  principle  from  that  of  crystalline  semiconductors. 

The  density  and  distribution  of  the  GB  trapping  states  strongly  depend 
on  the  polycrystalline  samples  preparation,  passivation  and  annealing  processes. 
The  typical  experimental  results  of  the  GB  trapping  states  observed  in  poly¬ 
silicon  are  summarized  in  Table  1.  As  mentioned  earlier,  the  Gaussian  distri¬ 
bution  is  capable  of  representing  different  types  of  distributions 

when  suitable  conditions  are  imposed  on  it.  With  the  approximation 

S„,  S.  «  kT,  the  GB  interface  state  distributions  given  by  Eqns.  (1)  and  (2) 

become  confined  to  two  discrete  energy  levels  (or  two  6-function  distributions 

in  energy).  With  the  approximation  EtD  -  Ey  -  0,  E^  ■  E£  -  Eg  and  SA,  SD  >  kT, 

Eqns.  (1)  and  (2)  approach  two  exponential  distributions  in  which  the  densities 

of  acceptor-like  and  donor-like  states  decrease  exponentially  into  the  energy 

gap  from  the  conduction  band  edge  and  the  valence  band  edge,  respectively. 

(12) 

In  fact,  in  this  case  the  distribution  becomes  a  U-shaped  distribution. 

For  S_,  S  >>  kT,  Eqns.  (1)  and  (2)  approach  two  uniform  distributions. 

U  A 
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The  various  GB  trapping  states  observed  in  polysilicons 


L_ 


The  following  major  conclusions  can  be  drawn  from  Che  work  reported 
herein.  Photoelec tronic  processes  at  grain  boundaries  (GBs)  In  polycrystalline 
semiconductors  have  been  studied  theoretically  and  experimentally.  A  gener¬ 
alized  Gaussian  distribution  of  GB  interface  states  has  been  proposed,  which 
can  be  specialized  by  suitable  conditions  to  express  other  distributions 
as  well.  The  photoconductivity  in  polycrystalline  silicon  is  observed  to 
be  strongly  dependent  on  the  optical  illumination  intensity,  showing  a  pro¬ 
nounced  increase  with  intensity  which  we  believe  to  be  correlated  with  a 
decrease  in  GB  diffusion  potential.  It  is  concluded  that  the  photoconductivity 
in  polysilicon  is  the  result  of  a  "majority”  carrier  mobility  modulation 
effect  and  hence  differs  in  principle  from  photoconductivity  mechanisms  of 
crystalline  semiconductors. 

*This  research  was  also  supported  by  the  Department  of  Energy  under  Contract 
DOE/SERI  #xw-l-1272-l. 
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B.  LASER  PROCESSED  SCHOTTKY  BARRIERS* 

(C.  M.  Mu,  E.  S.  Yang,  T,  0.  Sedgwick, +  R,  T.  Hodgson  ) 

(JSEPwork  units  16  and  17,  1979-82) 

(Principal  Investigator:  E.  S.  Yang  (212)  280-3120) 

Recently,  we  have  reported  a  thin-film  structure  using  polycrystalline 

silicon  as  the  energy  absorbing  layer  to  form  an  Al-nSi  Schottky  barrier, ^ 

In  contrast  with  the  p-n  junction  structure  without  the  capped  layer  reported 
(2) 

by  Parks  and  Rose,  the  required  laser  power  in  our  experiment  is  substan¬ 
tially  lower  and  our  diodes  have  low  turn-on  voltage,  indicating  they  are 
similar  to  a  Schottky  barrier  rather  than  p-n  junction.  In  our  experiments, 
both  a  Nd-YAG  pulse  laser  and  a  cw  argon  laser  have  been  used.  Because  the 
argon  laser  scanning  system  is  more  flexible  and  the  result  is  readily  repro¬ 
ducible,  we  have  used  it  to  investigate  the  electrical  behavior  of  diodes  as 
a  function  of  the  laser  power.  This  paper  reports  the  results  of  our  study. 
Aluminum  dots  with  a  thickness  of  2000  8  and  an  area  of  1.13  ram^  were 

deposited  on  (lll)n/n+  silicon  wafers.  The  epitaxial  layer's  resistivity  is 
14  -3 

7  &cm  (N^  •  7x10  cm  )  and  its  thickness  is  15  pm.  These  dots  were  covered 
with  an  E-beam  evaporated  polysilicon  film  of  approximately  100  8  thick.  In 
the  backside  of  the  wafer,  gold  was  evaporated  to  provide  a  good  ohmic  contact. 
Subsequently,  the  devices  were  subjected  to  laser  irradiation.  A  cw  argon 
laser  (4880  and  5145  8) ,  operated  at  a  power  level  between  12  and  22  W,  is 
focused  to  a  spot  size  of  50  y.  The  spacing  between  adjacent  scanning  lines 
is  20  urn,  so  that  complete  coverage  of  the  surface  by  the  laser  beam  is  en¬ 
sured. 

The  substrate  temperature  was  set  at  350°C  (the  minimum  temperature  in 

the  experimental  setup) ,  and  the  stepping  motor  movement  (with  a  speed  of 

2  cm/sec)  gives  rise  to  overlapping  spots  of  laser  irradiation.  The  average 

5  2  3  2 

power  reaching  the  target  is  (6-10)  x  10  W/cm  [(3.5)  x  10  J/cm  ]. 
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Laser-irradiated  dots  have  the  traces  of  the  laser  beam,  showing  that  the 
surface  has  been  melted  and  then  solidified.  The  surface  morphology  was  re¬ 
ported  previously. ^  When  the  sample  is  laser  irradiated  with  a  power  higher 
than  20  W,  the  evaporated  silicon  film  shows  ripping  between  the  aluminum  dots. 

In  the  areas  with  aluminum,  less  energy  is  absorbed  and  no  significant  damage 
is  observed.  The  ripping  of  silicon  film  around  the  periphery  of  aluminum 
dots  causes  a  large  leakage  current  and  unstable  I-V  characteristics.  Therefore, 
we  will  only  show  the  results  of  the  dots  prepared  by  power  levels  below  20  W. 

The  room-temperature  current-voltage  (I-V)  curves  of  the  dots  processed 
with  various  laser  powers  are  shown  in  Figs.  1  and  2.  At  high  forward  bias, 
the  dots  display  different  ideality  factors  which  may  be  used  to  identify  the 
various  conduction  mechanisms.  The  dots  of  12  W  has  an  I-V  curve  similar  to 
diodes  prepared  by  thermal  annealing ^ ^  and  the  ideal  factor  is  substantially 
less  than  two.  The  other  diodes  show  higher  saturation  current  and  larger 
ideality  factor.  The  reverse-biased  current  appears  to  be  a  generation  current, 
and  it  increases  as  the  irradiation  power  is  raised.  The  kink  in  the  reverse 
current  of  the  20-W  diode  is  not  understood,  but  it  is  likely  to  be  caused  by 
leakage  paths  through  the  diode  periphery.  From  the  generation  current  compo¬ 
nent,  the  lifetime  of  minority  carrier  in  the  16-W  diode  is  estimated  to  be 

14  -3 

200  nsec  which  corresponds  to  a  trapping-center  density  of  5x10  cm  . 

The  temperature  effect  of  the  I-V  characteristics  is  investigated  between 
200°  and  300°K.  The  reverse  saturation  current  Jq  is  found  by  extrapolating 
the  forward  current  at  high  bias  (>  0.5  V)  to  the  current  axis.  By  plotting 
JQ  versus  1/T,  we  deduce  the  activation  energy  from  the  slope,  and  the  result 
is  shown  in  Fig.  3,  The  activation  energy  of  the  12 -W  diode  is  0.807  eV  and 
its  current  is  dominated  by  the  thermionic  emission  across  a  barrier. 

A  value  of  activation  energy  equal  to  half  of  the  forbidden  gap  is  found 
on  the  16-W  diode.  Indicating  the  current  is  dominated  by  the  space-charge  layer 
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Figure  1.  The  forward  I-V  characteristics  of  the  diodes  irradiated  by 
12W,  16W,  20W. 


SATURATION 


Figure  3.  Saturation  current  density  vs  the  inverse  of  the  temperature 
The  diodes  are  separately  irradiated  by  12,  16,  and  20W. 
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recombination  process.  The  16-W  diode  has  a  higher  value  of  than  that  of 
the  12-W  diode  because  its  ideality  factor  is  equal  to  two;  however,  the  forward 
current  of  the  former  is  lower  than  that  of  the  latter  at  the  high-bias  regime. 
With  over  16  W  of  laser  irradiation,  the  diode  current  has  an  opposite  tendency; 
i.e.,  the  larger  laser  power  causes  a  higher  forward  current  because  both 
recombination  and  thermionic  emission  are  increased.  A  detailed  discussion  will 
be  presented  in  the  next  section. 

2 

In  the  C-V  measurements,  we  finu  that  the  plots  of  1/C  versus  the  bias 

voltage  are  all  straight  line  except  for  the  20-W  diode.  From  the  C-V  data, 

we  obtain  the  value  of  the  intersection  voltage  ,  the  doping  density  N,, 

int  d 

and  the  Fermi  potential  which  are  listed  in  Table  I.  These  parameters  for 
the  20-W  diode  are  estimated  by  using  the  C-V  data  at  the  regime  of  small  re¬ 
verse  bias  (<  IV).  It  is  evident  that  the  laser  power  of  12  W  does  not  change 
the  substrate  doping  density.  In  other  words,  the  irradiation  power  is  consumed 
in  the  process  of  melting  the  aluminum  and  the  interface  region.  The  contact 

barrier  is  formed  by  recrystallizing  a  p+  layer  between  aluminum  and  the  silicon 
(3) 

substrate.  However,  the  measured  substrate  doping  is  changed  when  the  laser 
power  is  higher  than  14  W.  The  measured  doping  density  is  lower  than  the  original 
doping  when  the  laser  power  is  between  12  and  16  W,  and  it  increases  when  the 
laser  power  is  higher  than  16  W. 

When  the  surface  layer  of  silicon  under  the  aluminum  is  irradiated  and 

its  temperature  is  near  the  melting  point  of  silicon  (1410°C) ,  the  aluminum 

atoms  may  easily  diffuse  into  the  silicon  substrate.  Because  of  a  very  short 

heating  period  (^2.5  msec),  the  aluminum  atoms  will  remain  in  a  shallow  region. 

At  the  temperature  of  1400°C,  the  diffusivity  of  aluminum  in  silicon  is 
-10  2 

2x10  cm  /sec.  Thus  the  thickness  of  the  aluminum  doped  silicon  layer  is 
approximately  70  X.  If  a  higher  power  or  longer  dwell  time  is  applied,  the 
substrate  will  stay  at  high  temperature  longer  presumably  allowing  more  aluminum 
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Table  I 


The  deduced 

values 

from  the 

C-V  data 

at  T  *  295 

°K. 

cw  laser 
(W) 

power 

12 

14 

16 

18 

20 

N,  (1014 
a 

-3. 
cm  ) 

6.9 

5.13 

4.42 

10.2 

11.6 

WV> 

0.523 

0.562 

0.567 

0.554 

0.543 

V  (V) 

0.270 

0.277 

0.281 

0.261 

0.257 

n 

diffusion.  The  aluminum-diffused  layer  forms  a  negative  space-charge  region 
which  enhances  the  contact  barrier  of  the  aluminum  on  the  n-type  silicon  sub¬ 
strate.  The  n-type  Schottky  diode  with  a  shallow  p  layer  between  metal  and 
semiconductor  has  been  studied  by  Shannon^  and  Van  der  Ziel.^  Recently, 
we  have  presented  both  theoretical  and  experimental  works  on  the  thermally- 
treated  Al-nSi  diodes.  Assuming  a  uniform  doping  in  both  p-  and 

n-type  region,  N  and  N, ,  the  capacitance  per  unit  area  is  given  by 
a  a 


q(N  +  N  ) 

h  «  <j>  -  V  -  v  +  — - - - —  t 

pint  9BnQ  n  T  2e  z  » 


(1) 

(2) 


where  q  is  the  electron  charge,  e  is  the  dielectric  constant  of  silicon, 

is  the  Fermi  potential  in  the  n  side,  t  is  the  thickness  of  p+  layer, 

is  the  thermal  voltage,  and  <p  is  the  barrier  height  determined  by  the  dif- 

Bn0 

ference  of  the  metal  work  function  and  the  electron  affinity  of  the  semi- 

(4) 

conductor.  <p_  was  experimentally  found  to  be  0.075  eV.  <p.  _  is  the 
on-  int 

2 

intersection  voltagd  by  extrapolating  the  1/C  curve  to  the  voltage  axis. 
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The  barrier  height  at  zero  bias  can  be  deduced  by  using 


q(N  +  N  )  , 

*Bn  *  *int  +  VT  +  Vn  *  *BnQ  + - 2^  C  *  (3) 

Therefore  the  C-V  data  yield  both  the  doping  density  and  the  barrier  height 
at  zero  bias.  From  the  value  of  the  doping  density,  we  obtain  the  Fermi  poten¬ 
tial  at  room  temperature  and  then  calculate  the  barrier  height  using  Eq.  (3). 

The  result  is  shown  in  Fig.  4.  It  should  be  noted  that  Eq.  (3)  is  valid  only 

(4) 

if  N  is  much  larger  than  N,.  The  barrier  height  of  the  diode  after  laser 
a  d 

irradiation  increases  with  the  laser  power.  However,  the  barrier  height 

5  2 

reaches  a  maximum  value  of  0.87  eV  when  the  laser  power  is  16  W  (%  4x10  J/cm  ) 
and  it  then  decreases  with  increasing  power.  Examining  Eq.  (3),  we  find  that 
the  barrier  height  strongly  depends  on  the  doping  of  p+-layer  N  and  the  thick- 

c 1 

ness  t.  If  the  p  layer  is  mainly  due  to  the  precipitation  of  silicon  from 

the  AL-Si  solution  formed  by  the  laser  irradiation,  the  thickness  t  should 

increase  monotonically  with  the  irradiation.  Thus  the  barrier  height  lowering 

when  the  power  is  beyond  16  W  is  due  to  the  decrease  of  the  impurity  density 

N  .  Under  a  higher-  power  laser  irradiation,  the  silicon  surface  region  is 

possibly  melted  such  that  the  final  solidified  silicon  layer  on  the  surface 

does  not  have  the  same  quality  as  before.  Apparently,  the  precipitated  silicon 

contains  donor-like  imperfection  centers ^  which  compensate  the  effect  of 

ionized  aluminum  atoms  and  causes  the  reduction  of  the  barrier  height. 

From  C-V  data,  it  is  found  that  the  substrate  doping  is  first  reduced 

because  of  the  aluminum  diffusion.  Subsequently,  N,  reaches  a  value  higher 

d 
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than  the  original  epitaxial  doping  density  (7  x  10  cm  ).  To  get  a  clearer 
picutre,  we  consider  these  data  together  with  the  generation  current  of  the 
diode  under  a  reverse  bias.  The  lifetime  of  minority  carriers  in  the  16-W 
diode  is  estimated  to  be  200  nsec  which  corresponds  to  a  trapping-center 
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CW  ARGON  LASER  POWER  (W) 


4.  Barrier  height  vs  the  laser  power  after  the 
diode  is  irradiated. 
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density  of  5  x  10  cm  ,  a  value  approaching  the  original  doping  density. 

When  the  laser  irradiation  power  is  increased,  a  large  number  of  defects  can 
be  formed  in  the  substrate  to  increase  the  generation  current.  Because 
the  defects  are  donor-like, ^  the  effective  substrate  doping  density  becomes 
higher  than  the  original  value.  A  thorough  understanding  of  this  phenomenon 
requires  a  detailed  measurement  of  the  traps  which  is  being  investigated. 

Since  the  recombination  current  is  a  non-negligible  term  in  the  I-V 
characteristics,  if  the  diode  is  irradiated  by  high  power,  the  forward  current 
is  the  sum  of  the  thermionic  emission  and  recombination  currents.  At  room 
temperature  or  below,  the  ideality  factor  is  equal  to  two  and  the  recombina¬ 
tion  current  dominates  the  forward  current.  But  at  higher  temperatures  the 
thermionic  emission  dominates  such  that  the  ideality  factor  decreases.  This 
is  consistent  with  the  fact  that  the  activation  energy  of  the  20-W  diode  under 
forward  bias  has  a  value  of  0. 723  eV,  which  is  between  the  barrier  height  and 
half  of  the  forbidden  gap. 
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c.  SCANNING  LASER-SPOT  PHOTOCONDUCTIVITY* 

(P.  Panayotatos,  E.  S.  Yang,  W.  Hwang) 

(JSEP  work  units  16  and  17,  1979-82) 

(Principal  Investigator:  E.  S,  Yang  (212)  280-3120) 

In  this  work  a  scanning  laser  spot  technique^  ^  has  been  used  for 
the  extraction  of  information  about  the  grain  boundary  (GB)  recombination 
kinetics.  From  the  experimental  measurements,  the  value  of  the  GB  recombination 
velocity  can  be  determined  using  practically  no  assumption  of  parameter  values. 

At  the  same  time,  the  laser  spot  provides  the  illumination  level  to  which  the 
GB  is  subjected  and  modifies  both  the  electrostatics  and  the  kinetics  there. 

The  intensity  of  the  laser  was  therefore  varied,  spanning  more  than  two  orders 
of  magnitude,  so  that  a  quantitative  relationship  between  illumination  and 
recombination  velocity  could  be  extracted.  Also  described  is  a  method  based 
on  the  transition  from  low-to  high-level  injection  of  minority  carriers  for 
the  estimation  of  the  quasi-Fermi  level  separation  corresponding  to  each  illumina¬ 
tion  level,  again  without  assumption  of  numerical  values  for  any  parameters 
other  than  the  minority  carrier  diffusion  coefficient. 

The  scanning  light  spot  technique  is  a  powerful  experimental  method  for 
the  extraction  of  information  about  the  electrostatics  and  the  recombination 
kinetics  at  the  GB.  A  derivation  of  the  correlation  between  the  macroscopic 
parameters  and  the  fundamental  properties  of  the  GB  is  given  in  Ref.  (3). 

Briefly,  under  the  assumption  stated  above,  i.e.  the  length  of  the  two  depletion 
regions  is  much  smaller  than  the  grain  size  and  for  low-level  injection  and  for 
a  voltage  drop  V  <kT/q  at  the  GB,  the  thermionic  emission  equation 


J 


*  2 
A  T 


exp 


Ef  +  qyD  \  v 

kT  /  kT 


(1) 
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is  usedv  '  to  correlate  the  measured  photovoltage  with  the  intrinsic  parameters. 
The  modified  form  for  the  photo conductance  is  given  by 


AG  -  bgT 


1  +  S 


where  x  is  the  distance  of  the  light  spot  from  the  GB.  b-  depends  on  several 

experimental  parameters  as  well  as  on  the  dark  GB  diffusion  potential  but  its 

value  does  not  enter  into  our  numerical  calculation.  Equation  (2)  predicts 

a  linear  log  AG  vs  x  relationship  for  as  long  as  remains  constant.  Furthermore 

and  always  appear  as  a  ratio  which  will  remain  constant  even  if  Lq  changes 

as  condition  change  from  low-  to  high-level  injection  of  minority  carriers. 

The  sanples  of  Wacker  polycrystalline  silicon  were  cut  into  2  mm  x  10  mm 

slabs  and  chemically  etched  to  a  mirror-like  surface.  Au  ohmic  contacts  were 

then  vacuum  deposited  at  the  ends  of  the  (p-type)  semiconductor  slab,  and  the 

sample  was  mounted  on  an  x-y  translation  stage  (see  insert.  Fig.  1).  The 

experimental  set-up  consisted  mainly  of  a  5  mW  He-Ne  laser  whose  beam  was  chopped 

and  then  focused  on  sample  to  a  10  ym  dia.  spot.  By  using  a  lock-in  amplifier 

and  a  reference  signal  from  the  copper,  a  signal  proportional  to  the  change  of 

the  sample  conductance  was  measured  and  plotted  vs  position  of  the  laser  spot 

on  the  sample.  A  tungsten  lamp  was  used  for  visual  observation  of  the  grain 

boundaries  under  consideration,  as  well  as  for  focusing  the  laser  spot.  The 

(4) 

schematic  for  the  experimental  set-up  is  shown  in  the  insert  of  Fig.  1. 

Scanning  of  the  sample  with  the  laser  spot  produced  conductance  vs 

position  curves  typically  as  those  of  Fig.  1,  where  it  was  verified  by  visual 

observation  that  the  peaks  (largest  increases  in  conductance)  coincided  with 

the  points  where  the  laser  spot  fell  exactly  on  a  grain  boundary.  Not  all  grain 

(3) 

boundaries  gave  rise  to  such  peaks,  in  agreement  with  earlier  observations. 


The  magnitude  of  the  largest  peak  also  varied  markedly  from  sample  to  sample. 
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ILLUMINATION  LEVEL : 


contacts)  and  the  experiment  set-up. 


The  grain  boundary  that  produced  the  largest  such  change  in  sample  conductance 
was  thus  singled  out. 

The  change  of  sample  conductance  vs  position  of  the  laser  spot  in  the 

neighborhood  of  the  grain  boundary  was  recorded  for  ten  different  illumination 

levels  chosen  by  the  introduction  of  filters  in  the  path  of  the  laser  beam, 

the  highest  level  corresponding  to  0.12  mW  of  light  energy  falling  on  the 

sample.  More  than  two  orders  of  magnitude  in  illumination  level  were  thus 

covered  and  plots  such  as  those  of  Fig.  2  were  obtained.  Care  was  taken  to 

keep  any  effect  of  the  DC  biasing  voltage  across  the  grain  boundary  negligible 

so  that  any  possible  symmetry  in  the  grain  boundary  electrostatics  was  not 

disturbed  and  that  the  assumptions  necessary  for  the  derivation  of  Eqn,  (2) 

(3) 

are  correct.  The  ratio  of  the  dark  sample  resistance  to  a  series  resistance 
was  such  that  the  voltage  across  the  whole  sample  was  of  the  order  of  kT/q. 

That  this  was  indeed  the  case  and  that  the  grain  boundary  under  consideration 
did  not  present  any  assymetry  was  verified  by  comparing  the  plots  of  laser  spot 
scans  for  both  polarities  of  the  DC  bias.  Curves  A'  in  Fig.  2  represents  such 
a  response. 

It  is  immediately  apparent  that  there  is  a  suppression  of  the  change  of 
conductance  when  the  laser  spot  is  very  close  or  exactly  on  the  grain  boundary. 
This  suppression  is  more  pronounced  for  higher  illumination  levels  (see  curves 
A,  B,  C  and  compare  to  curve  D  in  Fig.  2.  It  is  assumed  that  this  suppression 
is  due  to  the  radical  increase  of  the  grain  boundary  recombination  velocity 
Sq  that  takes  place  when  the  laser  spot  is  within  one  diffusion  length  away  from 
the  grain  boundary.  The  above  allows  the  effect  of  the  grain  boundary  to  be 
separated  from  any  bulk  effects.  Wacker  material  is  especially  suited  for  such 
measurements.  The  diffusion  lengths  in  Wacker  are  about  one  order  of  magnitude 
larger  than  the  diameter  of  the  laser  spot  and  the  grains  are  of  the  order  of 
mm  so  that  they  in  turn  are  much  larger  than  the  diffusion  length,  allowing 
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scan  In  the  ne tghborhood  of  the  (train  boundary  for  four  different  illumination 
nte  the  suppression  of  the  peak  at  higher  Illumination  levels.  Curve  A 
a  scan  Identical  to  that  of  curve  A  except  with  the  polarity  of  the 


rather  extended  scans  around  a  particular  grain  boundary,  with  no  interference 

from  other  grain  boundaries.  The  suppression  will  not  be  apparent  for  materials 

(3) 

where  the  diffusion  lengths  are  comparable  to  the  laser  spot  diameter  as  was 
verified  by  unfocusing  the  laser  beam  and  scanning  the  same  area  of  the  sample; 
much  smoother  curves  were  then  observed. 

We  can,  therefore,  recognise  two  areas  in  the  response:  one  with  the  laser 
spot  very  close  or  on  the  grain  boundary,  where  grain  boundary  effects  predominate 
and  where  there  is  a  significant  grain  boundary  recombination  velocity,  and  another 
with  the  spot  further  away  from  the  grain  boundary  where  bulk  effects  are  domi¬ 
nant  and  the  effect  of  the  recombination  velocity  is  negligible  so  that  the  re¬ 
sponse  should  follow  Gqn.  (2)  with  S  *  0,  or 

L 

AG  *  b  lit  exP(“x/Ln)*  (3) 

n  1‘ 

Indeed,  for  the  latter  region,  a  linear  log  AG  vs  x  plot  was  obtained  (see  Fig.  3) 

from  the  slope  of  which  the  diffusion  length  L  can  be  determined.  The  devia- 

n 

tion  from  linearity  for  the  values  of  x  near  the  grain  boundary  indicates  a 
nonconstant  recombination  velocity  S  which  suppresses  the  photoconductance  as 
represented  by  Eqn.  (2),  As  the  laser  spot  moves  closer  to  the  grain  boundary, 
photogenerated  carriers  modify  the  carrier  concentration  at  the  grain  boundary 
and  therefore  also  the  recombination  velocity  there. 

The  value  of  S  can  be  determined  from  the  plots  of  Fig.  3  in  the  following 
way:  Extrapolating  the  lineary  regions  (bulk  effect  dominated  part  of  the  plot) 
to  x  *  0,  which  corresponds  to  the  laser  spot  exactly  on  the  grain  boundary, 
determines  what  the  value  of  AG  would  be  if  there  were  no  interface  recombina¬ 
tion.  In  other  words  we  measure  AG^-(O)  *  bL  /2D  .  For  S  4  0,  on  the  other 

n  n 

hand,  from  Eqn,  (6)  we  obtain  that 
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CHANGE  OF  CONDUCTANCE  AG (mS) 


DISTANCE  OF  LASER  SPOT  FROM  GRAIN  BOUNDARY  x(mils) 


Figure  3.  Plot  of  the  logarithm  of  the  photoconductance  versus 
position  of  the  laser  spot  for  different  illumination 
levels  for  a  typical  sample.  Data  have  been  taken  from 
plots  such  as  those  of  Figure  3.  More  than  two  orders 
of  magnitude  in  illumination  were  covered.  Not  all 
illumination  levels  are  shown  here. 
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AG(O) 


2D  1  +  SL  /2D 
n  n  n 


so  that  the  ratio  of  the  observed  value  of  the  change  in  conductance  (AG(0)) 
to  the  value  obtained  by  the  extrapolation  (AG"(0))  has  the  value: 


AG(0) 

AG"(0) 


_ 1 _ 

1  +  SL  /2D 
n  n 


(4) 


Equation  (4)  can  therefore  be  used  for  the  experimental  determination  of  S 
independently  of  such  parameters  as  Richardson's  constant,  interface  state  den¬ 
sities  and  distributions  or  values  of  illumination  intensities.  The  recombina¬ 
tion  velocity  can  be  thus  directly  plotted  as  a  function  of  illumination.  However, 
since  the  values  of  the  diffusion  lengths  suggest  a  transition  from  low-  to 
high-level  minority  carrier  injection,  an  estimate  can  be  made  of  the  quasi-Fermi 
level  separation  AE^.  Thus  S  can  be  presented  as  a  function  of  Ef  as  well  as  of 
illumination.  In  Fig.  4  we  present  the  electron  diffusion  length  (obtained 
from  the  slopes  of  the  linear  segments  of  Fig.  3  as  above)  as  a  function  of  il¬ 
lumination.  It  is  apparent  that  there  is  a  transition  from  low-  to  high-level 

14 

injection  of  minority  carriers  corresponding  to  the  1.2  x  10  photons/sec  level. 

The  low  level  injection  value  of  is  slightly  larger  than  expected^  rising 

to  about  twice  its  value  at  high  illumination.  Given  that  the  minority  carrier 

concentration  n  at  this  illumination  level  is  about  equal  to  the  doping  concen- 

15  -3 

tration  (NA  *  2.5  x  10  cm  )  and  electron  concentration,  the  quasi-Fermi 
level  separation  AE^  can  be  calculated.  S  is  presented  as  a  function  of  AEf 
as  well  as  illumination  in  Fig,  5.  S  is  shown  to  increase  monotonically  with 
illumination  (or  quasi-Fermi  level  separation)  and  is  also  shown  not  to  saturate, 
even  for  the  very  large  AE^  produced  by  the  high  illumination  levels  used. 
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Figure  4.  Plot  of  the  electron  diffusion  length  vs  illumination  level. 

Ln  values  are  deduced  from  the  slopes  of  the  lines  like  those  of 
Fig.  3. 


Using  the  photoconductivity  scanning  technique,  we  have  been  able  to  measure 
the  minority-carrier  diffusion  length  within  the  grain  and  the  recombination 
velocity  at  the  grain  boundary.  The  diffusion  length  is  obtained  from  the  slope 
of  the  semilog  plot  of  the  photoconductivity  signal  vs  distance  from  the  grain 
boundary.  Extrapolation  of  the  semilog  plots  allows  us  to  determine  the  in¬ 
fluence  of  the  GB  recombination,  thus  enables  us  to  calculate  the  GB  recombina¬ 
tion  velocity.  This  method  is  attractive  in  that  it  does  nc  t  have  to  make  a 
priori  assumptions,  is  simple  and  self-contained.  The  numerical  values  we 
obtained  are  more  than  a  factor  of  ten  greater  than  Scholl's  results^  and 
an  order  of  magnitude  smaller  than  Seager's  estimation,^  An  important  con¬ 
sideration  in  our  method  is  the  spot  size  of  the  laser  beam  which  should  be 
small  compared  with  the  diffusion  length.  The  diffusion  length,  in  turn,  should 
be  much  smaller  than  the  grain  size.  It  was  found  that  the  GB  recombination 
velocity  increases  monotonically  with  illumination  and  does  not  saturate  even 
at  high  intensity, 
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D.  COULCMBIC  SCATTERING  AT  SEMICONDUCTOR  GRAIN  BOUNDARIES  * 

(C.  M.  Wu,  E.  S.  Yang) 

(JSEP  work  units  16  and  17,  1979-82) 

(Principal  Investigator:  E.  S.  Yang  (212)  280-3120) 

In  the  last  few  years,  research  activities  have  been  extensive  in  poly¬ 
crystalline  semiconductors  because  of  their  importance  in  very  large  scale 
integrated  structures  as  well  as  their  potential  in  low  cost  photovoltaic 
applications.  While  the  technology  of  polycrystalline  materials  is  reason¬ 
ably  well  established,  the  practical  significance  of  grain  boundaries  as  the 
dominating  influence  in  the  electronic  properties  of  these  semiconductors  is 
recognized  only  recently.  In  general,  a  grain  boundary  GB  has  two  different 
effects  on  the  material.  First,  it  has  interfacial  trapping  states  possibly 
resulting  from  mis-matched  dangling  bonds,  impurity  diffusion,  and  segregation, 
or  defects.  Thase  states  are  minority-carrier  sinks  such  that  the  overall 
minority-carrier  lifetime  and  diffusion  length  are  lowered.  Second,  a  potential 
barrier,  along  with  a  depletion  layer,  is  formed  at  the  GB  so  that  the  majority- 
carrier  current  is  impeded.  The  GB  potential  barrier  reduces  the  mobility 
and  increases  the  effective  resistivity  of  the  semiconductor.  This  is  un¬ 
desirable  for  most  applications  not  only  because  of  the  higher  resistance  but 
also  because  of  the  difficulty  in  the  reproducibility  of  the  resistance  value. 
This  work  addresses  the  second  problem  although  the  magnitude  of  the  GB  poten¬ 
tial  has  obvious  implications  for  minority-carrier  recombination  and  transport. 

A  single  energy-band  diagram  of  the  GB  consisting  of  two  back-to-back 
Schottky  barriers  is  shown  in  Fig,  1.  Most  papers  in  the  literature  make 
use  of  this  model  in  which  the  current  is  governed  by  the  thermionic  emission 
as  expressed  in  the  following  form:^ 

J^2  »  A*T2  exp[-q($B  +  V  )/kT]  (1) 
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where  A*  is  Richardson's  constant,  is  the  barrier  height  measured  from  the 
conduction  band,  Vn  is  the  Fermi  potential,  k  is  Boltzman's  constant,  q  is 
the  electronic  charge,  and  T  is  the  temperature  in  degree  Kelvin.  When  the 
appropriate  parameters  are  substituted  in  Eq,  (1),  however,  the  current  is 
found  to  be  significantly  different  from  experimental  data  such  that  an  attenua¬ 
tion  factor  is  required  in  Eq.  (1) .  ^  The  current  is  reduced  by  this  factor 
to  provide  more  reasonable  correlation.  The  experimentally  determined  value 
of  the  attenuation  factor  ranges  from  0.06  to  0.2.^  ^  The  large  varia¬ 
tion  is  probably  due  to  the  integrity  of  the  GB  interface  as  determined  by  the 
technology  and  processes.  A  model  has  been  proposed  by  Lu  et  al.^^  to 
interpret  this  phenomenon  by  adding  a  scattering  potential  which  gives  better 
agreement  between  theory  and  experimental  data.  However,  Lu  et  al.  did  not 
explain  the  origin  of  the  scattering  potential  except  that  they  did  refer  to 
the  possiblity  of  phonon  scattering.^ 

In  this  report,  we  show  that  a  two-dimensional  plane  of  fixed  charge 
can  contribute  to  significant  scattering  such  that  the  majority-carrier  trans¬ 
port  is  reduced.  This  Coulombic  scattering  is  found  to  be  much  more  important 
than  phonon-assisted  processes  at  room  temperature. 

The  scattering  of  an  electron  in  the  Coulomb  field  of  a  charged  center 
is  described  by  using  the  Rutherford  model.  Let  us  assume  that  the  GB  surface 
is  on  the  xy  plane  and  the  incident  electron  approaches  the  surface  with  a 
velocity  v  and  an  angle  8.  We  can  choose  a  system  of  coordinates  such  that 
the  electron  initially  travels  in  the  direction  parallel  to  the  xz  plane  as 
shown  in  Fig.  2.  The  initial  velocity  can  be  expressed  as 

■  v(sin8,  0,  cos8)  in  the  xyz  system,  (2) 
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For  simplicity,  an  ionized  charge  is  taken  as  located  at  the  origin  of  the 
system.  After  scattering,  the  electron  moves  in  a  direction  that  depends  upon 
not  only  the  incident  angle  8  and  the  initial  velocity  but  also  on  the  distance 
of  the  incident  path  and  the  center  of  the  ionized  charge,  b.  The  scattering 
angel  8  is  given  by^ 

cot  0/2  =  (4irejn*/q2)  v2b,  (3) 

where  £  is  the  semiconductor  permittivity  and  m*  is  the  effective  mass  of 
electron.  When  we  choose  the  x'y'z'  system  obtained  by  rotating  the  xyz 
system  around  the  y  axis,  as  shown  in  Fig.  2,  such  that  the  z'  axis  has  the 
same  direction  as  the  initial  velocity,  the  electron  approaches  the  ionized 
charge  and  is  headed  toward  the  point  (b  cos<J>,  b  sin<t>,  0)  on  the  x'y'  plane. 
Thus,  the  position  on  the  xy  plane  to  which  the  electron  points  is  (b 
cos  <>sec  8,  b  sin  <t>  ,  0).  If  the  collision  is  elastic,  the  final  velocity 
is  given  by 

vf  ■  v(sin8  cos<(>  cos8  +  cos0  sin0,  sin0  sin<)>  , 

cos8  cosQ  -  sin8  cos$  sin8  ),  (4) 

When  the  electron  is  scattered  back  from  the  xy  plane,  the  z  component  of  the 
final  velocity  must  be  less  than  zero.  The  cross  section  of  the  ionized 
charge  can  be  deduced  by  the  condition  of  <  0.  In  other  words, 

cot®  <  cos4>  tan8,  (5) 

Substituting  Eq.  (3)  into  Eq.  (5)  and  replacing  cos<f>  in  terms  of  the  variables 
in  the  xyz  system,  we  obtain 

(x2  +  y2)1^2  <  x  sin8  +  q2/4irem*v2  (6) 
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Figure  2.  (a)  x’y'z'  system  obtained  by  rotating  the  xyz 

system  along  the  y  axis  such  that  the  initial  velocity 
is  along  the  z'  axis,  (b)  Schematic  diagram  of  the 
relationship  between  the  initial  velocity  and  the 
final  velocity  due  to  the  charge  scattering. 


Thus,  the  cross  section  is  given  by 

o(9,v)  -/  dx  dy,  (7) 

s 

where  S  is  the  area  with  a  boundary  determined  by  Eq.  (6).  After  a  two- 
dimension  integration,  the  cross  section  is  given  by 

cr(0,v)  *  (q^/16ire^m  ^)v  ^cos  ^9.  (8) 


It  is  noticed  that  9  is  the  incident  angle  so  that  the  cross  section  in 
Eq.  (8)  is  different  from  the  cross  section  in  terms  of  the  deflection  angle  9 
in  the  calculation  of  the  impurity  scattering  in  the  bulk  single-crystal 
semiconductor . 

Therefore,  the  transmission  probability  T  of  the  electron  across  the 

grain  boundary  is  a  function  of  9  and  v.  When  the  interface  charge  has  a 
_2 

density  of  N(cm  ),  T  is  given  by 


T  =  0  if  o_  >  1/N, 

JJ 

-  1  -  No  0_  <  1/N, 

D  D 


(9) 


In  obtaining  the  thermionic-emission  current  across  the  GB,  it  is  necessary 
to  consider  this  scattering  effect  due  to  the  GB  charge.  By  assuming  that 
the  GB  surface  is  on  the  xy  plane,  the  electron  current  from  the  first  grain 
(z  <  0)  to  the  second  grain  (z  >  0)  is 


l->2 


2 q  f  r  ,  3  +  V  ) 

*3  /  /  /  dk  vz  exp< - - ~)  T(8,v), 


<2ir) 


(2ir) 


3  5*  /  dk  k  exP<- 


jkV 

2m*kT 


) 


(10) 


x  f  dd  alnd  coad  T  (9,v), 
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where  -K  is  h/2ir  and  h  is  the  Planck's  constant.  The  upper  and  lower  bounds 
of  the  intergration  are  determined  by  <  1/N.  When  the  electron  trajectory 

is  normal  to  the  GB  plane,  i.e. ,  0  *  0,  the  cross  section  is 

°B  ^V’  9  “  0)  M  (q^/16ire2m  2)v  *.  (11) 

The  transmission  coefficient  is  zero  when  the  electron  energy  is 

E  -  1/2  m*v2  <  q2/  N/8(/ir)e.  (12) 

In  other  words,  the  electron  cannot  pass  through  the  GB  under  the  condition 

of  Eq.  (12).  This  situation  is  corresponding  to  a  critical  energy  of  q<|>  , 

sc 

that  is  defined  as 

<(>  *  q/N/8(/ir)e  (13) 

SC 

when  the  electron  has  a  kinetic  energy  less  than  q4>  ,  no  thermionic  emission 

SC 

can  occur  (i.e.,  we  assume  tunneling  is  negligible).  On  the  GB  surface,  both 

positive  and  negative  charges  may  cause  the  scattering  of  the  approaching 

15  -2 

carrier.  The  maximum  surface-state  density  of  the  silicon  GB  is  10  cm 

which  is  approximately  equal  to  the  number  of  surface  atoms.  Thus,  the  maximum 

scattering  potential  is  340  mV.  The  net  GB  interface  charge  density  is  likely 

12  -2 

to  be  in  the  order  of  10  cm  . 

With  the  transmission  coefficient  given  by  (9),  Eq.  (10)  may  be  simplified 
and  we  can  obtain  the  thermionic  emission  current  with  a  scattering  factor 
n,  which  is  given  by 

Jl_*2  “  lA*T2exp[-q(<j)B  +  Vn)/kT], 

where 
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n  »  exp(- 


00 


In  Fig.  3,  we  show  n  versus  the  charge  density  N  at  300°K.  Note  that  the 

attenuation  of  the  thermionic  emission  strongly  depends  on  the  charge  density. 

12  13  -2 

For  the  charge  density  between  10  and  10  cm  ,  scattering  factor 
is  between  0.5  and  0.1.  These  values  are  a  factor  of  2  larger  than  reported 
experimental  values.  It  is  conceivable  that  there  are  both  positive  and 
negative  charges  at  the  GB  such  that  they  compensate  each  other.  Although 
the  resulted  net  charge  is  smaller  the  scattering  calculation  must  include 
both  types  of  charges.  For  this  reason  the  model  presented  in  this  letter  is 
reasonable. 

The  large  fixed  charge  is  not  surprising  in  that  they  must  exist  to 
establish  a  large  potential  barrier  at  the  grain  boundary.  Experimentally 
the  potential  barrier  has  been  determined  by  conductance  and  capacitance 
>■  *asurements,  ^  and  a  barrier  height  of  approximately  1/2  eV  at  a  doping 

level  of  10*^  is  consistent  with  the  results  reported  here.  As  for  the  origin 
of  the  fixed  interface  charge,  it  is  likely  that  they  are  introduced  by  the 
material  processing  and  discontinuity  at  the  GB.  In  particular,  impurity  atoms 
have  a  tendency  to  congregate  at  the  GB  since  they  diffuse  much  faster  along 
the  planes  of  defects.  In  addition  oxygen  precipitation  has  also  been  found 
to  enhance  the  electrical  activity  of  the  GB,^^  another  probable  origin  of 
charged  centers. 

*This  work  is  also  supported  by  DOE/SERI  Contract  xz-<j>-9226. 
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E.  INTERFACE  INTERACTION  OF  THE  Al-nSi  SYSTEM* 

(C.  M.  Wu,  E.  S.  Yang) 

(JSEP  work  units  lo  and  17,  1979-82) 

(Principal  Investigator:  E.  S.  Yang  (212)  280-3120) 

We  present  a  study  of  the  electrical  properties  of  Al-nSi  diodes  fabricated 
by  thermal  annealing  in  a  furnace  of  600°C  with  different  heating  periods  and 
cooling  rates.  Since  the  annealing  temperature  is  above  the  eutectic,  a  re¬ 
crystallized  layer  is  formed  in  the  aluminum-silicon  interface.  This  is  in 
contrast  with  previous  studies  where  a  low-temperature  (300-500°C)  annealing 
is  employed. ^  ^  Because  the  ideality  factor  is  far  from  unity,  the  I-V 
data  cannot  quantitatively  indicate  the  effects  of  the  recrystallized  silicon 
layer.  We  consider  theoretically  the  influence  of  a  thin  Al-doped  silicon 
on  the  capacitance  characteristics  of  the  diode.  Both  I-V  and  C-V  measurements 
are  made.  The  C-V  data  are  used  to  obtain  the  effective  barrier  height. 

Aluminum  dots  with  a  thickness  of  2000  X  are  deposited  on  (lll)n/n+ 
silicon  wafers  that  have  an  epitaxial  layer  with  a  resistivity  of  0.4  ft  cm 
and  a  thickness  id  *  15  pm.  These  dots  are  covered  with  an  e-beam-evaporated 
silicon  film  of  approximately  300  X  thick.  In  the  backside  of  the  wafer,  gold 
is  evaporated  to  provide  a  good  ohmic  contact.  The  structure  of  the  diode 
before  heat  treatment  is  shown  in  Fig.  1  Subsequently,  the  diodes  are  thermally 
annealed  in  a  hydrogen  environment  of  600°C  for  1  h.  It  is  important  to  note 
that  the  cap  silicon  layer  and  hydrogen  gas  are  employed  to  prevent  the  aluminum 
film  from  reacting  with  the  environment  during  the  heat  treatment.  To  study 
the  effect  of  cooling  rate,  we  cool  the  samples  either  with  a  temperature  rate 
l°C/min  (diodes  TR1)  or  quenchinp  h«  withdrawing  the  furnace  away  from  the  sample 
(diodes  TRQ,  the  cooling  rate  is  _<out  100°C/min),  A  control  sample  TRO  is  also 
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Figure  1.  Structure  of  the  Al-nSi  diode  before  heat  treatment. 

prepared  by  annealing  the  specimen  for  5  min  at  300°C  and  then  turning  off 
the  furnace.  This  cooling  process  takes  10  h  for  the  specimen  to  reach  room 
temperature. 

Another  series  of  diodes  were  made  to  test  the  dependence  of  the  diode 
performance  on  the  annealing  time.  T20,  T40,  and  T60  are  diodes  prepared, 
respectively,  by  20,  40,  and  60  min  of  thermal  annealing  on  the  same  structure 
as  shown  in  Fig.  1  The  substrate  of  these  diodes  are  (100)  oriented. 

The  metal-semiconductor  junction  with  a  heavily  doped  interfacial  layer 
has  been  studed  by  Shannon^  and  van  der  Ziel,^  The  models  they  proposed 
are  limited  to  a  description  of  the  I-V  properties  and  show  the  dependence  of 
the  ideality  factor  and  the  barrier  height  on  the  thickness  and  the  doping 
density  of  the  interfacial  layer.  However,  the  measurement  of  the  barrier 
height  using  the  I-V  curve  requires  the  ideality  factor  to  be  unity.  But  the 
space-charge  recombinations  current  gives  rise  to  a  2-kT-current  component. 
Therefore  it  is  ambiguous  to  deduce  the  barrier  height  only  from  the  1-V 


60 


data,  since  the  recrystallized  Al-nSi  Schottky  barriers  are  nonideal  diodes. 

The  space-charge  capacitance  is  derived, ^  and  it  is  used  to  determine  the 
effective  barrier  height. 

The  I-V  curves  of  the  diodes  T20,  T40,  T60  are  shown  in  Fig.  2.  We  find 
that  the  diode  characteristics  do  not  change  significantly  with  the  different 
annealing  times.  In  addition,  the  C-V  curves  yield  the  same  barrier  height. 

We  conclude  that  the  p+  layer  is  not  formed  by  the  diffusion  of  aluminum  into 
the  silicon  substrate,  otherwise  the  thickness  and  doping  density  would  have 
to  be  a  strongly  time-dependent  process.  Hence,  the  modulation  of  the  barrier 
height  is  more  likely  to  be  caused  by  the  recrystallized  silicon  layer  doped 
with  aluminum.  In  Ref.  7,  we  showed  that  the  barrier  height  increases  with  the 
thickness  of  the  evaporated  aluminum,  and  it  saturates  at  0.93  eV  when  the 
aluminum  layer  is  thicker  than  1500  X.  The  saturation  is  probably  due  to  the 
limitation  of  the  silicon  diffusivity  in  aluminum,  such  that  the  silicon  atoms 
dissolved  deep  in  aluminum  are  locked  in  and  unable  to  precipitate  out.  At 
the  annealing  temperature  (600°C) ,  the  maximum  silicon  molar  fraction  of  Al-Si 
liquid  solution  is  13.2%. ^  Thus  270  &  of  silicon  may  be  dissolved  in  1500  X 
of  aluminum.  However,  the  dissolved  silicon  is  not  completely  precipitated  out 
of  the  aluminum  after  the  cooling  down,  so  that  the  final  residue  aluminum 
contains  some  portion  of  silicon  atoms.  The  eutectic  composition  (11.3%  of  Si) 
is  most  likely  a  molar  fraction  in  aluminum  residue,  such  that  only  40  X  of 
silicon  is  recrystallized. 

Figures  3  and  4  separately  show  the  I-V  and  C-V  charactersitics  of  the 
diodes  TRO,  TR1,  and  TRQ,  As  mentioned  previously,  TRO  is  a  control  diode 
which  is  fabricated  without  the  silicon  recrystallization  at  the  interface. 

The  diode  has  a  barrier  height  determined  by  the  work  function  of  aluminum 
and  the  electron  affinity  of  silicon.  The  C-V  data  show  that  ^  is  equal 
to  -  0.187  V,  such  that  the  barrier  height  is  0.075  eV,  The  low  barrier 
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C-V  characteristics  of  TRO,  TR1,  and  TRQ. 


height  of  the  diode  is  correlated  to  the  poor  I-V  characteristics  (Fig.  3). 

We  define  the  barrier  height  of  TRO  to  the  4>_  ,  which  is  assumed  to  be  invariant 

bno 

in  the  aluminum-silicon  contact.  From  the  C-V  data,  we  obtain  the  intersection 
voltage  $  of  TR1  and  TRQ  to  be  0,693  and  0,734  eV.  The  barrier  heights  are 
calculated  to  be  0,96  and  0,99  eV,  respectively.  The  difference  is  clearly 
caused  by  different  cooling  rates.  The  faster  cooling  gives  rise  to  a  higher 
intersection  voltage.  It  is  reasonable  to  say  that  more  impurities  stay  in 
the  recrystallized  silicon  in  faster  cooling.  In  this  paper,  we  have  qualita¬ 
tively  correlated  the  cooling  rate  and  the  doping  density  of  the  recrystallized 
p+  layer  which  modulates  the  barrier  height  of  Al-nSi  diode.  In  fact,  with 

postannealing  at  a  low  temperature  ('''300°C),  the  oversaturated  aluminum  in 

+  (3) 

the  recrystallized  p  layer  will  precipitate  out.  Although  the  diode  pre¬ 
pared  by  quenching  process  has  a  large  barrier  height,  the  diode  may  not  be 
stable  if  it  operates  in  a  temperature  cycling  environment. 

The  thermal  annealing  of  the  aluminum-silicon  contact  yields  an  Al-nSi 
diode  with  a  p+-silicon  interfacial  layer  which  increases  the  barrier  height. 
Cooling  from  the  temperature  above  the  eutectic  point  (600°C) ,  the  diode  has 
a  barrier  height  which  does  not  vary  with  the  annealing  time,  but  changes 
with  the  cooling  rate.  To  avoid  the  inaccuracy  in  determinating  the  barrier 
height  from  the  I-V  curve,  we  have  used  the  C-V  measurements,  which  clearly 
show  the  effect  of  the  doping  density  and  the  thickness  of  the  p+  layer. 

It  is  concluded  that  the  diffusion  of  aluminum  into  the  silicon  is  not  signif¬ 
icant,  and  a  fast  cooling  process  causes  a  high  aluminum  concentration  in  the 
recrystallized  silicon  layer, 

*This  work  is  also  supported  by  DOE/SERI  Contract  xz -$-9226. 
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F.  LASER  PHOTOCHEMISTRY  FOR  MICROELECTRONICS* 

(R.  M.  Osgood) 

(JSEP  work  unit  5,  1982-1985) 

(Principal  Investigator:  R.  M.  Osgood  (212)  280-4462) 

Introduction 

This  report  covers  a  new  program  which  was  initiated  in  September* 

1981.  The  objective  of  the  research  is  1)  to  develop  an  understanding  of  laser  - 
initiated  interface  reactions  and  2)  to  apply  the  results  of  this  research  to 
problems  in  microelectronic  fabrication.  Recently  it  has  been  shown  that 
laser  chemistry  at  or  near  surfaces  can  form  the  basis  of  a  new  direct-write 

processing  technology  for  microelectronics.  In  addition,  an  understanding  of 
laser  surface  reactions  may  lead  to  the  realization  of  photochemical  catalysis 
reactions,  for  bulk  chemical  synthesis. 

In  the  first  year  of  this  program,  we  have  concentrated  on  treating 
theoretically  light-induced  surface  reactions  and  on  setting-up  a  laboratory 
which  can  analyze  both  the  theoretical  and  applied  aspects  of  laser  surface 
chemistry. 

Theoretical 

While  surface  chemistry  has  been  the  subject  of  a  vast  amount  of  pure  and 

applied  research,  there  have  been  only  limited  studies  of  photochemistry  at  fluid- 

solid  interfaces.  A  solid  surface  can  be  expected  to  modify  the  usual  homogeneous 

photochemical  reactions  through  changes  in  the  photodissociation  dynamics,  shifts 

in  the  spectral  features  of  adsorbed  species,  and  formation  of  new  chemical 

species,  The  surface  can  enter  in  a  more  subtle  but  equally  important 

manner  by  modification  of  the  local  optical  field  at  the  interface.  Depending  on 

the  surface  morphology  and  its  optical  properties,  the  local  electric  field  can 

be  enhanced  or  diminished  considerably  with  respect  to  its  value  in  the  gas 

(2) 


alone . 


The  enhancement  results  from  the  increase  in  the  local  optical  field  in 
the  vicinity  of  microscopic  metal  features  which  exhibit  plasma  resonances  for 

frequencies  close  to  that  of  the  incident  electric  field.  Recently,  related 

(3) 

studies  of  this  phenomena  for  Raman  scattering,  linear  scattering,  and  linear 

absorption^)  and  luminescence^)  have  been  reported.  For  Raman  scattering,  an 

increase  of  10^  is  seen  in  the  spontaneous  Raman  cross-section  of  molecules 

(3) 

adsorbed  on  metal  surfaces.  Surface-enhanced  Raman  scattering  (SERS)  has  been 
observed  for  photolithographically  fabricated  artificial  surfaces.^) 

We  have  shown  that  this  increase  in  the  local  field  can  also  cause  an 
enhancement  by  many  orders  of  magnitude  in  the  rate  of  an  interface  photochemical 
reaction.  The  effect  is  not  only  significant  for  species  adsorbed  directly  on  the 
surface,  but  it  can  also  be  an  important  factor  for  heterogeneous  reactions  which 
are  initiated  in  the  gas  phase.  The  increase  in  local  field  is  greatest  for  metal 
particles  with  the  appropriate  shape  and  material  characteristics  to  exhibit  a 
plasma  resonance  near  the  frequency  of  the  optical  source.  This  effect  has  many 
similarities  with  SERS;  however,  both  the  UV  operating  wavelength  and  the  polariza¬ 
tion  typical  of  photochemical  reactions  give  rise  to  different  conditions  in  cal¬ 
culating  the  enhanced  reaction  rate.  Under  typical  conditions,  reaction  enhance¬ 
ments  of  'v 100  are  found  for  adsorbed  layers  undergoing  single-photon  photodis- 

4 

sociatlon;  multiple-photon  reactions  exhibit  enhancements  of  ■vlO  .  The  random 

size  and  shape  which  can  be  expected  on  real  but  artificially  prepared  surfaces 

do  not  significantly  alter  these  ratios. 

Surface  enhancement  of  chemical  reactions  has  been  recently  considered 

(2) 

by  Brus  and  Nitzan.  They  point  out  that,  although  it  is  reasonable  to  expect 
an  Increase  in  the  rate  of  reactions  on  surfaces  due  to  the  enhancement  of  the  local 
optical  field,  radiationless  transfer  to  the  solid  substrate  may  be  expected  to 
reduce  the  quantum  yield  of  the  bound-free  event.  However,  the  nature  and  degree 
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of  energy  transfer  to  the  surface  during  a  photodissociative  event  is  currently 
not  well  understood  either  theoretically  or  experimentally.^  Recent  experi- 

/Q\ 

ments'  indicate,  for  example,  that  metal  alkyls  can  be  photodissociated  readily 
on  metal  surfaces.  To  simplify  our  treatment  of  the  influence  of  the  local 
field  enhancement,  we  have  not  included  energy  transfer  from  the  adsorbed  phase 
to  the  surface  or  the  reverse  process  in  our  calculation.  This  simplification 
does  not  influence  our  calculation  for  the  case  of  enhanced  fields  in  the  gas 
or  in  upper-multilayer-adsorbed  species,  where  energy  transfer  is  greatly  reduced. 

Experimental 

The  experimental  facilities  necessary  for  carrying  out  the  work  described 
above  should  consist  of  a  modern  optical  laboratory  including  laser  sources 
and  optical  spectroscopic  techniques.,  and  facilities  for  fabricating  and  testing 
microelectronic  devices.  During  this  year  we  have  set  up  an  optical  laboratory 
with  cw  and  pulsed  laser  sources  at  both  visible  and  infrared  wavelengths.  These 
sources  have  been  used  for  both  direct-writing  and  pulsed  laser  experiments. 

In  addition,  a  computer-controlled  UV  spectrometer  has  been  equiped  with  an 
orgometallic  gas  handling  system  and  temperature  controlled  cell.  The  apparatus 
permits  observation  of  the  UV  adsorption  spectrum  of  weakly  bound  surface  films 
of  adsorbed  material. 

Construction  of  a  clean  room  for  fabrication  and  testing  of  microelectronic 
devices  has  been  started.  This  facility  will  be  equipped  with  diffusion  furnances, 
wet  and  dry  clean  benches,  and  other  facilities  for  fabrication  of  small  scale 
electronic  devices. 

We  anticipate  that  the  room  will  be  completed  by  the  first  part  of  July. 

*This  research  was  also  supported  by  the  Army  Research  Office  under  Grant 
#DAAG29-82-K-0089  and  Air  Force  Office  of  Scientific  Research/DARPA  under 
Grant  F-49620-82-K-0008. 
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G.  A  NEW  CAPILLARY  NOZZLE  ION  SOURCE:  THE  CAPILLARITRON * 

(Toomas  Allik,  George  Flynn) 

( JSEP  work  unit  7,  1982-85) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

Within  the  past  few  years  the  use  of  ion  beams  in  the  field  of  micro¬ 
electronics  has  become  extremely  important.  Due  to  great  interest  in  ion 
beam  lithography,  ion  implantation  and  ion  sputtering,  immense 
technological  advances  have  been  achieved  in  the  field  of  ion  sources. 

For  microelectronic  devices  and  circuit  processing,  the  key  parameters  in 
determining  the  supremacy  of  a  given  ion  source  are:  brightness,  emittance, 
beam  divergence,  ion  current  production,  energy  distribution  and  ionization 
ef f iciency . ^  While  in  most  cases  the  construction  of  ion  sources  has  been 
tailored  to  the  semiconductor  industry,  these  ion  sources  can  be  utilized 
in  other  disciplines.  For  example,  ion-molecule  reactions  and  ion  spectroscopy 
are  two  fields  of  active  research  that  have  benefited  tremendously  from  the 
discoveries  of  new  ion  sources.  Our  objective  is  to  develop  new  and  improved 
methods  for  the  production  of  metastable  molecular  and  atomic  ion  beams ;  to 
characterize  the  translational,  vibrational,  rotational  and  electronic  energy 
distributions  and  relaxation  rates  of  such  beams,  and  to  study  the  chemical 
reactions  of  these  metastable  ion  species  with  both  gases  and  surfaces. 

Since  the  original  reports  of  a  new  capillary  nozzle  ion  source  called 
(2)  (3) 

the  "Capillaritron" ,  ’  we  have  concentrated  our  efforts  on  characterizing 

the  ion  beam.  The  Capillaritron  is  a  simple  device  consisting  of  a  metal 

(4) 

nozzle  and  an  extractor  electrode.  The  nozzle  is  electrically  isolated  by 
a  dielectric  from  the  source  gas  chamber  and  is  brought  to  voltage  levels  on 
the  order  of  several  kilovolts.  The  extractor  electrode  is  kept  at  ground 
potential,  and  is  placed  in  the  plane  of  the  nozzle  tip  or  slightly  in  front  of 
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(downstream)  the  nozzle.  When  gas  is  fed  to  the  nozzle  and  a  threshold  voltage 
is  applied,  a  plasma  is  formed  at  the  nozzle  tip,  which  is  visible.  The  theory 
of  operation  can  be  found  in  References  2  and  3. 

The  Capillaritron  was  installed  in  a  differentially  pumped  molecular 
beam  apparatus  so  that  backing  pressures  as  high  as  five  atmospheres  could  be 
applied  behind  the  nozzle.  Since  the  molecular  beam  contains  an  extranuclear 
quadrupole  mass  spectrometer,  a  home-made  Faraday  cup  and  time  of  flight 
chopper,  the  identity  as  well  as  the  energy  spread  of  the  ions  produced  can 
be  easily  determined. 

Various  gases  have  been  ionized  such  as  Ar,  Ne,  Xe,  He,  0^  and  N^. 

A  comparison  of  the  overall  ionization  efficiency  of  each  gas  shows  that  under 
equivalent  conditions  (backing  pressure,  orifice  diameter,  etc.),  the  number 
of  ions  produced  scales  inversely  as  the  ionization  potential.  With  regard  to 
flow  rates,  the  threshold  for  breakdown  ionization  in  all  cases  decreased  as 
the  backing  pressure  was  increased.  In  addition,  as  the  orifice  diameter  was 
increased  to  as  large  as  500  microns,  ion  production  increased.  Since  the 
supersonic  nozzle  was  placed  on  a  motorized  translational  stage  the  nozzle  to 
extractor  distance  could  be  varied  as  much  as  one  inch.  When  the  distance  was 
set  to  a  maximum,  the  only  effect  noticed  was  that  a  few  hundred  volts  more  were 

required  to  reach  threshold  ionization.  Mass  spectral  analysis  for  Ar  produced 

-t-  +2  + 

a  majority  of  the  Ar  ion,  but  a  significant  abundance  of  the  Ar  and  Ar^ 

species  was  present.  Neon  exhibits  similar  features  as  is  shown  in  Figure  1. 

Energy  analysis  of  the  ion  beam  by  two  techniques:  time-of-flight  and  floating 

the  potentials  on  the  quadrupole  rods  shows  a  large  energy  spread.  These 

results  are  in  good  agreement  with  recent  work  of  Goebel  et  al.^  Presumably 

the  large  energy  spread  is  due  either  to  charge  exchange  with  slow  neutrals  in 

the  ion  extraction  region, ^  or  non-uniform  acceleration  in  the  ionization 
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BACKGROUND  SPECTRA 


Figure 


region.  Change  transfer  collisions  of  the  type: 


Ar+(fast)  +  Ar(thermal)  -*•  Ar(fast)  +  Ar+(thermal) 

have  been  seen  since  the  time-of-flight  exhibits  a  double  peak.'  (see  Figure  2). 
Finally,  the  Capillaritron  has  the  ability  to  produce  negative  ions.  When  a 
negative  voltage  is  applied  to  the  nozzle  and  0^  is  sent  through  the  apparatus, 
a  plasma  is  visible  that  is  quite  different  in  appearance  from  that  observed 
with  a  positive  voltage.  In  addition,  no  positive  ions  could  be  detected  with 
the  mass  spectrometer.  Minor  modifications  to  the  Channeltron  will  make  it 
possibile  for  us  to  look  at  negative  ions  in  the  future. 

In  conclusion,  the  Capillaritron  is  a  powerful  ion  source  capable  of 
producing  large  ion  currents  (>  1  mA)  with  small  beam  divergence  (2-8  degrees) . ^ 
The  only  drawback  of  this  source  is  that  it  produces  beams  having  a  large  energy 
spread.  In  the  future  we  plan  to  use  several  novel  techniques  to  produce  re¬ 
active  ion  beams  and  to  characterize  the  chemical  and  energy  distribution  of 
these  beams.  Experiments  designed  to  probe  the  structure  and  relaxation  pro¬ 
perties  of  the  ion  beams  will  then  be  performed.  Finally,  the  chemical  reac¬ 
tions  of  these  ions  with  gaseous  neutrals  and  their  surface  etching  properties 
will  be  investigated  in  an  attempt  to  achieve  a  fundamental  understanding  of 
ion-neutral  interactions  and  the  reactive  ion  etching  process. 
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Figure  2.  Time  of  flight  of  argon  utilizing  the  Farday  cup. 

The  distance  of  travel  is  approximately  55  cm. 
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"Investigation  of  the  Capillaritron  Ion  Source  For  Electric  Propulsion", 
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(4)  Phrasor  Scientific,  Inc.  Manufactors  the  Capillaritron,  price  $1,950. 

Since  our  molecular  beam  appratus  already  included  gas  and  high  voltage 
feedthrus  and  translational  stages,  it  was  cost  effective  to  design 
our  own  Capillaritron. 

(5)  D.  M.  Goebel,  J.  Perel,  J.  P.  Mahoney,  and  A.  T.  Forrester,  "Characteristics 
of  the  Capillaritron  Ion  Source",  submitted  for  publication  in  the  Journal 
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III.  ENERGY  TRANSFER  AND  RELAXATION  IN  SMALL  POLYATOMIC  MOLECULES 

A.  PHOTOCHEMICAL  DYNAMICS  OF  79Br2  IN  Xe  MATRICES* 

(Paul  Beeken,  Mary  Mandich,  George  Flynn) 

(JSEP  work  unit  6,  1979-82) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

The  photochemical  dynamics  of  Br2  trapped  in  a  Xe  matrix  have  been 

3  +  1 

investigated.  Following  excitation  into  either  the  B  it  (0y)  or  *(lu)  states, 

emission  occurs  in  the  far  red  and  infrared  region.  Analysis  of  the  emission 

3  +  1  + 

allows  the  emitting  states  to  be  characterized  as  B  tf(0u)  v’«0  ■*  X  E^,  and 

3  1  + 

(tentatively)  the  A'  tr(2u)  -*•  X  E  .  The  spectrum  of  the  B  -*•  X  emission  was 
recorded  and  exhibited  unusually  broad  vibronic  emission  lines  when  compared 
with  Ar  spectra.  The  B  -*■  X  emission  lifetime  was  measured  to  be  3.6  +  .5  usee 
and  is  in  excellent  agreement  with  the  expected  value  for  the  radiative  life¬ 
time.  The  risetime  of  the  B  signal  was  found  to  be  limited  by  the  detection 

system  (y  15  nsec).  The  tentatively  identified  A'  X  emission  has  a  lifetime 

3  1  + 

of  4  +  1  msec.  The  A  ir(lu)  -+■  X  E^  emission  resolved  clearly  in  Ar  and  Kr 
matrices  is  altogether  absent  from  Br^  isolated  in  Xe.  Experiments  were  per¬ 
formed  to  measure  the  relative  intensities  of  the  B/A’  emissions  in  Xe  and 
compare  them  with  earlier  results  from  Ar  matrices.  The  Xe  B/A'  ratio  was 
found  to  be  nearly  100  times  greater  in  Xe  then  in  Ar.  The  observations  were 
found  to  be  independent  of  temperature  from  12K  to  60K  and  independent  of  con¬ 
centrations  from  1:500  to  1:10000,  Br^  in  Xe. 

This  work  has  been  accepted  for  publication  in  the  Journal  of  Chemical 
Physics.  Please  refer  to  the  article  (June  1982)  for  a  more  detailed  account. 

*  This  research  was  also  supported  by  the  Department  of  Energy  under 
Contract  DE-AS02-78ER04940  and  the  National  Science  Foundation  under  Grant 
NSF-CHE80-23747. 


B.  EMISSION  SPECTRA  AND  RELAXATION  DYNAMICS  OF  EXCITED  79Br  In  Ar 
AND  Kr  MATRICES*  Z 

(Mary  Mandich,  Paul  Beeken,  George  Flynn) 

(JSEP  work  unit  6,  1979-82) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

The  photochemical  dynamics  of  Br^  trapped  in  Ar  and  Kr  matrices  have 

3  +  1 

been  studied.  Following  laser  excitation  into  either  the  B  11(0^)  or  II(lu) 

states,  emission  occurs  in  the  near  infrared  spectral  region.  Analysis  of  this 

emission  Indicates  that  it  consists  of  three  electronic  transitions  which 

have  been  assigned  as  B9II(o+)  v'*0  -*•  X^E+  (previously  reported),  A9H(lu) 

u  8 

vf  «0  +  X3!*  and  (tentatively)  A,9H(2u)  -*■  X^T*.  The  B  -*■  X  emission  lifetime 

was  measured  to  be  8,0  +  .5  usee  in  Ar  matrices  (5.3  +  .5  usee  in  Kr)  and  is 

in  good  agreement  with  both  earlier  determinations  in  matrices  and  the  gas 

phase  B  X  radiative  lifetime.  The  observed  A  -*■  X  emission  lifetime  of 

170  +  10  usee  in  Kr  is  approximately  the  same  as  the  measured  A  X  gas  phase 

radiative  lifetime  but  appears  shortened  in  Ar  matrices  where  the  lifetime 

is  67  +  4  usee.  The  lifetime  of  the  A'  -*■  X  transition  is  11  +  1  msec  in  Ar 

and  6+1  msec  in  Kr  matrices.  The  spectral  and  temporal  data  are  consistent 

with  an  overall  energy  transfer  mechanism  involving  both  vibrational  and 

electronic  relaxation.  The  intramolecular  electronic  relaxation  is  somewhat 

3  + 

restricted  in  that  the  highest  lying  B  II (0^)  state  is  populated  only  by  direct 

absorption  into  its  bound  and  continuum  levels  whereas  the  two  lower  states, 

3  3  13  + 

A  H(lu)  and  A*  !I(2u),  acquire  population  after  either  the  TI(lu)  or  B  n(O^) 

state  are  initially  excited.  Both  the  initial  vibrational  and  electronic  relax. 

tlon  rapidly  (<  20  nsec)  routes  the  excited  population  into  the  lowest  v '  -  0 

levels  of  the  A  and  B  states  (and  tentatively  A')  where  the  subsequent  decay  occurs 
*•6  -2 

on  a  10  -  10  sec  timescale.  All  of  the  observed  emission  lifetimes  are  in- 

depedent  of  temperature,  laser  power  and  Br2:  matrix  dilutions  ranging  from 
1:500  -  1:10,000. 
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This  work  has  been  accepted  for  publication  in  the  Journal  of  Chemical 


Physics.  Please  refer  to  the  article  (July  1982)  for  a  more  detailed  account. 


*This  research  was  also  supported  by  the  Department  of  Energy  under 
Contract  DE-AS02-78ER04940  and  the  National  Science  Foundation  under  Grant 
NSF-CHE80-23747. 
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C.  PHOTOCHEMICAL  DYNAMICS  OF  I2  IN  RARE  GAS  MATRICES* 

(Paul  Beeken,  George  Flynn) 

(JSEP  work  unit  6,  1979-82) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

I^  Is  a  fascinating  molecular  system  which  has  stirred  considerable 
Interest  because  of  the  development  of  the  iodine  atom  laser. ^  ^  This 

device  has  shown  promise  as  a  high  powered  source  for  laser  fusion  experi- 

(5)  * 

ments.  One  possible  drawback  to  the  attainment  of  a  high  efficiency  I 

chemical  laser,  however,  might  be  interference  from  the  quenching  of  excited 

iodine  atoms  by  I^.  A  complete  analysis  of  the  intramolecular  relaxation 

* 

pathways  for  molecular  iodine  is  thus  essential  for  optimizing  the  I  laser. 

We  have  just  completed  a  detailed  analysis  of  the  spectra  and  dynamics 
of  Br2  in  rare  gas  matrices. ^ ^  Other  workers  have  investigated  the 
dynamics  of  Cl2  in  low  temperature  solids. v  '  The  common  elements  of  these 
studies  has  led  to  a  generalized  picture  of  halogen  molecule  relaxation  in 
a  rare  gas  case. ^  The  investigations  into  Br^  dynamics  include  characteriza¬ 
tion  of  the  lower  electronic  states. ^  In  iodine  the  low  lying  A  and  A' 

*  (9) 

electronic  states  are  thought  to  hinder  efficient  operation  of  the  I  laser. 
Studies  of  I^  would  not  only  serve  to  solidify  some  of  the  generalizations 
proposed  for  intramolecular  halogen  relaxation  but  would  also  contribute  to 
the  development  of  a  laser  system. 

Experimental 

The  apparatus  has  been  described  elsewhere  and  will  not  be  presented 
here.  The  only  change  in  the  experimental  setup  has  been  the  inclusion 

of  a  new  Lambda  Physik  excimer  pumped  dye  laser  system. 
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Results 


When  a  matrix  of  1:1000  I^Ar  was  deposited  on  a  Csl  substrate  and 
irradiated  with  568.2  nm  light  from  a  pulsed  dye  laser  two  fluorescences 
were  observed.  The  signals  occur  in  two  different  spectral  regions.  One 
signal  (visible  to  far  red  region  of  the  spectrum)  was  found  to  decay  ex¬ 
tremely  fast  (<  200  usee).  The  range  of  this  signal  extends  from  below 
700  nm  and  extends  to  just  beyond  1100  nm.  Using  an  SI  type  phototube  and, 
in  another  experiment,  a  more  sensitive  GaAs  phototube,  we  were  unable  to 
resolve  this  emission  into  distinct  spectral  features. 

The  other  fluorescence  signal  occurs  in  the  IR  and  was  found  to  have 
a  lifetime  of  ^250  usee.  The  spectral  range  of  this  signal  extended  from 
1.2  u  to  1.6  u.  Work  is  currently  in  progress  to.  resolve  this  emission 
using  a  monochromator. 


Discussion 


The  emission  spectrum  for  the  X1Eg  -*•  B3II(0+)  transition  has  been 

(10) 


reported  for  isolated  in  a  rare  gas  matrix.  The  experiment  employed 

+ 


a  cw  Ar  ion  laser  for  excitation  so  no  lifetimes  were  measured.  Dimer 
species  were  also  detected  which  complicated  the  otherwise  simple  emission 
spectrum.  The  spectral  bands  were  identified  as  emissions  from  the  lowest 


vibrational  level  of  the  excited  B  electronic  state.  This  observation  is 

(6)— (8) 


consistent  with  what  has  been  observed  in  matrix  isolated  Br^  and  Cl^. 


Measurments  of  the  gas  phase  lifetime  of  the  B  state  only  allow  us  to 


estimate  the  radiative  lifetime  since  the  v’«0  to  ground  electronic  state 


transition  has  not  been  characterized. 


(ID 


The  value  obtained  by  extrapola¬ 


tion  from  higher  vlbronlc  transitions  would  be  •n*  800  nsec. 


(ID 


If  the  matrix 


index  of  refraction  Is  used  to  correct  this  value  then  the  expected  radiative 
llfetlae  In  Ar  would  be  -\-650  usee. ^ ^  In  the  gas  phase,  the  rate  of  pre- 
dissociation  is  known  to  be  quite  rapid  as  excitation  moves  down  the  B  vibra¬ 
tional  ladder  from  v'  -  10.^^  Therefore,  as  was  suggested  for  Br2  in  various 
rare  gas  matrices,  non-radiative  curve  crossing  could  shorten  this  lifetime 
considerably. ^ 

The  IR  emission  has  not  been  spectrally  resolved  using  a  monochromator 
and  an  InSb  detector  yet.  But  the  general  spectral  features  of  the  emission 
were  determined  with  a  series  of  narrow  band  interference  filters.  The  resolu¬ 
tion  of  these  filters  is,  unf ortunately ,  not  sufficient  to  determine  the  origin 
of  the  emitting  species  by  spectroscopic  means.  A  crude  estimate  of  the  radia¬ 
tive  lifetime  for  the  A  state  places  the  expected  value  in  the  range  of  several 

(12) 

hundred  microseconds.  Since  the  lifetime  of  the  A'  state  is  expected 

to  be  much  longer  than  this  based  on  observations  of  other  halogen  species, 

3 

the  most  likely  source  for  the  250  usee  emission  would  be  the  A  n(lu)  electronic 
stage.  Further  experiments  are  required,  however,  to  confirm  these  conclusions. 

*This  research  was  also  supported  by  the  National  Science  Foundation  under 
Grant  NSF-CHE80-237A7  and  the  Department  of  Energy  under  Contract 
DE-AS02-78ER0A9A0, 
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D.  LASER  EXCITED  INFRARED  FLUORESCENCE  IN  OXALYL  FLUORIDE:  RELAXATION 

IN  THE  PRESENCE  OF  A  LOW  ENERGY,  ONE  DIMENSIONAL  QUASI-CONTINUUM* 

(Toomas  Allik,  George  Flynn) 

(JSEP  work  units  4  and  5,  1979-82) 

(Principal  Investigator:  G.  W,  Flynn  (212)  280-4162) 

Work  on  this  project  was  finished  in  January  of  1982.  A  paper  describing 
these  results  is  scheduled  for  publication  in  the  Journal  of  Physical  Chemistry, 
August  1982.  An  abstract  of  the  paper  is  given  below: 

Fluorescence  has  been  observed  from  various  fundamental  and  combination 
states  of  oxalyl  fluoride  excited  by  a  Q-switched  CO^  laser.  The  presence  of 
a  low  frequency  torsional  mode  in  this  molecule  provides  a  one  dimensional 
quasicontinuum  of  vibrational  levels  at  energies  above  the  barrier  to  rotation 
about  the  C-C  bond.  Intermode  vibrational  relaxation  rates  in  (COF)^  are 
extremely  rapid  (possibly  collisionless).  The  observed  temperature  dependence 
of  fluorescence  intensity  and  the  overall  rapid  equilibration  in  (COF) ^  suggest 
a  mechanism  in  which  many  rovibrational  transitions  that  possess  character 
(e.g.  0  -*■  -*■  2v1Q,  vx  vx  +  viq  etc.)  are  excited  by  the  laser. 

*This  research  was  also  supported  by  the  Department  of  Energy  under 
Contract  DE-AS02- 78ER04940  and  the  National  Science  Foundation  under  Grant 
NSF-CHE80-23747. 
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E.  INVESTIGATION  OF  CYCLIC  ENERGY  TRANSFER  PATHWAY  IN  OCS/He* 

(Steven  Goates,  George  Flynn) 

(JSEP  work  units  4  ant1  5,  1979-82) 

(Principal  Investigator:  G,  W.  Flynn  (212)  280-4162) 

Previously  in  this  laboratory  the  rates  of  collisional  vibrational 
energy  transfer  between  the  level  of  OCS  and  the  bending  mode  as  well 
as  the  overall  vibrational  relaxation  to  the  translational/rotational  degrees 
of  freedom  were  measured  for  mixtures  of  OCS  and  the  various  rare  gases. 

These  measurements  were  made  by  observing  fluorescence  from  the  state 

4 

following  absorption  of  a  laser  pulse  by  2\>2.  In  the  OCS/  He  system  the 
equilibration  of  the  and  v 2  modes  was  found  to  be  unexpectedly  rapid  with 

respect  to  the  rates  for  the  other  OCS/rare  gas  mixtures.  For  all  mixtures 
except  OCS/^He,  the  dominant  energy  transfer  pathway  connecting  the  and 
modes  is 

OCS(2v2)  +  M  ;— ■>  OCS(v^)  +  M  +  1 88  cm”1  (1) 

along  with  equilibration  of  \>2  and  2^  via 

OCS(2v2)  +  0CS(0)  < — ^  2  0CS(v2)  +  M  +  7  cm  1,  (2) 

the  intramode  equilibration  (process  2)  occurring  much  more  rapidly  than 
the  interoode  transfer  (process  1) .  It  was  proposed  that  in  the  case  of 
OCS  and  He  a  second  pathway,  namely 

0CS(v2)  +  M  OCS(v^)  +  M  -  339  cm"1,  (3) 

becomes  competitive  with  the  first,  accounting  for  the  rapidity  of  the  equilibra¬ 
tion.  It  is  to  be  noted  that  the  two  pathways  together  constitute  a  "cyclic" 
pathway.  Such  cyclic  routes  have  been  postulated  to  play  an  important  role  in 
the  relaxation  of  molecules  but  have  been  rarely  verified  experimentally. 

If  our  suspicions  are  well-founded,  that  the  pathway  directly  coupling 
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4 

and  is  important  for  collisions  of  OCS  with  the  light  He  atom,  this 

3 

route  should  be  further  enhanced  for  the  yet  lighter  collision  partner  He. 

3 

Accordingly,  rise  and  fall  rates  of  fluorescence  for  OCS/  He  mixtures 
have  been  measured.  While  the  analysis  of  these  data  have  not  yet  been 
completed,  it  is  clear  that  the  equilibration  of  and  in  this  mixture 

4 

proceeds  significantly  faster  than  in  OCS/  He.  This  result  is  certainly 
indicative  of  the  proposed  cyclic  pathway,  but  does  not  provide  definitive 
evidence. 

An  analysis  of  the  kinetics  for  the  cyclic  pathway  reveals,  however, 

that  its  presence  can  be  clearly  distinguished  by  observing  fluorescence 

(2) 

from  the  state  at  19.1  p.  Two  exponential  decays  of  this  fluorescence 

have  been  observed  in  pure  OCS  and  are  related  to  the  rates  of  v^/v^  equilibra- 

(3) 

tion  and  vibrational  to  translational/rotational  relaxation.  The  relative 

amplitudes  of  these  two  decay  rates  can  reveal  the  presence  of  the  cyclic 

pathway.  An  attempt  to  sufficiently  accurately  determine  the  amplitudes  from 

4  3 

the  weak  19p  fluorescence  for  OCS/  He  and  OCS/  He  mixtures  awaits  the  arrival 
of  a  detector  sensitive  to  this  long  wavelength, 

*This  research  w.  also  supported  by  the  Department  of  Energy  under 
Contract  DE-AS02-78ER04940  and  the  National  Science  Foundation  under  Grant 
NSF-CHE80-23747. 

(1)  M.  L.  Mandich  and  G.  W.  Flynn,  J.  Chem.  Phys.  73,  1265  (1980). 

(2)  M.  L.  Mandich  and  P,  Beeken,  unpublished  research. 

(3)  M.  L.  Mandich  and  G.  W.  Flynn,  J.  Chem,  Phys.  73,  3679  (1980). 
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F.  VIBRATIONAL  RELAXATION  DYNAMICS  IN  POLYATOMIC  GAS  MIXTURES: 

A  STUDY  OF  COF2  AND  NO* 

(J.  Subbiah,  G.  W.  Flynn) 

(JSEP  work  units  4  and  5,  1979-1982) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

In  order  to  achieve  control  of  energy  flow,  storage  and  release  in 
gas  phase  laser  systems,  electric  discharges,  shock  waves  and  gas  dynamic  flows,  it 
is  important  to  understand  the  dynamics  of  energy  flow  among  the  internal  degrees 
of  freedom  of  molecules.  As  of  now  an  empirical  study  of  the  experimental  data 
has  led  to  some  understanding  of  qualitative  trends.  In  addition,  semi-quanti¬ 
tative  theories  can  explain  the  experimental  data,  but  accurate  a  priori  pre¬ 
dictions  are  still  not  feasible.  This  study  of  intermolecular  vibrational 
energy  transfer  between  C0F2  and  NO  is  aimed  at  improving  our  understanding  of 
relaxation  dynamics. 

2 

NO  is  an  interesting  case  because  its  ground  electronic  state  (  it)  has 
an  unpaired  spin  which  may  result  in  unusual  relaxation  behavior.  In  fact,  its 
vibration-vibration  (V-V)  relaxation  is  of  the  same  order  of  magnitude  as  CO 
(which  is  similar  in  many  respects) ,  while  its  vibration-translation  (V-T) 
relaxation  rate  is  about  five  orders  of  magnitude  faster  than  CO.^  Most  of 
its  other  V-V  behaviour  is  normal,  but  its  excitation  from  laser-excited  C0F2 
may  still  hold  some  surprises. 

Fig.  1  shows  the  relevant  energy  levels  of  C0F2  and  NO.  The  energy 
flow  pattern  in  vibrationally  excited  neat  C0F2  has  been  studied. ^  There  it 
has  been  shown  that  following  excitation  of  the  v2  mode  by  different  C02  laser 
lines  (10. 6 p  P  branch),  C0F2  undergoes  intra  and  intermode  relaxation  as  shown 
below: 


Type  1 

excitation: 

cof2(0)  - 

- ►  COF2(v2) 

hv 

(1) 

Type  2 

excitation: 

COF2(v2)  ■ 

- 1  COF2(2v2) 

(2) 
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v-1 


2',6*“ 

LASER 
TYPE  2 
EXCITATION 


u6 


LASER 
TYPE  1 
EXCITATION 


C0F2  NO 

FIG.  1.:  ENERGY  LEVELS  OF  COF2  AND  NO. 
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(3) 


k 

u 

Ladder  Climbing:  COF2(v2>  +  COF2(v2)?=i  COF2(0)  +  COF2(2v  ) 

kd  " 

Equilibration  :  COF2(2v2>  +  M  ^  COF^V^  +  M  (4) 

of  Fermi-coupled  "  * 

levels 

Fig.  2  shows  typical  type  1  and  type  2  fluorescence  traces  in  COF2  alone.  (The 
same  eigenvalues  appear  in  both  types,  but  the  fastest  eigenvalue,  which  appears 
in  the  rise  of  type  1  turns  up  in  the  fall  of  the  type  2  signal.)  Subsequently 
the  excitation  remains  localized  in  the  v^,  v2  manifolds  for  about  35  psec  at 
1  torr  before  further  relaxation  sets  in. 

If  NO  is  added  to  this  system  one  possible  new  path  would  include  energy- 
transfer: 

Crossover:  COF2(v^)  +  NO  — ►  COF2(0)  +  NO  (v  ■  1)  (5) 

The  crossover  rate  in  COF2~NO  mixtures  was  determined  by  following  the  time- 
dependent  population  in  COF^v^  2v2)  via  its  5p  fluorescence  intensity. 


Experimental  and  Results 


A  Q-switched  C02  laser  operating  on  different  lines  of  the  10.6  micron 
P  branch  was  used  to  excite  the  v2  mode  of  C0F2  (C-F  stretch) .  Typical  pulse 
energies  were  0.5  to  1  mj /pulse,  and  typical  pulsewldths  about  1  microsecond. 

The  laser  was  run  single  line  to  better  than  95%.  The  C0F2  used  was  taken 
directly  from  a  high  purity  Synthatron  tank,  while  the  NO  was  purified  by  succes¬ 
sive  trap-to-trap  distillation  at  -160°  C.  In  all  cases  IR  fluorescence  was 
observed  through  a  4.6  micron  long  pass  filter  as  well  as  a  2  mm  thick  MGF2 
window.  This  isolates  the  5  micron  fluorescence  of  C0F2  corresponding  to  the 
»1»  2v2  band  (C  ■  0  symmetric  stretch),  as  veil  as  the  v  ■  1  fluorescence  from 
NO.  This  does  not  complicate  the  interpretation,  however,  as  the  same  eigen¬ 
values  appear  in  the  NO  and  C0F2  fluorescences,  although  with  different  amplitudes. 
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TYPICAL  5«  FLUORESCENCE  DECAY 
OF  cof2  ALONE  FOLLOWING  TYPE  1 
EXCITATION. 

COF2  PRESSURE  1.023  TORR 
O.N  NS  FULL-SCALE. 


TYPICAL  5»  FLUORESCENCE  DECAY 
OF  COF2  ALONE  FOLLOWING  TYPE  2 
EXCITATION. 

CQF2  PRESSURE  1.006  TORR 
0.02  TS  FULL-SCALE 
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Fluorescence  decay  rates  of  the  COF^  (v,  ,  2 )  level  were  studied  as 
a  function  of  pressure  at  different  mole  fractions  of  NO  in  the  COF^-NO  mix¬ 
tures.  The  results  for  both  type  1  and  type  2  excitation  are  shown  in  Figure  3. 
The  two  types  give  eigenvalues  that  agree  within  experimental  error.  Further, 
in  order  to  see  whether  NO  is  merely  acting  as  an  inert  collision  partner,  data 
taken  with  COF^-Ar  mixtures  at  identical  mole  fractions  is  also  tabulated  in 
Table  2. 

Fig.  3  shows  the  simplified  kinetic  scheme  used  to  model  the  intermolecular 
energy  transfer.  Under  the  approximations  listed,  there  are  two  eigenvalues, 
and  one  can,  in  this  case,  express  each  eigenvalue  in  terms  of  the  rate  constants. 
The  fast  eigenvalue  is  linear  both  in  the  total  pressure  of  the  mixture  at  con¬ 
stant  mole  fraction  as  well  as  in  the  pressure  of  NO  added  to  a  fixed  pressure 
of  COF-.  Although  k^  can  be  obtained  from  any  value  of  the  mole  fraction,  the 
low  mole  fraction  runs  are  not  as  accurate  as  the  here  turns  out  to  be  the 
difference  of  two  large  numbers  (see  Equation  (4),  Fig.  3).  Table  1  lists  the 
calculated  eigenvalues  based  on  this  scheme  while  Table  2  includes  the  experimen¬ 
tally  observed  rates.  As  can  be  seen  the  rates  converge  to  a  value  around 
170  ms  .^torr-'*'  for  k^. 

Discussion 

The  measured  k^  (the  rate  constant  for  energy  transfer,  see  Equation  (5)) 

corresponds  to  50  collisions  or  a  vibrational  energy  transfer  probability  per 
_2 

collision  of  2  x  10  at  room  temperature,  which  is  a  large  value.  This  is  about 
an  order  of  magnitude  larger  than  expected  even  for  exact  resonance  on  the  basis 
of  a  theory  of  energy  transfer  assuming  exponentially  repulsive  interaction 

(3) 

potentials.  It  seems  likely,  however,  that  long-range  attractive  interactions 

(4) 

involving  dipole-dipole  potential  (as  modelled  by  Sharman  and  Brau  )  should 
in  agreement  with  the  results.  Work  along  these  lines  is  being  pursued. 
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COF2  NO 


C0F2d/2)  +  C0F2(v2)5:=tC0F2(2^)  +  C0F2(0)+14cm-l  (i  ) 

COFgt^.a^J+NOtvOJ^siCOFgtO+NOtvD+SScm-1  (Li) 

Kb 


THE  RESULTING  RATE  MATRIX  A  IS 


4rCOF2(v,,2l/2)  -Kd(l-x)-Kf  *-X  Kb(1-x)  NCOp2 

a'  -  ‘•v**1  mn 


#NO(«»l) 

a  1 


-KJ  l-x)-X  NnO 

JL  (v»u 


RESULTS:  • 

*  0.99  (1^+  Kb><  1"»)  +  Kf  * 

X"  *  kf  0.99  kb  ftor2+  0.99 1^  PCOf2 

X"  «  [0.04  kf  RCOf2  pno  +  kb  Pno  1 

ffcorj-PARTfALPRESSURE  OF  COF2 
Pnq  'PARTIAL  PRESSURE  OF  NO 


FIG.  3  KINETIC  SCHEME  FOR  COF2  -  NO  ENERGY 
TRANSFER 
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TABLE  1:  Numerical  Solution  for  the  Time-dependence  of  Population  in 
2^)  using  the  kinetic  scheme  in  Fig.  4. 

Following  is  the  calculated  time-dependent  variation  of  the  population 
deviation  (from  thermal  equilibrium)  of  C0F2(v^,  2vj)  in  a  mixture  of 
CO2  and  NO  at  an  NO  mole  fraction  of  0.9. 

The  guess  values  for  the  bimolecular  rate  constants  used  are:  (in  units 
of  ms  torr  ) 


C0F2(vif 


k  *  620 
u  . 

kd  -  620 

kf  -  170 

k20  ”  30 

■  4-5 

k10  "  2*5 

klo  *  2'° 

2v  )  (t)  *  A  e  *1C  +  A  e  +  A  e  +-  B 
2  12  3 


0.93 

X  - 

242  ms" 

1 

0.04 

X2  - 

13.1  ms 

0.01 

A3“ 

5.4 
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TABLE  2s  Results  of  Type  1  and  Type  2  excitation  of  CC^-NO  mixtures 


Type  2  Excitation 

Slope  of  *£ast  (®s  S  versus  pressure  of  added  NO:  154  ms-*  torr 


Gases  Mix 

cof2-m 

Mole  fraction 

of  M 

^fast 

(experimental) 

(ms-1) 

kf 

ms-1  torr- 

cof2-no 

0.986 

164 

171  to  137 

cof2-no 

0.9 

240 

195  to  136 

cof2-no 

0.7 

317 

161  to  71 

C0F,-Ar 

0.9 

95 

46  to  16 

* 

-ikf  • 
ms  torr 

155 

170 


150 


Type  1  excitation 


The  slope  of  ^fagt  (ms-1}  versus  pressure  of  added  NO:  170  ms-1  torr-1 


‘The  range  of  k^  got  using  the  experimental  lfa8t  and  eqn.  iv  of  Fig-  3. 

The  guess  value  of  which  when  used  in  numerical  solution  using  this 
kinetic  scheme  results  in  the  best  agreement  with  *fast* 
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The  expected  additional  rate  process  following  addition  of  NO  to  COF^  is 
not  experimentally  observed,  but  this  is  understandable  in  view  of  the  small 
amplitude  expected  for  it  based  on  our  kinetic  model  (the  rate  processes  included 
are  indicated  in  Fig.  4  and  the  results  shown  in  Table  1) .  The  rates  with  argon 
are  found  to  be  non-zero  but  still  much  smaller  than  the  rates  observed  with  NO. 

Table  3  lists  the  known  transfer  probabilities  from  NO  to  other  molecules. 

It  is  seen  that  the  NO  to  D^S  probability  is  very  large,  as  might  be  expected 
for  its  small  energy  defect.  The  NO  to  N02  is  also  large  and  one  possible  explana¬ 
tion  has  been  a  possible  reactive  channel: 

NO  +  N02  - *  N203  +  9.7  Kcal/mole 

There  is  the  possibility  that  a  similar  channel  may  be  available  for.  the  COF2~ 

NO  system  as  well.  The  C  atom  in  COF2>  flanked  by  two  highly  electronegative 
F  atoms,  as  well  as  the  polarizable  cloud  connected  to  the  electronegative  0 
atom,  is  likely  to  be  an  electron  deficient  center.  NO,  on  the  other  hand,  has 
an  odd  electron,  and  its  attack  on  the  nucleophilic  C  atom  seems  plausible. 

In  the  case  of  the  HF-COF2  and  DF-C02  systems,  where  similar  strong  intermolecular 
forces  operate,  it  was  found  that  agreement  with  experiment  could  be  improved 
by  invoking  the  possiblity  of  'orbiting'  collisions.  ^  That  is,  when  the 
kinetic  energy  is  less  than  0.8  times  the  well  depth,  for  certain  values  of  the 
impact  parameter  the  trajectory  can  no  longer  be  described  by  a  straight  line 
path,  but  executes  multiple  revolutions  (i.e.  scattering  angle  greater  than  2  ir ) . 
Inclusion  of  such  trajectories  in  the  case  of  the  HF-C02  system  resulted  in  in¬ 
creasing  the  calculated  V-V  probability  by  a  factor  of  2.5.  Another  system 
that's  probably  reactive  is  C02-H20,^  with  a  probability  of  energy  transfer 
of  5  x  10  2. 

A  mixture  of  C0F2  and  NO  in  a  stainless  steel  cell  develops  a  new  IR 
absorption  band  which  grows  as  a  function  of  time.  The  pressure  dependence 
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kf 


kf 

C0F2(i/i  )  +  N0(0)  5==*  C0F2(0)  +N0  (v*1) 


COF2(i/2)  +  NO  — C0F2(0)  +  NO 
k?n 

C0F2(i/2)  +  C0F2-=i  C0F,(0)  +  COF2 


NO (v»1)  4-  NO  — ^  NO  (0)4' NO 


It 

NO  (v *1)  +  C0F2—^-  NO  (0)  +  COF2 

THIS  RESULTS  IN  A  SYSTEM  OF  THREE  INDEPENDENT 
LINEAR  FIRST  ORDER  COUPLED  DIFFERENTIAL 
EQUATIONS. 


FIG.  4  DETAILED  KINETIC  SCHEME  FOR  THE  COF2-NC 
SYSTEM  INCLUDING  V-T  RELAXATION 
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TABLE  3:  Interne  lecular  Vibrational  Energy  Transfer  Probabilities. 
(Largely  taken  from  Ref.  5a) 


Prom  M1 

To  M2 

AE*  cm-1 

B+ 

NO 

CO 

267 

10*4 

10,000 

NO 

n2 

455 

6  x  10"6 
-4 

166,667 

NO 

co2 

491 

1.7  x  10 

5,882 

NO 

d2s 

-16 

1.1  x  10'2 

91 

NO 

h2o 

281 

6  x  10~3 

167 

NO 

h2s 

586 

3.2  x  10“3 

313 

NO 

d2o 

697 

10~3 

1,000 

NO 

ch4 

343 

9  x  10'4 

1,111 

NO 

cof2 

-68 

2  x  10'2 

50 

HP 

C02(00°1) 

1612 

8.5  x  10-3 
-2 

118 

DF 

C02(00"1) 

557 

3  x  10 

33 

co2 

H20 

55 

5  x  10“2 

20 

8  (collision  number) 


where  Og^  is  the  gas-kinetic  collision  cross- 


section,  and  is  the  vib  energy  transfer  cross  section. 
*Energy  defect  AE:  M^(v»l)  +  M2  -*■  M^(0)  +  M2(v«l)  +  AE 


**. 


Probability  of  vibrational  energy  transfer  per  collision. 
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of  the  intensity  of  this  band  suggests  a  1:1  complex  formation  or  reaction. 
However  the  reaction  appears  to  be  slow  enough  that  on  the  timescale  of  a 
fluorescence  measurement  the  fraction  of  molecules  that  have  reacted  is  not 
large.  It  would  be  interesting  to  study  the  temperature  dependence  of  the 
intensity  of  the  1800  cm  ^  band.  It  would  be  expected  that  in  the  case  of 
complex  formation  the  equilibrium  would  be  shifted  in  the  direction  of  the 
complex  at  low  temperature.  Once  that  is  known  the  temperture  dependence  of 
the  crossover  probability  may  yield  information  about  the  role  of  the  complex 
in  the  excitation  transfer. 

Ultimately  the  most  definitive  evidence  for  energy  transfer  would  be 
to  directly  monitor  the  population  of  NO  in  v  *  1.  This  however,  is  hard  to 
do  by  IR  fluorescence,  both  because  of  the  relatively  low  intensity  of  the 
NO  fundamental,  and  because  its  rotational  envelope  overlaps  that  of  COF^- 
One  alternative  would  be  to  directly  excite  NO  by  a  CO  or  doubled  CO^  laser, 
and  monitor  COF^  5  micron  fluorescence. 

This  system  is  interesting  in  the  context  of  using  C02  as  an  energy  storage 
medium  or  sensitizing  agent.  For  example,  the  relaxation  time  of  the  population 
in  NO  (v  *  1)  in  the  gas  mixture  is  about  77  usee  at  1  torr  total  pressure 
(NO  mole  fraction  *  0.9).  It  may  now  be  possible  to  use  this  metastable  popula¬ 
tion  to  study  the  vibrational  enhancement  of  the  chemical  reactivity  of  NO  with, 
say,  Cl^.  In  addition,  this  is  an  example  where  a  population  of  vibrationally 
excited  molecules  can  be  prepared  even  when  no  coincidence  with  a  convenient 
laser  can  be  found. 

The  intermolecular  vibrational  relaxation  probability  from  C0F2  to  NO 

-2 

has  been  measured  to  be  2  x  10  .  This  can  be  understood  in  terms  of  long-range 

attractive  Interactions  and  possibly  reactive  scattering.  This  system  has 
interesting  applications  Involving  sensitization  and  energy  storage. 
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G.  OBSERVATION  OF  VIBRATIONAL  EXCITATION  IN  PHOTOFRAGMENTS  OF  S-TRIAZINE* 

(Steven  Goates,  Jack  Chu,  George  Flynn) 

(JSEP  work  unit  7,  1979-82) 

(Principal  Investigator:  G.  W.  Flynn  (212)  280-4162) 

Last  year's  progress  report  contained  very  preliminary  results  of  an 
experiment  which  seeks  to  determine  excess  energy  in  the  vibrations  of  HCN 
molecules  produced  in  the  photodissociation  of  s-triazine,  (HCN) ^ .  Work  on 
this  project  has  recently  been  advanced. 

s-Triazine  is  a  benzene-like  molecule  which  has  elicited  much  attention 

in  recent  years  becuase  of  the  importance  of  Jahn-Teller  distortion  of  its 

excited  states. ^  Molecular  beam  experiments  have  demonstrated  that  triazine 

is  cleanly  dissociated  into  three  HCN  fragments  when  exposed  to  193  nm  radia- 
(2) 

tion.  The  distribution  of  excess  energy  among  the  degrees  of  freedom  of 
the  HCN  fragments  is  expected  to  provide  information  about  the  dissociating 
states  of  the  triazine  molecule  and  the  dissociation  process.  Considering 
the  different  geometries  of  the  HCN  units  in  triazine  and  free  HCN  molecules 
(bent  versus  linear)  a  high  degree  of  vibrational  excitation  in  the  HCN  pro¬ 
ducts  is  not  unexpected. 

We  have  used  the  pulsed  193  nm  output  (ArF)  from  a  rare  gas-halide 
excimer  laser  to  irradiate  a  flowing  sample  of  triazine  in  a  glass  cell  and 
viewed  Infrared  fluorescence  following  each  pulse  at  a  right  angle  to  the 
UV  beam  path.  Combinations  of  filters  were  employed  to  select  particular 
wavelength  regions  and  the  infrared  fluorescence  was  detected  with  either  a 
HgCdTe  or  InSb  detector;  temporal  response  of  the  detectors  was  recorded  with 
a  transient  signal  digitizer /averager  combination.  We  have  made  crude  mea¬ 
surements  of  the  absorption  coefficient  of  triazine  at  193  nm  ('v  0.025  torr  ^cm  *) 
and  calculated  that  less  than  a  few  hundredths  of  a  percent  of  the  triazine 
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\  in  Che  beam  path  are  dissociated  with  each  pulse  of  the  laser  at  the 

2 

fluences  20  mJ/cm  )  employed.  Confirmation  of  the  low  dissociation 
provided  by  the  very  small  increase  in  pressure  observed  after  irradia- 
riazine  vapor  in  a  closed  cell. 

espite  the  low  dissociation  yields  we  have  been  able  to  observe  infrared 

ence  from  many  vibrational  transitions  in  HCN  3v^  and  v^, 

v^  +  v^,  2Vj),  and  the  HCN  fragments  appear  to  be  highly  vibrationally 

The  rise  of  fluorescence  from  all  HCN  states  is  essentially  instan- 

and  preliminary  analysis  indicates  that  excess  vibrational  energy 

y  evenly  distributed  among  all  the  modes.  Time-of -flight  measurements 

tolecular  beam  study  of  reference  2  revealed  that  very  little  of  the 

tcess  energy  resides  in  the  translational  motion  of  the  fragments.  Rise 

:  fluorescence  which  we  have  observed  from  undissociated  triazine,  excited 

Lsions  with  the  HCN  fragments,  seem  to  match  the  molecular  beam  results. 

to  quantitatively  determine  the  amount  of  energy  residing  in  each  mode 

ition  of  HCN  immediately  after  triazine  dissociation  by  employing  the 
(3) 

s-f liter  technique  for  determining  vibrational  temperatures.  The 

3f  rotational  excitation  cannot  be  established  from  our  present  mea¬ 
ts;  experiments  with  a  circularly  variable  filter  are  planned  to  probe 
itional  envelope  of  the  observed  transitions. 

)ur  results  for  193  nm  dissociation  of  s-triazine  contrast  sharply  to 

rom  a  study  of  dissociation  of  a  similar  molecule,  s-tetrazine,  where 

(4) 

vibrational  excitation  of  the  fragments  (2  HCN  +  N^)  was  found.  These 
suggest  a  much  different  dissociation  mechanism  for  these  two  molecules 
their  structural  similarity.  Molecular  beam  work  cited  above  indicates, 

,  that  a  much  different  distribution  of  excess  energy  occurs  in  the 
ition  fragments  of  triazine  when  it  is  irradiated  with  248  nm  rather 
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(2) 

than  193  nm.  We  will  extend  our  studies  to  248  nm  dissociation  in  the  coming 
months . 

*This  research  was  also  supported  by  the  Department  of  Energy  under 
Contract  DE-AS02-78ER04940  and  the  National  Science  Foundation  under  Grant 
NSF-CHE80-23747. 
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H.  VIBRATIONAL  EXCITATION  OF  CARBON  MONOXIDE  BY  COLLISIONS  WITH 

TRANSLATIONALLY  "HOT"  H  ATOMS  PRODUCED  VIA  LASER  PHOTO FRAGMENTATION* 

(C.  F.  Wood,  G,  W,  Flynn) 

(JSEP  work  unit  7,  1979-82) 

(Principal  Investigator:  G.  W,  Flynn  (212)  280-4162) 

Although  the  theory  of  translational  to  vibrational  energy  transfer 
processes  is  quite  well  developed  in  the  literature, ^  direct  experimental 
confirmation  of  these  theories  is  not  easily  obtainable.  The  reason  for 
this  paucity  of  experimental  information  has  been  largely  due  to  the  small 
excitation  cross  sections  for  T  •+•  V  processes  at  easily  obtainable  energies. 
Associated  with  this  problem  is  the  additional  difficulty  resolving  state -to- 
state  cross  sections  at  collisional  energies  large  enough  to  vibrationally 
excite  a  polyatomic  molecule;  at  translational  temperatures  high  enough  to 
cause  collisional  excitation,  the  higher  vibrational  levels  of  a  molecule 
are  also  likely  to  be  populated  causing  difficulties  in  Interpreting  data. 

Recently  a  bulb  technique  has  been  developed  that  overcomes  the 

problem  of  producing  translatlonally  energetic  atoms  in  the  presence  of  thermal- 

(2) 

ly  populated  vibrational  and  rotational  states  of  molecules.  An  excimer 
laser  is  used  to  photolyze  simple  molecules,  producing  translatlonally  fast 
photofragments.  The  subsequent  energetic  collisions  of  the  fast  atoms  with 
other  polyatomics  present  causes  IR  emission  from  the  polyatomic.  The  emis¬ 
sion  can  be  both  wavelength  and  time  resolved.  Since  most  of  the  theories 
of  T  -*■  V  processes  best  describe  an  atom  colliding  with  a  diatomic  molecule 
it  is  of  particular  interest  to  study  such  a  system.  We  have  begun  a  study 
of  the  vibrational  and  rotational  excitation  of  carbon  monoxide  caused  by 
collisions  with  fast  hydrogen  atoms  produced  in  the  photolysis  of  HBr. 
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H*  +  C0(v-0)  CO(v-n)  n-1,2,... 

CO(v“n)  -*■  C0(v=n-l)+hv  hv  ■  4.6p 


A  pulsed  Lumonics  TE-860S  laser  was  used  to  produce  the  193  nm  UV 
light  needed  to  dissociate  the  HBr.  The  laser  was  pulsed  at  low  repetition 
rates  (M.  Hz)  and  at  low  energies  (^  10  mJ/p).  Infrared  emission  from  the 
system  was  monitored  using  an  InSb  photovoltaic  infrared  detector  (77°K) . 

The  emission  was  wavelength  resolved  using  either  narrow  band  pass  filters 
O'*  120  cm  1  B.P.)  or  a  circular  variable  filter  (30  cm  1  B.P.).  The  signal 
from  the  detector  is  amplified,  digitized  with  a  Biomation  model  8100  transient 
recorder,  summed  in  a  signal  averager  (Nicolet  1170  or  Northern  575  A),  stored 
in  a  PDP  11/34  or  VAX  11/780  for  analysis.  A  flowing  stainless  steel  vacuum 
line  was  used  in  order  to  reduce  buildup  of  fragmentation  products  in  the 
cell. 

We  have  observed  T  V  energy  transfer  in  this  system  up  to  CO  (v  *  2) 
and  perhaps  higher.  The  emission  profiles  from  the  R  branch  (4.5  p)  and  from 
the  P  branch  (4.8  p)  are  quite  different.  (See  Figure  1.)  The  R  branch  profile 
is  dominated  by  a  fast  "spike"  that  decays  at  nearly  a  gas  kinetic  rate 
(6x10**  sec  1  torr  1),  followed  by  a  slower  rise  and  then  a  very  slow  decay. 

The  spike  may  be  caused  by  "hot"  rotational  levels  rapidly  relaxing  to  lower 
levels.  The  "spike"  is  unattenuated  by  passing  through  a  2  cm  long  cold  gas 
filter  containing  up  to  4  torr  CO,  (See  Figure  2.)  The  intensity  of  the  spike 
as  a  function  of  wavelength  was  studied  using  the  30  cm  1  band  pass  circular 
variable  filter.  This  study  indicates  that  the  H  atom  projectiles  are  ini¬ 
tially  exciting  the  CO  (v»l)  up  to  at  least  J-20.  The  spike  occurs  as  the 
high  J  states  rapidly  relax  to  what  is  probably  a  Boltzmann  rotational  dis¬ 
tribution.  The  final  rotational  temperature  has  not  been  determined  and,  due 
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FIGURE  2 


to  the  overlap  of  the  higher  vibrational  levels  with  the  R  branch  CO  (v  - 1) 
transitions,  cannot  easily  be  measured. 

The  Initial  fast  fall  in  the  P  branch  profile  decays  with  a  rate  cor¬ 
responding  to  the  process  CO(v-O)  +  C0(v-2)  -*•  CO(v- 1) +CO(v*  1).  The  second 
slower  rise  in  the  R  branch  rises  at  this  sane  rate.  Cold  gas  filter  studies 
were  perforaed  using  a  4.0  to  5.0  y  filter.  The  results  indicate  that  about 
702  of  the  excited  carbon  nonoxide  is  In  v»2  or  higher  levels.  (  See  Figure  3.  ) 

In  order  to  study  the  efficiency  of  energy  transfer  as  a  function  of  mass, 

13 

some  isotope  studies  have  been  performed.  CO  was  studied  under  the  same 

conditions  as  ^ CO .  Use  of  the  spectator  model ^  to  calculate  the  efficiency 

of  energy  transfer  shows  that  the  difference  in  energy  transfer  probability 

*  12  *  13 

would  be  on  the  order  of  12  greater  for  H  +  CO  than  H  +  CO.  Mo  significant 

differences  in  the  data  were  observed.  DBr  was  also  used  as  a  source  of  hot 

(3) 

atom  excitation.  The  same  energy  transfer  calculation  indicates  that  the 
D  atoms  can  transfer  up  to  302  of  their  energy  per  collision  as  opposed  to  the 
162  calculated  for  H  atoms  from  HBr.  Cold  gas  filter  studies  comparing  CO  excita¬ 
tion  with  both  D  atoms  and  H  atoms  do  not  show  appreciable  differences  in  the 
per  cent  excitation  in  the  v«2  or  higher  vibrational  levels.  However,  since 
it  is  impossible  to  distinguish  the  extent  of  excitation  into  v  ■  2  versus  higher 
levels  with  cold  gas  filter  studies,  the  possibility  remains  that  the  hot  D 
atoms  might  be  exciting  CO  into  the  v  >2  levels  more  efficiently  than  the  hot 
H  atoms. 
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FIGURE  3 
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I.  TRANSLATIONAL- TO-VIBRATIONAL  (T-V)  EXCITATION  OF  C02  BY  COLLISIONS 
WITH  "HOT"  HYDROGEN  ATOMS  PRODUCED  FROM  LASER  PHOTOLYSIS* 

(J.  0.  Chu,  G.  W.  Flynn) 

(JSEP  work  unit  7,  1979-82) 

(Principal  Investigator:  G,  W.  Flynn  (212)  28-4162) 

Energy  transfer  in  the  collislonal  excitation  or  relaxation  o.  c  e  internal 

states  of  a  molecule  is  a  topic  of  fundamental  importance  which  plays  a  signiii- 

cant  role  in  the  study  of  chemical  reactions. ^  Such  processes  involve  the 

simplest  type  of  energy  transfer  collision  event,  namely  the  energy  exchange 

between  the  translational  and  the  internal  vibrational  and  rotational  degrees 

of  freedom  of  a  molecule.  However,  in  contrast  to  the  numerous  investigations 

(2) 

of  the  collislonal  vibrational  relaxation  of  neutral  molecules,  there  have 
been  very  few  studies  of  the  reverse  process  of  collislonal  excitation.  This 
can  in  part  be  attributed  to  the  expected  small  excitation  cross  section  for 
such  energy  exchange  events  involving  particles  with  velocities  near  those  at  room 
temperature.  Relatively  difficult  molecular  beam  techniques  must  be  employed 
in  order  to  obtain  collislonal  energies  above  the  endothermic  threshold.  Fur¬ 
thermore,  beam  scattering  experiments  involving  neutral  collision  partners 
suffer  from  difficulties  associated  with  resolving  individual  vibration-rotation 
transitions.  Only  one  report  of  a  molecular  beam  scattering  experiment  has 

been  made  in  which  a  single  vibrationally  excited  state  has  been  resolved  under 

(3) 

single-collision  conditions, 

A  new  technique  has  recently  been  developed  for  studying  the  process 
of  translation  to  vibration  (T-V)  energy  transfer  and  translatlonally  enhanced 
reactive  events.  The  method  involves  excimer  laser  photolysis  of  bound  molecules 
to  produce  translatlonally  "hot"  (super thermal)  atoms  (or  radical  fragments)  in 
a  gas  sample,  followed  by  direct  time-resolved  detection  of  infrared  fluorescence 
from  vibrationally  excited  product  species.  This  technique  essentially  provides 
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a  sensitive  and  direct  probe  of  the  translation-to-vibration,  rotation  (T-V/R) 
energy  transfer  dynamics  that  occur  during  each  highly  energetic  collision 
event.  Furthermore,  infrared  luminescence  detection  of  T-V  excited  products 
can  readily  be  made  wavelength  state  specific,  and  subsequent  vibrational 
(and  possibly  rotational)  distributions  can  be  obtained  from  the  direct  wave¬ 
length  resolved  emission  profile.  In  cases  where  "hot",  highly  reactive  radical 
fragments  are  generated,  dynamical  competition  between  T-V  energy  transfer  and 
reactive  exchange  reactions  can  be  studied  in  detail. 

In  a  general  application  of  this  technique  to  the  study  of  "hot-atom" 
excitation,  initial  results  focussing  on  the  collisional  T-V  excitation  of  the 
molecules  HC1,  DC1,  HBr,  DBr,  CO^,  and  N^O  by  hot  H  and  D  atoms  produced  by 
ArF  laser  photolysis  of  H61/DC1,  HBr/DBr  or  HjS  have  been  reported. ^ ^ 

The  present  study  extends  this  work,  in  particular,  to  wavelength  resolution 
of  the  4-5  pm  IR  emission  observed  from  excited  by  translation-to-vibration 
(T-V)  energy  transfer  collisions  with  hot  H  atoms  produced  by  the  ArF  laser 
photolysis  of  HBr. 

The  output  of  a  pulsed  Lumonics  TE-860S  excimer  laser  operated  at  193  run 
(ArF,  pulse  width  10  ns  FWHM)  was  used  to  irradiate  a  cyclindrical  Pyrex  sample 
cell,  which  was  fitted  with  Suprasil  windows.  A  conventional  stainless  steel 
vacuum  flow  system  was  used  to  flow,  from  a  pre-mixed  sample  in  a  glass  bulb, 
fresh  gas  samples  through  the  fluorescence  cell  to  prevent  gas  depletion  and 
the  build  up  of  reaction  products.  Typically  all  experiments  were  performed 
at  low  fluences  in  order  to  avoid  multiple  photon  events.  Infrared  emission 
was  detected  through  a  NaCl  window  at  right  angles  to  the  laser  beam  with  a 
Judson  Infrared  InSb  detector  (1  ys  risetime,  77°K) .  Either  broadband  inter¬ 
ference  filters  or  a  narrowband  OCLI  circular  variable  filter  (%  30  cm  \  FWHM) 
were  used  to  view  only  the  desired  wavelength  regions.  When  working  with  the 
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circular  variable  filter  (CVF) ,  a  KC1  or  CaF^  lens  was  employed  to  focus  the 
IR  emission  through  the  narrow  entrance  slit  (.125")  of  the  CVF  housing  and  onto 
the  detector.  FTIR  spectral  analysis  was  used  to  provide  an  absolute  calibra¬ 
tion  of  the  wavelength  and  resolution  function  of  the  circular  variable  filter. 


CC>2  hot  band  emission  was  detected  by  cold  gas  filtering  techniques  whereby 
a  gas  filter  cell  (3.5  cm  optical  length,  filled  with  CO^)  was  placed  between 
the  detector  and  the  cell  window  to  remove  emission  from  the  00 °1  -  00°0  and 
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transitions. 


Wavelength  resolved  emission  from  vibrationally  excited  CO^ 


molecules 


produced  via  collisional  (T-V)  energy  exchange  with  "hot"  hydrogen  atoms  have 


been  observed  for  the  first  time.  The  relative  rovibrational  populations  of 


the  excited  CO^  product  were  estimated  from  the  amplitudes  of  the  initial  fast 
rise  component  of  fluorescence  (see  Fig.  la).  Figure  2  shows  the  CVF  (band¬ 
width,  %30  cm  *  FWHM)  spectrum  of  the  amplitude  of  the  fast  rise  emission 
associated  with  excited  CO^  product  superimposed  over  a  FTIR  (IBM  85,  0.5  cm  * 
resolution)  absorption  spectrum  of  room  temperature  CO^.  The  amplitudes  for 
the  fast  rise  emission  process  were  determined  by  fitting  decay  times  and 
amplitudes  to  latter  portions  of  the  fluorescence  signals  and  then  extrapola¬ 


ting  these  amplitudes  to  zero  time  to  give  an  amplitude  for  the  initial  fast 


rise  component.  The  CVF  emission  profile  does  show  coarse  rotational  structures 
which  are  noticeably  red  shifted  from  the  (00°1  -00°0)  fundamental  transi¬ 
tion  indicating  considerable  hot  band  emission. 

Cold  gas  filtering  (CGF)  techniques  were  used  as  a  probe  to  study  and 


identify  the  state  and  nature  of  the  highly  excited  CC^  product.  A  typical 
fluorescence  observed  through  an  empty  filter  cell  exhibits  three  characteristic 


exponential  components  in  its  emission  profile  as  shown  in  Fig.  la.  To  isolate 
the  emission  from  highly  excited  CO^  states,  the  cold  gas  filter  is  used  to 
effectively  remove  all  emission  terminating  on  the  00°0  and  01*0  levels 
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j-V  EXCITED  C02_CVF  WAVELENGTH  RESOLVED  SPECTRUM 
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Figure 


filter  pressure  7  torr,  length  3.5  cm).  The  results  obtained  using 
gas  filter  (Fig.  lb)  show  a  small  attenuation  of  the  initial  fast 
onent  whereas  the  slowly  rising  component  is  completely  absorbed  by 
r  cell.  This  implies,  in  agreement  with  the  circular  variable  filter 
data,  that  the  fast  component  has  contributions  largely  from  hot  band 
ns  and  possibly,  from  highly  excited  rotational  states.  Subsequent 
and  R-branch  filter  results  indicate  that  this  emission  is  not  due  to 
cited  J  states  but  rather  to  vibrational  states  higher  in  energy  than 
e  decay  rate  of  the  4.3  pm  CO^  fluorescence  observed  through  the  filled 
filter  cell  (see  Fig.  lb)  was  also  measured.  Figure  3  shows  a  Stern- 
lot  of  the  decay  rate  versus  pressure  which  yields  a  relaxation  rate 
12  msec  ^  per  torr  of  CC^.  Similarly  at  a  constant  pressure, 
ate  of  approximately  20  msec  ^  per  torr  of  HBr  was  observed.  Without 
liter  cell,  rates  of  127  +  16  msec  1  per  torr  of  CO,,  and  144  +  11  msec 
IBr  were  also  obtained  for  the  slowly  rising  portion  of  the  fluorescence 

s  large  red  shift  in  the  wavelength  resolved  CVF  spectrum  of  the  col- 
y  excited  CC^  molecule  is  strong  evidence  that  the  majority  of  the 
fluorescence  emission  is  not  due  directly  to  the  fundamental  vibra- 
ansition  of  the  antisymmetric  stretch  00°1  -  00°0  of  CO^  centered  at 
Rather  the  emission  corresponds  to  hot  band  transitions  nmP£  -  nmP (2.-1) 
In  particular,  this  suggests  that  during  collisional  excitation,  the 
:rgy  transfer  which  produces  CO^  molecules  having  at  least  one  quantum 
iferentially  excites  combination  or  overtone  levels  in  the  target 
,  An  estimate  of  the  collisional  excitation  to  the  mode  was  pre- 
lade  and  considered  for  the  case  with  H-atoms.  Bass^  has  derived 
expressions  for  the  collisional  vibrational  excitation  of  the 
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mode  of  CO which  indicate  that  from  the  available  center  of  mass  translational 
energy  the  maximum  vibrational  energy  available  for  the  molecule  is  4290  cm  ^ . 
This  suggests  that  single  collisions  cannot  excite  levels  higher  than  v*  1  in 
the  manifold  (i.e.  nm^l),  and  the  minimum  translational  energy  necessary  for 
pure  excitation  (i.e.  00°1)  is  11186  cm  Moreover,  the  combination  states 
or  hot  band  transitions  of  CC^  which  have  energies  less  than  4290  cm  ^  can  be 
collisionally  excited. 

Vibrational  energy  transfer  in  CO2  molecules  has  been  studied  extensively 
and  a  number  of  resonant  vibration-vibration  energy  transfer  ates  have  either 
been  well  established  or  recently  measured.  Weitz  et.  al.  ^  have  studied  the 
equilibration  of  the  (10°0)  and  (02°0)  lower  laser  states  in  CO2  and  concluded 
that  the  coupling  rate  between  these  states  is  very  rapid,  _>  1100  msec  ^torr  ^ , 

In  addition  .these  Fermi  mixed  levels  are  also  coupled  to  the  bending  modes  of 


C02(02°0/0220/10°0)  +  C02(00°0)  ~ — ^  C02(0l10)  +  C02(0110) 

at  a  slower  forward  rate  of  389  +  75  msec  ^torr  ^ .  By  comparison  with  the 
present  cold  gas  filter  results,  this  (V-V)  transfer  rate  is  also  observed  as 
a  relaxation  process  following  T-V  excitation.  This  suggests  that  the  states 
emitting  through  the  cold  gas  filter  are  tightly  coupled  by  collisions  to  the 
bending  manifold.  Moreover,  in  consideration  of  the  estimated  maximum  4290  cm 
excitation  at  the  available  collisional  energies,  the  possible  excited  C02 
"hot  band"  transitions  are  (10°1  -  10°0),  (02°1  -  02°0)  and  (02^1  -  0220) . 

Under  the  present  conditions,  we  are  unable  to  resolve  emission  from  these 
closely  coupled  Fermi  mixed  or  degenerate  combination  states  of  CC^. 

The  studies  we  plan  to  pursue  to  monitor  the  product  vibrational  state 
populations  of  CO2  involves  time  resolved  absorption  spectroscopy  of  the  target 
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molecules.  This  method  employs  a  tunable  narrowband  IR  diode  laser  to  monitor 
and  discriminate  the  absorption  of  a  specific  v,  J  -*■  v+1,  J  +  1  transition 
during  and  immediately  following  the  T-V  excitation  process.  In  addition,  time 
dependent  absorption  changes  resulting  from  collisional  process  will  provide 
information  about  some  of  the  energy  rearrangement  and  relaxation  events  which 
cannot  be  determined  with  interference  filters.  We  have  also  obtained  consider¬ 
able  data  in  the  study  of  energy  variation  effects  (i.e.  mass  and  velocity  de¬ 
pendence  of  projectiles)  on  the  (T-V)  energy  transfer  dynamics  that  occur 
during  the  collisional  excitation  events.  These  studies  involve  laser  photo¬ 
lysis  at  193  nm  (ArF)  and  248  nm  (KrF)  of  DBr,  ^S,  HC1  and  Cl^  molecules  as 
sources  of  light  and  heavy  "hot"  atoms  to  determine  the  effect  of  mass  and 
electronics  structure  on  the  excitation  of  CO^.  Further  energy  variation  is 
also  available  by  using  a  doubled  excimer  pumped  dye  laser  pulse  in  the  region 
of  250-300  nm  for  photolysis. 

*This  research  was  also  supported  by  the  Department  of  Energy  under  Contract 
DE-AS02-78ERQ4940  and  the  National  Science  Foundation  under  Grant  NSF-CHE80-23747. 
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IV.  PICOSECOND  ENERGY  TRANSFER  AND  PHOTOFRAGMENTATION  SPECTROSCOPY 

A.  GENERATION  OF  SINGLET  OXYGEN,  102,  FROM  LASER  INDUCED  PHOTOREACTIONS 

(S.  Y.  Hou,  C.  G.  Dupuy,  M.  J.  McAuliffe,  K.  B.  Eisenthal) 

(JSEP  work  units  8  and  9,  1979-82) 

(Principal  Investigator:  K.  B.  Eisenthal  (212)  280-3175) 

Singlet  oxygen  is  significantly  involved  in  chemical,  physical,  biological, 

atmospheric,  and  in  recent  years  medical  phenomena.  Its  spectroscopy,  its 

novel  chemistry,  its  role  in  energy  and  materials  degradation,  have  made 

this  field  an  extremely  important  and  fertile  one. ^  We  have  been  engaged 

in  using  picosecond  lasers  to  study  the  photodissociation  pathways  of  molecules 

which  yield  ^C>2  as  one  of  the  photofragments.  These  processes  are  an  important 

(2) 

class  of  excited  state  adiabatic  photo react ions  as  well  as  providing  a  con¬ 

venient  photochemical  source  for  producing  *02  and  studying  its  photophysics 
and  chemistry. 

(3) 

We  have  recently  discovered  that  one  of  the  paths  in  the  photodissocia¬ 
tion  of  an  anthracene  endoperoxide  produces  not  only  102  but  also  an  excited 
aromatic  fragment;  this  represents  one  of  the  few  cases  of  a  photodissociation 
producing  two  electronically  excited  fragments. 


The  importance  of  this  pathway,  relative  to  other  decay  routes,  has  a  direct 
bearing  on,  and  can  be  used  to  test  correlation  diagram  predictions  of  the 
photodissociation  of  the  endoperoxides  as  well  as  the  reaction  of  *C>2  with  the 
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An  important  aspect  of  this  photodissociative  route  is  the  particular 
state  in  which  the  oxygen  is  formed.  By  considering  the  energetics  of  the 
reaction  we  have  been  able  to  determine  the  particular  excited  state  in  which 
^C>2  is  produced.  From  Fig.  1  we  can  see  that  it  is  only  the  *A^  state  which 
is  populated  in  this  photoreaction,*  the  ^1*  state  being  too  high  in  energy  to 
be  reached.  This  finding  represents  the  first  time  that  the  state  of  ^0^ 
resulting  from  a  photodissociation  has  been  established. 


Figure  1.  Energetics  of  endoperoxlde  photodissociation, 
X  “  266  nm. 


*This  work  is  also  supported  by  the  Air  Force  Office  of  Scientific  Research  under 

Grant  AFOSR-81-0009-A  and  the  National  Science  Foundation  under  Grant  NSF-CHE81- 

21945. 

(1)  For  a  theoretical  and  experimental  discussion  of  ^2 ,  its  spectroscopy, 
and  reactions,  see:  Vasserman,  H.H.;  Murray,  R.  W.,  Eds.  "Singlet 
Oxygen",  Academic  Press:  New  York,  1979. 

(2)  Fflrster,  Th.  Pure  Appl.  Chem.  34,  225  (1973) . 

(3)  S.  Y.  Hou,  C.  G.  Dupuy,  M.  J.  McAuliffe,  D.  A.  Hrovat,  and  K.  B.  Eisenthal, 
"Picosecond  Laser  Study  of  the  Adiabatic  Photodissociation  of  an  Endo- 
peroxide".  Jounral  of  the  American  Chemical  Society,  103,  6982  (1981). 


B.  INTRAMOLECULAR  CHARGE  REDISTRIBUTION  AND  ANOMALOUS  FLUORESCENCE 

(Ying  Wang,  K.  B.  Eisenthal) 

(JSEP  work  units  8  and  9,  1979-82) 

(Principal  Investigator:  K.  B.  Eisenthal  (212)  280-3175) 

Photoexcitation  can  lead  to  significant  charge  separation  in  a  molecule 
thereby  producing  a  large  molecular  dipole/^  ^  The  phenomenon  of  photo- 
induced  charge  separation  is  one  of  the  most  important  primary  processes  in 
photochemistry  and  photobiology.  For  the  molecule  p-(dimethylamino)  ben- 
zonitrile,  DMABN, 


DMABN 


this  charge  redistribution  process  leads  to  the  appearance  of  an  anomalous 
dual  fluorescence  in  polar  media. ^  ^  A  bond  twisting  event  between  the 
amino  group  and  the  benzene  ring  and  solvation  of  the  excited  molecule  are 

involved  in  the  charge  rearrangement. 

(8) 

In  our  picosecond  laser  studies v  we  have  addressed  the  key  issues  of 
(1)  the  relationship  between  the  two  fluorescences,  one  in  the  ultraviolet  at 
350  nm  and  the  other  in  the  visible  at  480  nm,  (2)  the  nature  of  the  short 
wavelength  emission  at  350  nm  and  (3)  the  role  of  the  solvent  on  the  dynamics 
of  the  charge  transfer  process. 

Our  results  definitively  show  that  an  equilibrium  is  rapidly  established 
between  the  short  wavelength  emission  from  S*  and  the  visible  emission  from 
CT*,  (Figs,  1  and  2),  and  can  be  depicted  as  follows 

S*  Z~,  -  CT* 

/  \ 
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and  the  points  are  experimental  data. 


The  emitting  state  S*  we  conclude  is  due  to  the  excited  singlet  state 
which  borrows  intensity  from  by  vibronic  coupling.  The  measured  time 
of  20  ps  obtained  from  the  reaction  rates  is  interpreted  as  the  time  required 
for  the  molecule  to  relax  conformatlonally  to  the  twisted  charge  transfer 
geometry  as  shown  below. 

*This  work  is  also  supported  by  the  Air  Force  Office  of  Scientific  Research 
under  Grant  AFOSR— 81— 0009— A  and  the  National  Science  Foundation  under  Grant 
NSF-CHE81-21945 . 


(1)  E.  Lippert,  W.  LUder,  and  H.  Boos,  Advan.  Mol.  Spectrosc.  443  (1962). 

(2)  0.  S.  Khalil,  J.  L.  Meeks,  and  S.  P.  McGlynn,  Chem.  Phys,  Lett.  39, 
457  (1976). 

(3)  H.  Dodiuk  and  E.  M.  Kosower,  Chem.  Phys.  Lett.  34,  253  (1975). 

(4)  E.  A.  Chandross  in  "The  Exclplex",  M.  Gordon  and  W.  R.  Ware,  Ed., 
Academic  Press,  New  York,  1975,  p.  187. 

(5)  K.  Rotklewicz,  K.  H.  Grellmann,  and  2.  R.  Grabowski,  Chem.  Phys.  Lett. 
19,  315  (1973), 

(6)  E.  Kir ko-Kaminska ,  K.  Rotklewicz,  and  A.  Grabowska,  Chem.  Phys.  Lett. 
58.  379  (1978). 

(7)  Z.  R.  Grabowski,  K.  Rotklewicz,  A,  Siemiarczuk,  D.  J.  Cowley,  and 
W.  Baumann,  Nouv.  J.  Chim.  _3 ,  443  (1979). 

(8)  Ying  Wang,  M.  McAuliffe,  F.  Novak,  and  K.  B.  Eisenthal,  "Picosecond 
Dynamics  of  Twisted  Internal  Charge-Transfer  Phenomena",  The  Journal 
of  Physical  Cehmistry,  85,  3736  (1981). 


125 


C.  SMALL  CHAIN  RELAXATION  OBSERVED  BY  INTRAMOLECULAR  EXCITED  STATE 

CHARGE  TRANSFER 

(Ying  Wang,  M.  Crawford,  K.  B.  Eisenthal) 

(JSEP  work  units  8  and  9,  1979-82) 

(Principal  Investigation:  K.  B,  Eisenthal  (212)  280-3175) 

The  nature  of  rapid  internal  rotational  motions  can  be  obtained  by 
determining  the  dynamics  of  reaction  of  the  end  groups  on  the  opposite  sides 
of  the  chain.  ^  In  our  earlier  studies^*  ^  we  determined  that 

exciplex  (excited  charge  transfer  complex)  formation  in  the  small  chain 
system. 


requires  a  favorable  geometry  in  nonpolar  solvents,  most  likely  a  sandwich 
type  structure. 

We  have  used  picosecond  laser  excitation  and  streak  camera-optical  multi¬ 
channel  detection  of  exciplex  fluorescence  at  different  emission  wavelengths 
and  as  a  function  of  solvent  viscosity  to  characterize  the  dynamics  of  the 
chain  motions. 

For  the  model  compound  we  have  studied  we  conclude  that: 

(1)  The  dynamics  of  exciplex  formation  are  influenced  by  the  starting  ground 
state  conformations.  For  the  model  compound  studied,  it  was  found  that  ex- 
ciplexes  with  different  emission  energies  were  formed  independently  from  two 
distinct  groups  of  ground  state  conformers  at  different  rates,  Table  1, 

The  structures  of  these  conformers  are  not  known  but  possibly  differ  in  the 
relative  orientation  of  the  anthracene  and  the  dimethylaniline  chromophores . 

(2)  The  observed  viscosity  dependence  we  find,  Fig,  1,  combined  with  litera¬ 
ture  values  for  other  three  chain  systems,  points  to  the  methylene 

chain  motions  as  the  key  factor  in  determining  the  viscosity  dependence  of 
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Figure  1.  Plot  of  log  of  rotational  relaxation  time  vs.  log 

of  solvent  viscosity  for  two  conformers  of  anthracene- 
(CH_),-DMA  in  isopentane,  hexane,  decane  and  detradecane 


The  attached  end  groups  alter  the  formation  rates  but  not  the 


dependence  of  the  rates  for  most  of  the  cases  considered.  We  observe 
diverse  cases  that  the  lifetime  of  chain  relaxation  scales  as 


experimentally  observed  exponential  rise  for  intramolecular  exciplex 
suggest  that  a  single  chain  relaxation  mode,  perhaps  the  large  scale 
d  relaxation  mode,  is  dominant,  at  least  in  the  picosecond  to  nano- 
me  regime  we  have  examined. 
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V.  GENERATION  AND  CONTROL  OF  RADIATION 

A.  OPTICAL  COHERENT  TRANSIENT  SPECTROSCOPY 

(R.  Beach,  E.  Whittaker,  E.  Xu,  R.  Kichinski,  M.  Glick,  S.  Hartmann) 

(JSEP  work  units  10  and  11,  1979-82) 

(Principal  Investigator:  S.  R.  Hartmann  (212)  280-3272) 

During  the  last  year,  our  laboratory  has  undergone  a  major  overhaul. 

We  acquired  a  second  Quanta  Ray  YAG  laser  which  motivated  us  to  develop  a 
centralized  laser  facility  in  which  a  pair  of  YAG  lasers  could  be  used  simul¬ 
taneously  in  several  different  laboratories.  Stable  beam  steering  optics  had 
to  be  developed  to  separately  image  each  YAG  into  the  several  laboratories. 
Currently,  the  two  YAGs  are  located  next  to  each  other  with  ^30%  of  their 
output  split  off  for  use  in  the  room  in  which  they  are  housed  and  the  other 
70%  steered  into  a  duct  suspended  just  below  the  ceiling  so  as  to  be  out  of 
the  way  of  traffic  within  the  rooms.  The  duct  work  containing  the  laser 
beams  passes  into  two  other  laboratories  each  of  which  uses  half  of  the  remaining 
70%  of  the  laser  light.  We  are  now  in  the  process  of  extending  the  laser  beam 
line  into  a  fourth  laboratory. 

The  YAG  lasers  pump  tunable  dye  lasers  in  the  different  labs.  Here  again 
we  revamped  our  system  by  building  a  set  of  ten  new  dye  lasers  and  amplifiers 
of  the  Littman  design  as  Improved  by  Frank  Tompkins  of  Argonne  Laboratories. 

These  new  lasers  are  easy  to  use,  stable,  have  a  narrow  linewidth  of  the  order 
of  1  GHz,  and  are  a  significant  improvement  over  our  previous  Hansch  design. 
There  is  also  an  advantage  in  several  workers  using  identical  lasers  as  exper¬ 
tise,  new  approaches,  good  ideas,  etc.  are  readily  transferred  and  adapted. 

The  feasibility  of  the  system  was  proven  when  we  saw  stable  photon  echoes  in 
the  laboratory  furthest  removed  from  the  YAG  lasers.  In  fact,  no  retuning 
or  refocusing  were  necessary  to  generate  echoes  wh®  the  system  was  turned  back 
on  a  week  later.  This  observation  was  quite  conforting  as  it  tested  the  stability 
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not  only  of  the  YAG  beam  transport  system  but  of  the  dye  lasers  and  the  echo 
generating  set-up  as  well.  Since  the  echoes  were  generated  by  using  two  dif¬ 
ferent  dye  lasers  to  provide  the  excitation  pulse  pairs,  rather  than  an  optical 
delay  line  with  a  single  dye  laser,  the  test  was  quite  demanding.  The  very  narrow 
pulsed  dye  laser  linewidth  is  achieved  by  pumping  the  dye  lasers  lightly,  so  laser 
amplifiers  are  necessary  to  attain  sufficient  laser  pulse  power.  These  also 
have  been  totally  successful,  as  demonstrated  by  the  above-mentioned  echo  experi¬ 
ment.  We  now  routinely  and  simultaneously  achieve  pairs  of  100  kW  1  GHz-wide 
laser  pulses  in  each  of  three  laboratories.  Before  our  change  over,  all  pumping 
was  done  by  separate  homemade  Nj  lasers,  each  of  which  put  out  no  more  than  100 
to  200  kW  of  pump  power  when  and  if  they  worked.  Our  dye  lasers  now  provide  as 
much  power  as  our  pump  lasers  previously  provided  and  the  reliability  and 
stability  have  been  significantly  increased. 

The  great  advantage  of  having  two  YAG  lasers  is  that  we  can  create  and 
control  a  continuously  variable  time  separation  between  them  -  with  remote  elec¬ 
tronics  -and  thereby  do  experiments  to  directly  study  echo  modulation  and  relaxa¬ 
tion.  Implementation  of  the  electronic  system  to  achieve  this  control  was  itself 
a  major  accomplishment.  Finally,  the  problem  of  a  residual  one  or  two  nanosecond 
jitter  in  pulse  separation  is  circumvented  by  a  detection  scheme  which  directly 
measures  t  pulse  separation  to  within  a  few  hundred  picoseconds  and  labels  the 
corresponding  data  accordingly. 

At  the  tail  end  of  the  laboratory  overhaul  we  are  replacing  our  obsolete 
PDP8  computer  system  with  a  multi-user  LSI  11/23.  Separate  stations  will  be 
placed  in  the  individual  laboratories  to  allow  for  a  more  efficient  operation 
than  was  possible  when  our  one  PDPS  had  to  be  physically  wheeled  from  lab  to  lab 
as  demand  required.  ^ 

During  the  laboratory  overhaul  research  progress  proceeded  as  manifested 
in  the  Billiard  Ball  Echo  model. ^  Photon  echoes  in  gaseous  media  have  up 
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to  now  been  treated  in  a  rather  formal  manner.  This  is  unfortunate  particularly 
because  photon  echo  phenomena  display  a  wide  variety  of  forms.  Aside  from  the 
two  pulse  photon  echo  and  the  three  pulse  stimulated  photon  echo  there  are  Raman 
echoes,  tri- level  echoes,  grating  echoes,  two  photon  echoes,  etc.  Each  echo 
has  its  particular  application  and  a  simple  uniform  treatment  would  facilitate 
the  proper  general  understanding  and  application  of  the  echo  technique.  Such 
a  treatment  is  now  available  with  the  Billiard  Ball  Echo  model 

The  Billiard  Ball  Echo  model  takes  an  unconvential  view  of  nonlinear 
optics.  The  key  is  atomic  recoil.  When  an  atom  is  excited  to  a  higher  atomic 
state  by  absorbing  light  it  must  recoil  with  the  momentum  of  the  absorbed  photon. 
An  atom  which  has  been  subjected  to  irradiation  by  a  light  pulse  thus  finds 
itself  in  a  linear  superposition  of  its  ground  and  excited  states  with  its 
excited  state  recoiling  away  from  its  ground  state.  An/2  light  pulse  has  a 
502  probability  of  putting  an  atom  into  its  excited  state  so  the  atom  has  equal 
probabilities  of  being  left  alone  and  not  recoiling  or  being  excited  and  re¬ 
coiling  with  the  momentum  of  the  absorbed  photon.  Since  the  coherently  excited 
atoms  are  in  a  linear  superposition  of  these  states  they  effectively  split 
apart.  But  the  atom's  ability  to  radiate  coherently  proportional  to  its  radia- 
tive  dipole  moment  <pe  r>  (where  p  is  the  dipole  moment  operator  and  k  is 
the  photon  wave  vector)  which  depends  on  the  overlap  of  these  ground  and  ex¬ 
cited  states.  Thus  as  the  "atom"  splits  apart  the  coherent  radiation  it  emits 
diminishes. 

The  universal  trick  to  generate  photon  echoes  of  whatever  type  is 
to  irradiate  the  atoms  with  a  such  a  series  of  pulses  as  to  split  its  component 
parts  apart  and  then  to  bring  them  back  together  again.  The  analysis  is  just 
a  simple  matter  of  mechanics.  Each  laser  pulse  resonant  with  an  optical 
transition  splits  the  atom  apart,  its  component  parts  recoiling  as  the  laser 
pulse  stimulates  the  atom  to  emit  or  absorb  a  photon.  We  need  only  keep  track 
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of  atomic  recoil  to  determine  when  and  if  an  echo  will  occur.  Since  all  atoms 
recoil  in  identical  manner  we  need  only  consider  one  atom.  We  illustrate  this 
idea  with  the  well  known,  basic,  two-pulse  case.  In  Figure  1  we  display  the 
atomic  recoil  for  the  simple  case  where  the  two  excitation  pulses  are  collinear. 
The  ground  state  trajectory  is  represented  by  a  solid  line  and  the  excited 
state  trajectory  by  a  dashed  line.  The  atom  is  initially  in  its  ground  state 
and  since  it  is  therefore  not  recoiling  its  trajectory  is  horizontal  as  in¬ 
dicated.  At  t  *  0  a  n/2  laser  pulse  excites  the  atom  into  a  linear  superposi¬ 
tion  of  its  ground  and  excited  states  so  the  trajectory  behaves  as  shown  in 
Fig.  1.  Only  the  excited  state  trajectory  indicates  recoil  as  only  the  excited 
state  has  absorbed  a  photon.  The  trajectory  slope  is  given  as -fik/m  where  m 
is  the  atomic  mass  to  represent  the  recoil  velocity.  The  short  wiggly  line 
at  the  vertex  at  t  *  0  indicates  the  absorption  of  a  photon.  The  atom  proceeds 
to  separate  at  a  uniform  rate  until  t  =  t  when  a  second  collinear  excitation 
pulse  arrives.  This  pulse  is  adjusted  to  have  area  it  so  its  effect  is  as 
shown.  The  ground  state  atom  component  has  100%  probability  of  absorbing  a 
tr  pulse  (the  definition  of  a  n  pulse)  and  so  it  is  transformed  entirely  into 
an  excited  state  component  with  a  recoil  trajectory  of  slope-fik/m.  The  excited 
state  component,  on  the  other  hand,  is  stimulated  with  100%  probability  to  emit 
a  photon  of  momentum 'hk.  It  thus  recoils  backward  and  since  it  was  originally 
created  by  the  absorption  of  an  identical  photon,  it  stops  in  its  tracks. 

We  thus  have  the  situation  of  the  newly  created  ground  state  component  standing 
motionless  while  the  newly  created  excited  state  component  approaches  it  with 
velocity "hk/m.  When  overlap  reoccurs  the  radiative  moment  reforms  and  a  burst 
of  radiation  is  emitted.  This  is  the  photon  echo.  The  short  wiggly  lines  in 
the  diagram  indicate  the  absorption  or  emission  of  radiaton  due  to  the  action 
of  the  laser  excitation  pulses,  or  emission  from  the  renewed  overlap  of  the 
component  parts  of  the  atoms. 
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Figure  1.  Recoil  diagram  for  a  two  pulse  photon  echo  formed  by  applying 
a  ir/2  pulse  at  t  *  0  followed  by  a  ir  pulse  at  t »  t.  The 
echo  forms  at  t*2x.  The  ground  and  excited  state  trajectories 
are  indicated  by  solid  and  dashed  lines  respectively.  The 
short  wiggly  lines  at  the  vertices  represent  photons  that 
are  either  absorbed  or  emitted. 
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The  above  analysis  requires  no  formulas  or  complicated  mathematical 
analysis  -  no  SchrHdinger ' s  equation!  Nevertheless,  the  analysis  is  readily 
extended  to  much  more  complicated  scenarios  where  noncollinear  laser  pulses 
separately  resonant  to  different  optical  transitions  are  applied. 

This  analysis  is  rather  unconventional  and  experts  in  the  field  of  non¬ 
linear  optics  have  been  uncomfortable  with  it  at  first  since  it  doesn't  im¬ 
mediately  seem  to  oorollate  with  well-known  phenomena.  Consider  for  example 
the  free  radiation  decay  from  a  gas  of  atoms  excited  by  a  single  resonant  laser 
pulse.  In  its  simplest  form  a  tt/2  pulse  coherently  excites  the  atoms  in  a 
gaseous  sample  into  a  linear  superposition  of  their  ground  and  excited  atomic 
states  and  the  atoms  then  proceed  to  radiate  coherently.  The  convential  view 
is  that  since  the  atoms  of  a  gas  are  moving,  they  radiate  at  Doppler  shifted 
frequencies.  The  atoms  are  in  a  Maxwell  Boltzmann  distribution  for  the  given  gas 
temperature  T,  so  the  coherent  radiation  disappears  due  to  Doppler  dephasing 
in  a  time  depending  on  T. 


The  dephaslng  time  is  the  time  it  takes  atoms  to  move  away  from  each 
other  by  some  critical  distance.  When  atoms  initially  excited  in  phase  separate 
by  the  order  of  an  optical  wavelength  they  destructively  interfere.  Writing 


vDoppler  C°  rePresent  the  Doppler  velocity  spread,  the  dephaslng  time  T<jep|usjI1g 
is  determined  from 


vDopplerTdephasing 


optical 


"here  *opticai  is  the  optical  wavelength  of  the  radiative  transition  excited. 

This  can  be  cast  in  a  more  obvious  form  by  writing  v  ,  in  terms  of  a 

Doppler 

d.  Broglie  «wle»gth  *charm4^  of  the  .to.,  .vD,pH<ir  -  'SAth.nn,r  The 


dephasing  time  is  then  expressed  as 
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Tdephasing  \jptical\hermal 


The  thermal  wavelength  ^thermai  0311  be  seen  33  the  size  of  a  wave  packet 
constructed  from  the  Maxwellian  distribution  of  momentum  states  which  describe 


the  gas  at  temperature  T.  In  the  context  of  the  Billiard  Ball  model  it  is 
what  we  might  consider  the  size  of  the  atom.  Thus  the  time  t 
for  the  atom  to  split  apart  is  given  by 


splitting 


vrecoil  ^splitting  ^thermal 

where  vreco^  is  the  atomic  recoil  velocity.  This  recoil  velocity  is  obtained 
by  equating  atomic  recoil  momentum  to  photon  momentum. 


so 


mv 

recoil 

mv  . , 
recoil 


=  -fik 

“  -KVA 

optical 


and  t  ,  .  _ „  .  becomes 
splitting 


Splitting  \jptical*thermal 


which  is  identical  to  the  expression  for  t 


dephasing’ 


T  m  T 

splitting  dephaslng 

and  the  unconventional  point  of  view  taken  by  the  Billiard  Ball  model  leads 
to  conclusions  in  accord  with  previous  experience.  A  thorough  discussion  of 
the  Billiard  Ball  Echo  model  appears  in  Ref.  1. 
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Until  now  photon  echoes  have  been  confined  to  those  cases  where  the  laser 
excitation  is  within  a  wave  numbers  of  exact  resonance.  We  now  report  pre¬ 
liminary  experiments  which  lift  this  restriction.  Using  a  laser  working  at 
6854  R  we  have  succeeded  in  generating  second  harmonic  radiation  (SHG)  by 
applying  a  transverse  magnetic  field.  The  second  narmonic  of  6854  R  is  at 
3427  R  the  frequency  of  the  3D-3S  transition  in  Na;  the  intermediate  3P  state 

at  5890  R  is  far  removed  from  the  6854  R  laser  frequency.  At  a  number  density 
16  —3 

of  order  10  cm  and  in  a  magnetic  field  of  ^100  gauss  we  obtain  a  harmonic 
generation  efficiency  of  2  x  10  ^  Our  analysis  suggests  that  much  higher 
efficiencies  are  possible  at  higher  values  of  the  transverse  magnetic  field. 

The  detection  of  SHG  means  that  the  atoms  have  been  excited  in  a  linear 
superposition  of  the  3S  and  the  3D  states  by  radiation  of  half  the  transition 
frequency.  We  plan  next  to  generate  a  two-photon  echo  directly  on  the  3S-3D 
transition  using  laser  excitation  at  half  the  transition  frequency  and  detect¬ 
ing  the  echo  at  the  full  transition  frequency. 

This  experiment  has  not  been  done  but  we  have  tried  another  experiment 
which  suggests  that  the  direct  two-photon  echo  is  feasible.  The  experiment 
is  to  irradiate  Na  atoms  with  two  laser  pulses,  one  at  5890  R  and  one  at  6854  R. 
We  observe  an  echo  at  8200  R.  This  is  an  entirely  new  kind  of  two  pulse  photon 
echo  experiment.  It  is  readily  analysed  by  means  of  the  Billiard  Ball  Echo 
model.  Let  the  5890  R  pulse  be  applied  at  t  =  0.  This  excites  the  3P  state 
which  recoils  as  indicated  in  Fig.  2.  The  6854  R  pulse  then  excites  the  3D  state 
at  t~T,  each  atom  absorbing  two  6854  R  photons  and  recoiling  as  indicated  in 
Fig.  2.  Following  the  recoil  trajectories  in  Fig.  2  we  note  that  component 
parts  of  the  atoms  corresponding  to  the  3P  and  the  3D  states  should  cross  at 


T 

e 


5890  R 

5890  R  -  i  6854  R 


2.39  t  . 
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Recoil  diagram  for  a  two-pulse  -  two-photon  trilevel  echo. 

A  tt/2  pulse  exciting  a  single  photon  transitions  is  followed 
by  a  ir/2  pulse  exciting  a  two-photon  transition.  The  echo  is 
emitted  when  the  recoiling  excited  states  overlap. 
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In  our  experiment  when  t  -  20  nsec  we  see  an  echo  at  =  47.8  nsec.  This  experi¬ 
ment  demonstrates  the  viability  of  the  basic  two  photon  technique. 

From  the  two  photon  echo  experiment  on  the  3S-3D  transition  we  will  be 
able  to  determine  the  relevant  collision  cross-section.  In  our  previous  study 

of  3S-nD  collisional  relaxation  by  the  tri-level  echo  technique  we  went  from 

(2) 

principle  quantum  number  n  =  4  to  n  =  40.  The  tri-level  echo  for  n  =  3  was 
not  accessible  because  the  technique  depends  on  the  energy  separation  of  the 
3S-3P  transition  being  less  than  the  3P-nD  which  is  not  true  for  n  =  3. 

Another  interesting  preliminary  experiment  in  Na  vapor  involves  the 
generation  of  what  we  must  now  call  ordinary  two-pulse  photon  echoes  excited 
on  the  3S-4P  transition.  This  is  the  first  time  two-pulse  ground  state  echoes 
have  been  observed  on  other  than  the  3S-3P  transition.  This  was  accomplished 
by  frequency-doubling  a  laser  excitation  pulse  to  coircide  with  the  3302  X, 

3S-4P  transtiion.  This  transition  is  interesting  because  its  lifetime  is 
about  seven  times  longer  than  the  3S-3P  lifetime.  It  is  therefore  an  interest¬ 
ing  candidate  for  study  of  the  effect  of  velocity  changing  collisions  in  Na 

(3) 

vapor,  to  complement  our  recent  work  in  Li  vapor. 

In  the  process  of  generating  the  3S-4P  echo  we  observed  cascade  superradiance 

from  the  4P  to  4S  to  3P  and  the  4P  to  3D  to  3P  levels .  This  has  also  been  seen 
(4) 

by  other  workers;  what  is  new  is  that  we  also  saw  coherent  radiation  in  the 
3P-3S  transition.  Whereas  the  superradiant  bursts  seen  by  others  could  be  seen 
at  both  the  input  and  output  parts  of  the  Na  sample  cell  the  3P-3S  radiation  was 
observed  only  at  the  output  end.  Here  we  have  a  true  four-photon  process.  Our 
present  set-up  limits  us  to  excitation  pulses  no  shorter  than  6  nsec.  By  shorten¬ 
ing  the  excitation  pulse  to  1  nsec  we  hope  to  observe  the  coherent  3S-3P  radia¬ 
tion  after  the  excitation  pulse  is  gone,  thereby  demonstrating  a  superradiant 
state  phase  transfer. 
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In  another  experimental  set-up  we  have  been  improving  on  the  basic 

echo  measurement  technique.  The  idea  is  to  use  non-parallel  excitation  pulses 

to  permit  non-elec tronic  preliminary  separation  of  the  exciting  beams.  It 

is  well  known  that  an  echo  is  radiated  with  a  polarization  density  whose  k 

vector  is  given  by  k  =  2k2~k^,  where  k^  and  k^  are  the  wavevectors  of  the 

first  and  second  excitation  pulses.  With  high  intensity  lasers  broad  sample 

areas  can  be  irradiated,  so  that  it  is  possible  to  switch  an  unattenuated  echo 

far  out  of  the  beam  divergence  of  initial  laser  pulses.  For  a  beam  of  diameter 

D  the  diffraction  limited  beam  divergence  0  is  given  by  0  ■  1.22  X/D  where  X 

is  the  wavelength  of  the  coherent  light.  For  a  sample  of  length  L,  on  the 

other  hand,  it  is  necessary  that  at  least  AkL  <  tt/4  for  the  polarization  density 

to  contribute  uniformly  throughout  the  sample  without  destructive  interference. 

If  makes  an  angle  0  with  k^  then  as  long  as  0  <<  1  it  is  correct  to  write 
2 

Ak  -  0  k.  Thus,  choosing  AkL  *  ir/4,  we  find 
0/0  -  DA/8r  XL 

which  is  generally  much  larger  than  one. 

Our  present  angled  beam  system  achieves  enough  rejection  of  the  excita¬ 
tion  pulses  to  allow  echo  detection  100  nsec  after  the  first  excitation  pulse. 
For  a  sample  gas  of  lithium  this  corresponds  to  about  four  lifetimes.  Further 
excitation  pulse  rejection  is  achieved  by  using  an  acoustic-optic  modulator 
switch.  Overall  excitation  pulse  rejection  already  exceeds  10^  and  we  expect 
to  be  able  to  improve  that  considerably.  One  last  trick  is  to  apodize  the 
initial  laser  beam,  taking  advantage  of  the  fact  that  in  the  focal  plane  of  a 
lens  placed  in  the  laser  beam  the  light  field  is  the  Fourier  transform  of  the 
light  field  before  the  lens.  Apodizing  in  a  Gaussian  manner  achieves  a 
Gaussian  profile  which  decreases  beam  intensity  outside  the  beam  divergence. 

This  particular  apparatus  will  be  very  valuable  for  all  echo  techniques,  but 
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especially  those  involving  sub-Doppler  linewidth  measurements .  We  already,  for 
example,  are  able  to  "beat"  the  uncertainty  principle  by  observing  echoes  in  Li 
vapor  15  lifetimes  after  the  exicted  state  population  has  died  away  -  the  "life 
after  death"  effect. 

We  have  also  completed  a  theoretical  analysis  of  the  effects  of  velocity 
changing  collisions  in  optical  experiments.^  To  further  aid  in  our  study 
of  velocity  changing  collisions  we  have  set-up  a  cw  dye  laser  to  burn  holes  in 
the  Doppler  line.  With  this  we  will  be  able  to  modify  the  effective  Doppler 
profile  and  thus  achieve  more  experimental  control. 

Further  measurements  and  analysis  of  photon  echo  modulation  in  the 
1  3  3+ 

i>2  “  optical  transition  in  Pr  :LaF^  have  resulted  in  a  more  precise 
evaluation  of  the  spectroscopic  and  relaxation  characteristics  of  this  system. 
Fourier  transformation  of  the  modulation  yields  the  hyperfine  splitting  fre¬ 
quencies  and  the  related  linewidths. 

3+ 

The  resonance  widths  of  Pr  iLaF^  reflect  the  total  (inhomogeneous  and 
homogeneous)  nuclear  linewidths  of  the  various  hyperfine  transitions.  Pre¬ 
vious  measurements^  ’  ^  have  shown  the  ground  state  widths  to  be  about  200  kHz. 
This  width  is  mostly  inhomogeneous  and  arises  from  the  interaction  of  the  en- 
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hanced  ground  state  Pr  nuclear  moment  of  11.5  kHz/G  with  the  surrounding 
distribution  of  static  ^F  nuclear  moments.  Macfarlane  and  Shelby ^  have 
measured  the  g-tensor  for  the  state  and  found  the  enhancement  to  be  roughly 
one-fifth  as  much  as  in  the  ground  state.  This  would  lead  one  to  expect  a  much 
narrower  linewidth  for  the  nuclear  transitions  in  the  state  and  our  measure¬ 
ment,  while  very  nearly  resolution-limited  by  the  overall  echo  decay  envelope, 
establishes  a  new  upper  limit  of  30  kHz  [FWHM]  for  the  hyperfine  line  at  8.51  Mhz, 
consistent  with  this  expectation.  The  linewidths  of  the  hyperfine  at  3.72  and 
4.79  Mhz  are  somewhat  larger,  being  70  kHz  and  60  kHz  respectively.  The  results 
for  the  ^2  excited  state  parameters  are  summarized  in  Table  1. 
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Table  1 


Full  width 

Resonance 

Frequency 

at  half  max 

e 

U1 

3.72  +  0.02  MHz 

70  +  20  kHz 

e 

*2 

4.79  +  0.02  MHz 

60  +  20  kHz 

e 

w12 

8.51  +  0.02  MHz 

30  ±  30  kHz 

able  1:  Values  for  the  ^D2  excited  state  hyperfine  frequencies 
id  nuclear  linewidths  obtained  from  the  Fourier  transform  of  the 
ita.  The  linewidths  have  been  corrected  for  the  effect 


f  the  finite  truncation  interval 


The  effective  hyperfine  Hamiltonian  is  given  by 


Hg(e)  "  Pg(e)  ^z(Z) 


qg(e) 

3 


(I 


x(X) 


-  L  1 
y(Y)J 


(l) 


inhere  P  ,  .  and  n  .  are  interaction  constants.  The  upper  and  lower  case 
g(e)  g(e; 

subscripts  emphasize  that  the  ground  (g)  and  excited  (e)  state  principal  axes 
need  not  be  the  same,  although  Neumann's  principle  requires  that  they  have 
at  least  one  common  axis  parallel  to  the  local  symmetry  axis.  This  constraint 
significantly  reduces  the  number  of  possible  orientations  and  we  may  represent 
these  possibilities  on  a  two-dimensional  map/^  The  coordinates  of  such  a 
map  give  the  particular  choice  of  common  axis  ,  and  the  angle  6  through 
which  one  set  of  axes  is  rotated  with  respect  to  the  other  set  about  the  common 
axis.  Since  the  Hamiltonian  in  Eq.  1  is  insensitive  to  the  sign  of  the  co¬ 
ordinates  an  orientation  of,  say,  x  | |  X  is  equivalent  to  x  | |  -X.  We  have 
calculated  the  echo  modulation  in  increments  of  five  degrees  in  the  rotation 
angle  9.  In  the  map  presented  in  Pig.  3,  for  each  orientation  coordinate  a 
box  is  drawn  whose  width  is  proportional  to  the  minimum  mean  square  deviation 
between  our  data  and  the  calculated  modulation.  Only  the  first  microsecond  of 
data  were  used  in  the  fitting  procedure  for  Figure  3>  Beyond  one  usee,  the 
data  are  dominated  by  the  excited  state  hyperfine  splittings, due  to  the  relative 
linewldths  of  ground  and  excited  state  levels, and  Including  these  data  in  the 
fit  considerably  reduces  the  selectivity  of  the  map. 

Figure  4  shows  in  expanded  scale  the  first  three  ysec  of  the  data, 
and  Figure  5  shows  the  theoretical  curve  that  was  calculated  using  the  orienta¬ 
tion  coordinate  giving  the  best  fit  to  the  data  (as  indicated  by  an  arrow  in 
the  lower  left  of  Figure  3) . 

The  map  in  Figure  3  shows  a  high  degree  of  symmetry  and  concommitant 
lack  of  selectivity  for  determining  the  orientation  parameter.  To  explore  the 
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Figure  3. 
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extent  to  which  this  symmetry  is  intrinsic,  we  have  substituted  for  the  data 

used  to  calculate  the  map  of  Figure  3  the  values  of  the  theoretical  modulation 

corresponding  to  Figure  5  ,  producing  the  map  shown  in  Figure  6.  The  scale 

in  Figure  6  was  chosen  to  make  the  box  corresponding  to  the  coordinate  with 

the  worst  fit  in  Figure  6  approximately  the  same  size  as  the  box  with  the 

worst  fit  in  Figure  4.  It  is  apparent  that  the  high  degree  of  symmetry  per- 

sists.  To  estimate  the  extent  to  which  our  signal-to-noise  would  have  to 

improve  to  allow  us  to  select  one  orientation  with  certainty,  we  replotted 

Figure  6  with  the  box  width  scale  expanded  by  a  factor  of  ten,  as  shown  in 

Figure  7,  It  appears  that  an  improvement- in  signal-to-noise  by  a  factor  of 

ten  or  twenty  would  have  made  possible  a  more  definitive  statement  about  the 

orientation.  This  improvement  will  probably  come  from  the  use  of  frequency 

and  amplitude  stabilized  dye  lasers.  Frequency  jitter  in  particular  causes 

considerable  fluctuation  in  echo  intensity,  presently  averaged  by  the  data- 

taking  system, and  such  averaging  tends  to  wash  out  some  of  the  fine  details 

of  the  modulation  pattern.  Despite  the  high  symmetry,  we  can  draw  a  number  of 

conclusions  about  the  relative  orientation.  All  configurations  with  all  axes 

parallel  are  ruled  out,  as  are  those  with  the  excited  state  Z-axis  parallel  to 

C£.  These  conclusions  are  similar  to  the  measurement  of  Chen  et  al./^ 

3  3 

on  the  relative  orientation  of  the  and  Pq  states.  Since  the  g-tensor  is 
nearly  isotropic  in  the  state^  and  absent  completely  in  the  ^Pq  state 
the  relative  orientation  parameter  must  be  primarily  determined  by  the  rela¬ 
tive  orientation  of  the  pure  quadrupole  tensor  with  respect  to  the  ground  state 
hyperflne  tensor, 

*Thi8  research  was  also  supported  by  the  National  Science  Foundation  and  the 
Army  Research  Office  under  Grant  NSF-DMR80-06966  and  the  U.  S.  Office  of 
Naval  Research  under  Contract  No.  N00014-78-C-0517. 
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Figure  6.  Map  of  mean  square  deviation  between  the  theoretically  calculated 
data  of  Figure  S  and  the  theoretical  curve  of  Equation  1. 

Notice  the  similarity  to  Figure  3. 
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SIGNIFICANT  ACCOMPLISHMENTS  AND  TECHNOLOGY  TRANSITION  REPORT 

* 

Significant  Accomplishments 

I.  Quantum  Detection  and  Sensing  of  Radiation 

Luminescence  radiation  has  been  characterized,  for  the  first  time, 
in  terms  of  the  underlying  physical  processes  giving  rise  to  it.  This  formula* 
tion  shows  how  differently  luminescence  radiation  behaves  from  thermal  light. 

A  new  method  has  been  developed  for  incorporating  dynamics  into  particle, 
photon,  and  pulse  multiplication  processes. 

A  new  technique  for  characterizing  fiber-optic  system  performance  and 
for  designing  single-threshold  fiber-optic  receivers  has  been  developed. 

II.  Physical  Properties  and  Effects  in  Electronic  Materials 

A  technique  for  enhancing  photochemical  reactions  through  use  of 
photolithographically  prepared  substrates  has  been  developed.  This  approach 
has  potential  application  to  the  development  of  new  chemical  processes  through 
the  use  of  catalyzed  photochemical  reactions.  In  addition,  such  surface  chemical 
reactions  can  be  used  to  produce  novel  microstructures  for  microelectronics 
devices . 

A  physical  theory  has  been  developed  using  a  charge  scattering  model  to 
describe  the  grain  boundary  transport  in  polycrystalline  semiconductors.  The 
calculated  result  explains  the  reduction  of  majority  carrier  current  in  agree¬ 
ment  with  experimental  data.  It  provides  a  theoretical  basis  for  use  of  poly¬ 
crystalline  semiconductors  in  VLSI  interconnects  and  solar  cells. 
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III.  Energy  Transfer  and  Relaxation  in  Small  Polyatomic  Molecules 

Relaxation  processes  involving  the  electronically  excited  states  of 
the  1^  molecule  have  been  discovered  which  may  shed  light  on  the  operation 
of  the  iodine  chemical  laser. 

A  laser  photofragmentation  technique  has  been  applied  to  the  study 
of  fast  hydrogen  atom  interactions  with  molecules.  These  studies  are  of 
potential  importance  in  plasma  research,  chemical  reaction  diagnostics,  and 
infrared  laser  development. 


IV.  Picosecond  Energy  Transfer  and  Photofragmentation  Spectroscopy 

We  have  discovered  that  one  of  the  paths  in  the  photodissociation  of 
an  anthracene  endoperoxide  produces  not  only  ^0^  but  also  an  excited  aromatic 
fragment.  This  represents  one  of  the  few  cases  of  a  photodissociation 
producing  both  fragments  in  excited  electronic  states  and  furthermore  is 
the  first  time  that  the  state  of  resulting  from  a  photodissociation  has 
been  established.  The  relative  importance  of  competing  dissociative  channels 


for  generating  0?  is  determined  and  the  role  of  a  short-lived  (45  ps)  chemical 
intermediate  in  process  is  proposed. 


V.  Generation  and  Control  of  Radiation 


A  simple  theoretical  model  explaining  photon  echo  formation  is  Gaseous 
media  has  been  developed.  It  is  called  the  Billiard  Ball  Echo  Model  and 
works  equally  well  with  two-pulse,  three-pulse,  two-photon,  Raman,  Grating, 
standing-wave,  etc.,  photon  echoes.  The  key  is  atomic  recoil.  Echo  genera¬ 
tion  can  now  be  understood  without  recourse  to  lengthy  mathematical  analysis 
thereby  facilitating  the  use  of  echoes  to  solve  problems  in  science. 


Technology  Transition  Reports 

1.  A  new  nonlinear  signal  processing  scheme  to  selectively  provide  improved 
signal-to-noise  ratios  for  systems  including  multiplication  noise  has  been 
conceived.  This  may  be  used  as  part  of  the  optical  guidance  system  in  the 
NASA/JPL  Gallileo  orbital  mission  scheduled  to  orbit  Jupiter  in  1988. 

2.  A  frequency-shift  keyed  fiber  optic  data  link  that  may  be  used  in  con¬ 
junction  with  wavelength  division  multiplexing  for  high-capacity  multi-user 
data  networks  has  been  demonstrated. 


3.  Although  too  numerous  to  mention,  there  are  many  examples  of  laser 
applications  to  defense,  industrial,  and  research  uses.  One  of  the  most 
potentially  interesting  is  the  use  of  lasers  in  photochemical  direct  writing 
of  microcircuits.  The  original  research  leading  to  the  development  of  the 
laser  was  begun  under  the  JSEP  program  at  Columbia  in  1951  when  work  on  the 
maser  started  under  the  direction  of  Professor  Townes. 
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Naval  Research  Laboratory 

AITN:  Mr.  B.  D.  McCombe,  Code  6800 

Washington,  DC  20375 

Commanding  Officer 
Naval  Research  Laboratory 
ATTN:  Mr.  W.  L.  Faust,  Code  6504 
Washington,  DC  20375 

Technical  Director 

Naval  Underwater  Systems  Center 

New  London,  CT  06320 

Naval  Research  Laboratory 
Underwater  Sound  Reference  Detachment 
Technical  Library 
P.  0.  Box  8337 
Orlando,  FL  32856 

Naval  Ocean  Systems  Center 
Aim  Dr.  P.  C.  Fletcher,  Code  92 
4 an  Diego.  CA  92152 

Naval  Ocean  Systeas  Center 
ATTN:  Mr.  W.  J.  Dejka,  Code  8302 
San  Diego,  CA  92152 

Naval  Ocean  Systeas  Center 
ATTN:  Dr.  Alfred  K.  Nedoluha, 

*  Code  922 

San  Diego,  CA  92152 

Naval  Weapons  Center 

ATTN:  Dr.  G.  H.  Winkler,  Code  381 

China  Lake,  CA  93555 

Dr.  Donald  E.  Kirk  (62) 

Professor  and  Chairman,  Electrical 
Engineering 
SP-304 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Dr.  D.  F.  Dence 

Naval  Underwater  Systeas  Center 
New  London  Laboratory 
ATTN:  Code  34 
New  London,  CT  06320 

Director,  Technology  Assessment 
Division  (OP-987) 

Office  of  the  Chief  of  Naval  Oper. 
Navy  Department 
Washington,  DC  20350 

Mr.  J.  W.  Willis 

Naval  Air  Systems  Coomand 

AIR-310 

Washington,  DC  20361 

Naval  Electronics  Syscems  Command 
NC  n 

ATTN:  Code  61R 

2511  Jefferson  Davis  Highway 

Arlington,  VA  20360 

Department  of  the  Navy 
Naval  Sea  Syscems  Consand 
ATTN:  W.  W.  Blaine  ( SEA-6 2R) 
Washington,  DC  20362 


David  Taylor  Naval  Ship  Research 
and  Development  Center 
ATTN:  Mr.  G.  H.  Gleissner,  Code  18 
Bethesda,  MD  20084 

Mr.  Martin  Mandelberg 
Coast  Guard  RAD  Center 
Avery  Point 
Groton,  CT  06340 

Naval  Underwater  Systeas  Center 

New  London  Laboratory 

ATTN:  10 IE  (Dr.  Edward  S.  Eby) 

New  London,  CT  06320 

Mr.  Thomas  J.  Manuccla,  Head 
Chemistry  and  Application  Section 
Code  6543 

Naval  Research  Laboratory 
Washington,  DC  20375 

Dr.  Stephen  G.  Bishop,  Head 
Semiconductor  Branch 
Code  6870 

Naval  Research  Laboratory 
Washington,  DC  20375 

Dr.  William  F.  Gabriel 
Antenna  Systems  Staff 
Code  5342 

Naval  Research  Laboratory 
Washington,  DC  20375 

Dr.  John  W.  Rockwav 
Comm.  Technology  Prog.  Off. 

Code  8105 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Dr.  Barry  P.  Shay 
Systeas  Integration  and 
Instrumentation  Branch 
Code  7522 

Naval  Research  Laboratory 
Washington,  DC  20375 

Dr.  Sydney  R.  Parker 
Professor,  Electrical  Engineering 
Code  62PX 

v*val  Postgraduate  School 
Moncerey ,  CA  93940 

Dr.  George  B.  Wright 
Office  of  Naval  Research 
Code  427 

Arlington,  VA  22217 

OTHER  GOVERNMENT  AGENCIES 

Dr.  Ronald  E.  Kagarlse 
Director 

Division  of  Materials  Research 
National  Science  Foundation 
1800  G  Street 
Washington.  DC  20550 

Director 

Division  of  Electrical,  Computer 
and  Systems  Engineering 
National  Science  Foundation 
Washington,  DC  20550 


Dr.  Dean  L.  Mitchell 
Section  Head 

Condensed  Matter  Sciences  Section 
Division  of  Materials  Research 
National  Science  Foundation 
1800  G  Street,  N.  W. 

Washington,  DC  20550 

Judson  C.  French,  Director 
Center  for  Electronics  and  Electrical 
Engineering 

3  358  Metrology  Building 
National  Bureau  of  Standards 
Washington,  DC  20234 

SON-GOVERNMENT  AGENCIES 

Director 

Columbia  Radiation  Laboratory 
Columbia  University 
538  West  120th  Street 
New  York,  NY  L0027 

Director 

Coordinated  Science  Laboratory 
University  of  Illinois 
Urbana,  IL  61801 

Associate  Director  pf  Materials 
and  Electronics  Research 
Division  of  Applied  Sciences 
McKav  Laboratory  107 
Harvard  University 
Cambridge.  MA  02138 

Director 

Electronics  Research  Center 
University  of  Texas 
P.  0.  Box  7728 
Austin.  TX  78712 

Director 

Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 

Director 

Electronics  Sciences  Laboratory 
University  of  Southern  California 
L09  Angeles,  CA  90007 

Director 

Microwave  Research  Institute 
Polytechnic  Institute  of  New  York 
333  Jay  Street 
Brooklyn,  NY  11201 

Director 

Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 

Director 

Stanford  Electronics  Laboratory 
Stanford  University 
Stanford,  CA  94305 

Director 

Edward  L.  Glnzton  Laboratory 
Stanford  University 
Stanford,  CA  94305 
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Or.  Outer  Cutiu 

School  of  Electrical  Engineering 

Cornell  University 

316  Phillip*  Hell 

t theca,  ST  14850 

Or.  Carlton  Halter 
Electro  Science  Laboratory 
The  Ohio  State  Univaraity 
1320  Kinnear  Road 
Coluabua.  08  43212 

Or.  Richard  Saaka 
Dept,  of  Eleccrical  Engineering 
Texas  Tech  Unlvarslcy 
Lubbock,  TX  79409 

Director 

School  or  Eleccrical  Eaginuring 
Georgia  Inatltute  of  Technology 
Atlanta,  GA  30332 

Or.  John  F.  Walkup 
Dept,  of  Eleccrical  Engineering 
Texaa  Tech  Univaraity 
Lubbock.  TX  79409 

Mra.  Renata  D 'Arcangelo 

Editorial  Office 

130  Fierce  Ball 

Olvlaion  of  Applied  Sclencea 

31  Oxford  Screec 

Cambridge,  MA  02138 
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